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, ABSTRACT 

>«✓ We present deep 450/xm and 850/xm observations of a large, uniformly covered 

394 arcmin^ area in the COSMOS field obtained with the ScuBA-2 instrument on the 
James Clerk Maxwell Telescope (JCMT). We achieve root-mean-square noise values of 
cr45o=4.13 mJy and iT850 =0-80 mJy. The differential and cumulative number counts are 
presented and compared to similar previous works; this region of COSMOS is found 
to have a counts deficit at both wavelengths. Individual point sources are identified at 
^ >3.6o' significance at 450/xm and 850/im (a threshold corresponding to a 3-5% sample 

—i contamination rate). We identify 78 point sources at 450/im and 99 at 850/um, rang- 

ing in flux densities from ^450 =13-37 mJy and ^sso =2-16 mJy. Approximately half 
of both samples are detected at both wavelengths. Since the positional uncertainties 
at 450/im are small (1-2.5"), this allows a precise identification of multiwavelength 
counterparts without reliance on detection at 24/im or radio wavelengths; we find 
that only 40% of 450/xm-sources have 24/im or radio counterparts. Other submillime- 
^J!^ ter work reliant on such counterparts could be significantly biased. We use the vast 

ancillary data available in the COSMOS field to characterise the ScuBA-2-selected 
galaxies using optical/near- infrared photometric redshifts, and we find redshift distri- 
butions that peak at {z) = 2.0±1.1 for 450/im-selected galaxies and (z) = 2.2 ±1.6 for 
850/im-selected galaxies. 850/im-identified galaxies and 450/im-identified galaxies oc- 
". • cupy similar parameter space in redshift, luminosity, dust temperature and mass with 

^ surprisingly few population differences; the primary difference between the two is their 

k>( far-infrared S850/S450 colour. This colour does not evolve with redshift. The 450/im 

sub-sample which is not detected at 850/Ltm are not warmer-dust systems, as previously 
suggested in the literature, but instead they are lower luminosity sources which just 
fall under the detection limit at 850/Ltm. All Scuba-2 sources have UV/optical/near- 
infrared derived star formation rates which underpredict their infrared-based SFR by 
factors ^10; this has implications of how and if infrared-luminous systems sitting 
along the "main sequence" can be interpreted. The total contribution of these submil- 
limeter sources (both 450/xm and 850/tm) to the cosmic star formation rate density 
is measured as roughly flat at wO.Ol M© yr~^ Mpc~^ across < 2; < 5; we end by 
discussing implications for future, deeper observations of short-wavelength single-dish 
bolometers. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: infrared — 
galaxies: starbursts 
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1 INTRODUCTION 



The arrival of the Submilhmeter Common User Bolometer 



Array (ScuBA; Holland et al. 19991 on the James Clerk 



Maxwell Telescope (JCMT) fifteen years ago ushered in a 
new age of galaxy evolution studies by highlighting the im- 
portance of distant infrared-bright, ultraluminous galaxies 
to the buildup of the Universe's stellar mass and the forma- 
tion of massive, local elliptical galaxies. Submilhmeter galax- 
ies (SMGs, canonically selected at 850/im with Ssso ~ 2- 
5mJy;|Smail et al |[T997| [Barger et aL][T998l [Hughes et ah] 
1998[ [Eales et al. 1999 1 have since been shown to have a 
peak volume density at z « 2.2 — 2.5 ( Chapman et al.|2003 



2005 Yun et al. 20121. Follow-up studies of their molecu- 



lar and ionised gas properties (|Swinbank et al. |2004 


Neri 


et al.ll2003l |Greve et al. 2005| |Tacconi et al. 2006| 


2008 


BothwoU ct al.||2010| Engel et al.||2010| Banerji et al.| 


2011 


Alaghband-Zadeh et al.|2012 1 indicate that most SMGs are 



triggered by major mergers of gas-rich disk galaxies, similar 
to local (ultra-) luminous infrared galaxies, (U)LIRGs (e.g. 

[Sanders fc Mirabel|[T996l [Armus et aL] 



Sanders et al.||1988 



2009 U et al.[[2012 K Furthermore, several detailed studies 
of the X-ray and mid-infrared properties of SMGs have led 
to the conclusion that the majority of the SMGs' bolometric 



energy is star formation dominated (Alexander et al. 2005 



Pope et al.pOOSl [Coppin et al.[[2008[ [Menendez-Delmestre[ 
et al.|2009||Coppin et al.|2010[ [Laird et al.[[2010[ ). 

The limiting factors in the analysis and characterisa- 
tion of the SMGs is small number statistics and the lack 
of secure counterpart identifications. The vast majority of 
the studies mentioned above— although thorough— are based 
on small samples of SMGs, anywhere from 2-75 galaxies 
(but mostly on samples with less than ten sources). Some 
of these SMGs could have been mischaracterised due to 
incorrect counterpart identifications, which may be a con- 
sequence of the large beamsize of submm single-dish ob- 
servations (the JCMT beamsize at 850/im is ~15"). The 
problem of limited statistics can be alleviated by wide-field 
submillimetre mapping. Mapping large regions of sky at 
long wavelengths is technically challenging, particularly if 
the maps are intended to be sufficiently deep and sensitive 
enough to detect unlensed galaxies in the early Universe. 
New submillimeter/far-infrared (FIR) observatories, such as 
AZTEC ( [Wilson et aL]|2008| l on the L arge MiUimeter Tele - 
scope and Herschel Space Observatory ( [Pilbratt et al.|2010[ ), 
have mapped larger sky areas in recent years. Herschel has 
mapped hundreds of square degrees at 250-500/:im. However, 
such surveys are more limited in sensitivity than SCUBA, 
primarily due to the smaller primary dish, which results in 
large beam sizes and much higher confusion limits. 

The new Submilhmeter Common User Bolometric 
Array-2 (ScuBA-2) instrument on JCMT presents exciting 
new avenues to explore infrared-bright starbursts in the 
early Universe. Besides the four-fold increase in bolome- 
ter arrays and large field-of-view scans, which make SbOfim 
mapping more efficient than it was with ScuBA, the SCUBA- 
2 450/im bolometers provide the first opportunity to identify 
450/im-bright sources with remarkable resolution in the FIR; 
the JCMT beamsize at 450/^m is ~7", compared to ~36" 
at 500/im on Herschel. Direct detection in the FIR with a 
small beamsize is a large step forward, since multiwavelength 



counterpart identification has been a key limiting factor in 
the interpretation of submillimeter-bright sources over the 
past decade. Deep ScuBA-2 maps, which are much less lim- 
ited by confusion noise, are the scientific complement to the 
large-scale mapping done with Herschel. 

This paper presents 450/im and 850/xm observations 
from ScuBA-2 in a central region of the Cosmic Evolu- 
tion Survey (COSMOS) field fSco ville et al.|2007i ) surround- 
ing the Cosmic Assembly Near-infrared Deep Extragalac- 



tic Legacy Survey (CANDELS) area. Recently, Geach et al 



(20121 presented deep ScuBA-2 450/im observations in the 
central ~100''^ of this area. Our observations have a wider. 



more uniform coverage than the data presented in [Geach[ 
et al. (20121, with the motivation of selecting a homoge- 
neous sample, while the [Geach et al.[ ( |2012[ ) data are deeper 
and patterned to find fainter sources. The pointing on the 
CANDELS area is motivated by the perceived high correla- 
tion of submillimeter emission with near-infrared emission 
(e.g. [Wang et aL][2006[ [Serjeant et al.[[2008[) and the abil- 



ity to morphologically characterise SMG counterparts using 
the Hubble Space Telescope's high spatial resolution. How- 
ever, we leave the detailed morphological analysis of these 
galaxies to a future work. This work is also complemented 
by the deeper, yet smaller area ScuBA-2 coverage and anal- 



ysis in the lensed cluster field A 370 by Chen et al. (20121, 
which currently provides the deepest 450/im counts. 

Section |2| describes our ScuBA-2 observations and data 
reduction as well as ancillary data used, section [3] discusses 
the identification and flux density measurements of point 
sources in the maps and differential and cumulative num- 
ber counts, section [4] describes the multiwavelength proper- 
ties of the ScuBA-2-selected galaxies, and section|5]discusses 
their implications, relation to similar starburst galaxy pop- 
ulations, and implications for future, deeper ScuBA-2 obser- 
vations. Throughout we assume a A CDM cosmology wit h 
//o=71kms-^Mpc-^ and n„,^0.27 ( [Hinshaw et al.f2009 l. 



2 OBSERVATIONS & DATA 

Observations were taken with the ScuBA-2 instrument on 
the JCMT on 2011-Dec-26, 2011-Dec-28, 2012-Feb-04, 2012- 
Feb-05, 2012-Feb-07, 2012-Apr-09, 2012-Apr-30, 2012-May- 
01, 2012-May-02, and 2012-Dec-21. Conditions were opti- 
mum, with the optical depth varying from 0.02< r225GHz < 
0.05 and averaging r225GHz ~ 0.04. At the wavelengths 
of our observations, the optical depths were approximately 
T450fim ~ 0.55 and rsso^ ~ 0.16. The integration times on- 
field per night were 3.5 hrs, 3.6 hrs, 3.1 hrs, 5.5 hrs, 5.5 hrs, 
4.5 hrs, 2.9 hrs, 3.8 hrs, 0.6 hrs and 5.0 hrs respectively, to- 
talling 38.0 hrs on-field. We centred our map at position 
10:00:28.0, -1-02:24:00. The pong-900 mapping pattern was 
used to achieve a uniform RMS over a large area «15'xl5' 
in preference to the Daisy mapping which has substantial 
patterning and non- uniformity and covers a smaller area. 

Raw data were downloaded via the Canadian Astron- 
omy Data Centre (CADC) JCMT Science Archive and pro- 
cessed with the ScuBA-2 software packages Smurf, the Sub- 
Millimetre User Reduction Facility, and Picard, Pipeline 
for Combining and Analyzing Reduced Data. Data were 
split up by night, wavelength, and sub-array for process- 
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Arc Minutes 

Center: R.A. 10 00 24.08 Dec +02 24 04.8 

Figure 1. The COSMOS AbOfim ScuBA-2 map. The colour scale is linear with upper and lower ranges of -15< 5450 <25 mjy/beam. 
Contours mark the 2, 3, 4, 5, and 6x the central RMS value, measured to be 4.13mjy/beam at 450/im. We search for sources which are 
detected with RMS less than double the central value (roughly corresponding to the inner most contour). The area probed by this cut 
has a radius '~11.2' and an area ~394'^. The 450^m-identified >3.6(7 sources are identified with white boxes (a justification of the 3. 60" 
cut can be found in the text). 



ing. Smurf processing was first executed on each night's 
sub-array data (whicli is equivalent to combining each sub- 
array per Minimum Scheduable Block (MSB), and then com- 
bining all MSBs). Partial maps were constructed using the 
MAKEMAP routine in Smurf with the default deep extra- 
galactic configuration file optimising background subtrac- 
tion. The MAKEMAP uncalibrated maps were combined at 
each stage (by sub-array, then nights) using the Picard 
recipe mosaic_,tcmt_images. Each of the four sub-arrays 
at 450/^m and 850/im were combined for each night's worth 
of data, and then the nights were combined to produce the 
unprocessed, raw maps at 450/xm and 850/im in units of pW. 

We applied the matched-filter Picard recipe to sup- 
press pattern noise and increase the signal from individ- 



ual sources, which are generally unresolved. At 450/im, 
the map RMS noise drops substantially from 14.08 mjy to 
4.13 mjy. A similar yet less dramatic drop is seen at 850/im, 
from 2.18 mJy to 0.80 mJy. Flux calibration is done by ap- 
plying the measured flux conversion factors (FCF) from 
the first eight months of calibrator data, primarily from 
Uranus and Mars (FCF85o=556±45 Jy pW"^ beam"^ and 
FCF45o=606±55 Jy pW~^ beam~^). Calibration data taken 
during the observations for this project agree with these 
FCF values to ~10%. Final reduced maps, in units of 
mJy/beam, are shown in Figures [1] and [2] The units are 
such that the flux density, in mJy, for an unresolved source 
is simply the sources' peak value in the map. Signal-to-noise 
maps are constructed using the Picard recipe makesnr; 
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Figure 2. The COSMOS 850^m Scuba-2 map, plotted in the same fashion as Figure ^ The colour scale is linear and runs from 
-5< 5850 <10mjy/beam, and the measured central RMS is 0.80 mjy/beam. All >3.6a- sources are marked with white boxes. Some 
sources appear outside the nominal 2xRMS contour, but locally have noise that dips below the 2xRMS threshold. 



these are not shown here but are quite similar to the final 
reduced maps within the central area. Signal-to-noise maps 
are used for the identification of point sources. 

To avoid contamination from the less sensitive edges of 
the map in calculating an overall characteristic map noise, 
we first measured the effective area of the sensitive region 
of the map, where the RMS was less than or equal to a 
factor of two times the RMS in the centre. We determined 
this boundary by measuring the RMS in concentric rings 5" 
wide from a jackknife map, having subtracted out the signal 
from point sources. The jackknife map represents the instru- 
mental noise of our observations with sources removed and 
is constructed by subtracting one half of our data from the 
other, then scaling the noise by the square root of the in- 
tegration time. This 2xRMS boundary lies ~11.2' from the 
map centre in both the 450pm and 850/im maps, correspond- 



ing to an area of 394'^, slightly larger than the anticipated 
15'xl5' PONG-900 target area. The pixel fiux distributions 
for the ~400'^ central region are shown in Figurejsjwith the 
distributions of the jackknife pure-noise maps, in the same 
~400'^ area, highlighted in gray. 



2.1 COSMOS Ancillary Data 

We use the vast ancillary data available in the COSMOS 



field ( Scoville et al^^2007 ) to investigate the physical na- 



ture of the ScuBA-2 submillimeter-bright sources and com- 
pare with larger samples of submillimeter sources previ- 
ously analysed in the literature. Photometric redshifts were 
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Figure 3. Pixel flux distributions (per unit area) within the central ~400'^ area at 450^m and 850/im. The solid black lines indicate 
the distributions of the data, showing excesses at high flux densities where point sources contribute to the distributions. At lower flux 
densities, the distributions are dominated by instrument noise. The negative excesses are due to residual troughs around bright sources (a 
characteristic of the matched filter technique). The gray shaded histograms show the pixel noise distributions within the same region for 
the jackknife maps. As discussed in the text, the jackknife maps were constructed by differencing two halves of the data, which removes 
point sources and leaves only residual noise. 




345678 345678 
Measured S/N ot 450/.im Measured S/N at 850ju.m 



Figure 4. Our results from Monte Carlo simulations, showing more than 10000 source realisations and the measured boosting factor 
and positional offset for each source (gray points) at 450/im (left) and 850/im (right). The boosting factor is measured as the ratio of 
output measured flux density over input flux density. We deboost our sources' flux densities using the median boosting factors (red lines). 
The 90% confidence intervals (dashed blue lines) and 95% confidence intervals (dotted blue lines) are overplotted. We assign a positional 
uncertainty to sources according to the 90% confidence interval. The 3.6(t signal-to-noise threshold is marked by dotted vertical lines. 



fit using the Le Phare^ code (Ilbert et al. 



2009 



20101 



to multi-band ultraviolet, optical, near-infrared (NIR) and 
mid-infrared (MIR) pliotometry. Stellar population tem- 



^ http:/ /www. cfht.hawaii.edu/~arnouts/LEPHARE/cfht_lephare/ 
lephare.html 



plates ( [Polletta et al.||2007| [Bruzual fc Charlot|[2003l ) were 
assumed as input to Le Phare. In this paper we use the 
measured output parameters of stellar mass and star forma- 
tion rate. Both are measured from stellar template match- 
ing to UV/optical/near-infrared photometry only, and in 
this paper, we refer to the template star formation rate as 
"SFRuv." The imaging and associated data reduction for 
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Spitzer Space Telescope IRAC and MIPS data is described 
in|Sanders et al.| ([200YI), |Le Floc'h et al.|([2009t and|Frayer 



et al. ( 2009 1 



3 IDENTIFICATION OF POINT SOURCES 

Point sources are extracted in eacli map by isolating high 
signal-to-noise pixel groups (e.g. adjacent pixels) in the 
signal-to-noise maps. Within each pixel group, the highest 
signal-to-noise pixel marks the object centre and point at 
which the raw flux density is measured. The raw flux den- 
sities are measured as the peak flux density of the source 
in mjy/beam, assuming the sources are unresolved (a safe 
assumption for high-redshift infrared galaxies). 

Here we describe our use of Monte Carlo simulations 
to determine an appropriate detection threshold for report- 
ing detections, and for measuring immber counts, boosting 
factors and completeness and contamination rates. 



3.1 Monte Carlo Simulations 

The estimation of cumulative number counts at both 450/im 
and 850^m, along with the completeness and contamination 
rates of our samples, deboosting factors, and adequate detec- 
tion thresholds requires the use of Monte Carlo simulations. 
This is done by injecting fake sources into a pure-noise, or 
jackknife, map. Injecting sources into our map with known 
flux densities and positions, we can measure the accuracy by 
which we recover those sources using the above extraction 
method. 

The injection of fake sources into noise maps requires 
an assumption of the intrinsic surface density of submillime- 
ter sources and their intrinsic flux densities. First, we inject 
our jackknife map with sources following the prescriptions of 
|Chen et al!] ( |2012| ) and |Geach et al.| ( |2012| at 450/xm. How- 
ever, injecting with their flt parameters does not predict our 
measured raw number counts optimally. Keeping the mod- 
els the same, we vary the flt parameters from [Chen et al.| 
(20121 and Geach et al. (20121, which use a double-power 



law and Schechter form respectively, until the raw number 
counts extracted from our Monte Carlo simulations accu- 
rately reflects the extracted raw number counts from our 
real data. The optimal parameters we measure are given in 
Table |2] and described later in the text. At 850/im we run a 



similar test but use Knudsen et al. ( 2008 1 and Chen et al 
( 2012[ ) formulations as initial input. In all cases, sources are 



injected down to and beyond the flux densities anticipated 
to correspond with 100% of the far-infrared Extragalactic 
Background Light (EBL), which are at Ril mjy at 450/im 



and ~0.5mJy at 850^m (estimates from 


Fixsen et al. 


1998 


Barger et al.T999!|Coppin et al.|2006||Knudsen et al. 


2008 


Oliver et al. 2010,, Chen et al.|2012 


1. To generate substantia] 



statistical samples of input and extracted sources, we gener- 
ate 200 maps at both wavelengths (which generates ~10,000 
extracted sources). Boosting, positional offset, completeness 
and contamination are all estimated using the known input 
sample, position and flux density, and the measured output 
sample, position, flux density and signal-to-noise. 

The flux densities of sources we measure from our maps 
need correction for flux boosting. Sources' flux densities are 
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Figure 5. The completeness curves at 450/tm and 850/im. Com- 
pleteness is measured as a function of input flux density as tlie 
number of sources in a given flux bin which are recovered at the 
signal-to-noise threshold (3fT and 4cr are plotted here). The un- 
certainty on the completeness curve is bimodal. 



expected to be boosted in two different ways. First, there is 
statistical variation around sources' true flux densities. More 
sources with intrinsically fainter flux densities means more 
sources are likely to 'scatter up' and contaminate higher 
flux density bins where there are intrinsically fewer sources. 
This Eddington boosting (described by (Eddington|19'T3 \ as- 
sumes that fainter sources are more numerous, but overall, 
it is largely insensitive to the number counts distribution 
assumed at flux densities much lower than the detection 
threshold. The second form of flux boosting comes from con- 
fusion. This relates to the beamsize of our observations and 
the sky density of low S/N sources, whereby multiple sources 
fall into one beamsize and contribute flux to a measured 
source. We can estimate the magnitude of these two boost- 
ing factors together using our simulations results, which we 
show as a function of detection S/N in the top panels of 
Figure [1] 

Figure|4]also shows the median positional offset between 
the measured output position of a 450/tm- or 850/im- source 
and its input position in the bottom panels. The positional 
uncertainty quoted in Tables [l] and [5] is the 90% confldence 
interval, which we then use for matching to radio, 24/im, 
and optical counterparts as described in sections |3.4| Higher 
signal-to-noise detections can yield very small positional un- 
certainties. 

We note that our simulations results have a slightly dif- 



ferent outcome than those in Chen et al. (20121, who flnd 
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that the boosting factor (i.e. log(Smcasurcd/Sinput)) increases 
more rapidly as the measured S/N decreases. We do not 
measure such a rapid increase. This difference is caused by 
two effects. First, our jaclcknife map is more spatially uni- 
form, therefore sources of S/N~ 3 are approximately of the 
same flux density whereas the simulated sources in |Chen| 
et al. (20121 are at a wider range of measured flux densi- 



450/im at >5a 
at >4a 



ties, with more of them at higher flux densities than input. 
Second, we observe a deficit in submillimeter counts in this 
region of COSMOS, an effect which is described fully in the 
next section. With a deficit, source at lower S/N will not 
be boosted as much as they would in a more spatially dense 
(or populated) area. At high boosting factors, where sources 
are intrinsically more rare and less impacted by the faint 
end of the number counts, our results are consistent with 
other work. Similarly, our measured positional uncertain- 
ties, at both 450/im and 850/im, are consistent with other 
work, although the 850pim positional uncertainties are a bit 
smaller than would be expected from 850/im sources in a 
more crowded field. 

Before measuring the cumulative number counts or de- 
ciding on an appropriate S/N detection threshold, it is essen- 
tial to understand our sample's completeness and contami- 
nation rate. Completeness deteriorates at low fiux densities. 
We estimate completeness from our Monte Carlo simula- 
tions as the fraction of sources recovered >3(j or >4(j in our 
simulated maps and is shown in Figure [5] As is expected, a 
more conservative (4ct) cut on the sample results in a more 
incomplete sample at moderate fiux densities than a more 
liberal cut {3a). We will take these completeness curves into 
consideration with the calculation of the cumulative number 
counts. 

Deciding on an appropriate detection threshold requires 
an understanding of our sample's contaminants. We measure 
contamination rate as a function of detection signal-to-noise. 
At a given detection S/N, the contaminating fraction is the 
number of sources which are generated by fluctuations in 
the map due to noise or from sources well below the nominal 
detection limit. Since our maps are roughly uniform, we fix 
this limiting fiux density at the lowest deboosted fiux density 
of our 3cr sample. Any source with input fiux density below 
this threshold would not be expected to be measured in our 
data even after Eddington and confusion boosting, so here 
is considered a contaminant. The contamination rate as a 
function of both fiux density and signal-to-noise is shown in 
Figure |6] 

Figure |6] tells us the probability that a source with a 
given signal-to-noise ratio is spurious. Accordingly, we are 
able to assign a "probability of contamination" (Pcontam) 
for each source. We use this source-by-source probability to 
estimate total contamination rates for samples defined by 
different signal-to-noise detection thresholds. For example, 
there are 37 sources in a 4a cut in our 450/^m map. We esti- 
mate that the 4cr sample has a 76% likelihood of having no 
contaminants, and a 97% chance of having a contamination 
rate of <3%. This cut is quite conservative, so we push the 
detection threshold lower until we estimate contamination 
rates of 3-5%. This is our target contamination rate since it 
is low enough not to significantly impact the global proper- 
ties of the sample, while it allows for a far larger sample size 
strengthening the analysis of the population. At both 450/im 
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Figure 6. The contamination rates at 450/^m and 850/jm as a 
function of fiux density (top two panels) and signal-to-noise (bot- 
tom panel) . The contamination rate is computed as the number of 
sources which are extracted from Monte Carlo simulations which 
are (a) truly spurious, caused by noise fluctuations, or (b) at flux 
densities lower than the nominal 3cr detection limit, after cor- 
rection for deboosting. We use the measured contamination rate 
curves in the bottom panel, as a function of signal-to-noise, to 
infer the likelihood that a source of a given signal-to-noise is a 
contaminant. We plot our final signal-to-noise threshold of 3.6(t 
on the bottom panel, which shows that the maximum individual 
source contamination likelihood is ~15% in both 450/im (green) 
and 850/tm (magenta) samples. 



and 850/im, we estimate that a 3.6a detection threshold will 
yield a 3-5% contamination rate. As is seen in Figure |6j 
these limits correspond to similar 'worst' individual source 
contamination rates of ~0.14 at 450/im and ~0.17 at 850/im 
(i.e. the lowest S/N sources in the 3.6a samples have a ~14 
or 17% likelihood of being spurious). 



3.2 Number Counts 

We calculate the cumulative number counts (as has been 
done at different wavelengths in the literature: [Barger et al.| 



1999 


Smafl et al.|2002 Cowie et al.|2002||Scott et al.|2002 


Web 


b et al. 


20031 IBorys et al. 2003| |Barnard et al.| 


2004 


Coppin et a 


.||2006| IKnudsen et al.|2008 Oliver et al.] 


2010 


Chen et al.|2012|) for our sample in three stages and plot the 



cumulative number counts of the raw flux densities measured 
from our maps above >3.6a at 450/im and above >3.4a at 
850/im (flrst panel on Fig [Tj labelled 'raw'). The second 
calculation incorporates the correction for flux boosting, as 
shown in Figure |4] Here we correct the flux density of each 
individual source according to its detection signal-to-noise 
ratio. The uncertainty in the deboosting factor is propagated 
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Table 1. 450^jm-idcntificd >3.6o- point sources in COSMOS. 



Short Name RA450 DEC450 ^/N 



S45O 
Raw 
[mjy] 



S45O 
Deboosted 
[mJy] 



SOL'RCE 



Ss50 ^850 
Raw Deboosted 
[mJy] [mJy] 



SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 
SMM 



J100028 
J100033 
J100109 
JIOOOOS 
J100019 
J100023 
J100039 
J100014 
J100025 
J100016 
J100022 
J100026 
J095957, 
JlOOOOl 
J100044, 
J100004 
J100057, 
J100017 
J100007 
J09594S 
J100004 
J100017, 
J100030 
J100004, 
J100060 
J10002S 
J100016 
J095942 
J100066 
JIOOOOO 
J10004S 
JIOOOOS 
JlOOOlO 
J100021 
J095945 
J100014 
JlOOlOl 
JlOOOOl 
JlOOOlO 
JlOOOll 
J100033 
JIOOOOS 
J100036 
J10002S 
J100064 
J100056 
J100032 
J100025 
JIOOOOO 
JIOOOOO 
J095952 
J100044 
J100007, 
JlOOOlO 
J100026 
J100023 
J095959 
J100032 
J100036 
JlOOOlO 
JIOOOOS 
JIOOOIS 
J100047 
J100020 
JlOOOlO 
J100027, 
J100104 
JIOOOOO 
J095946 
J100047, 
J095962 
J095946 
J095950 
J100106 
J100041 
J100034 
J100021 
J095952 



6+023201.. 
3+022559.. 
0+022255.. 
0+022611.. 
7+023205.. 
6+022156.. 
2+022221.. 
2+023017.. 
1+02150S.. 
6+022637.. 
2+022354.. 
4+022542.. 
0+022153.. 
7+022426.. 
9+021919.. 
2+022059.. 
2+022009.. 
2+022521.. 
2+021803.. 
1+022014.. 
5+023042.. 
2+02213S.. 
8+023104.. 
6+021820.. 
1+022116.. 
5+021927.. 
6+022000.. 
9+022144.. 
7+022014.. 
5+022503.. 
3+022926.. 
4+022241.. 
2+022624.. 
3+023055.. 
2+022253.. 
4+021706.. 
0+022800.. 
5+021939.. 
3+022223.. 
2+021554.. 
1+023046.. 
6+021527.. 
9+02193S.. 
8+023336.. 
5+021919.. 
1+021831.. 
4+022148.. 
2+021930.. 
4+022223.. 
0+022524.. 
2+022133.. 
3+022313.. 
8+022306.. 
4+022026.. 
8+022318.. 
8+021751.. 
3+023000.. 
4+021802.. 
1+022152.. 
4+022024.. 
1+023059.. 
7+022813.. 
3+022049.. 
0+022129.. 
7+022707.. 
2+022448.. 
6+022633.. 
3+022902.. 
7+021837.. 
5+022520.. 
6+022711.. 
3+022931.. 
5+022016.. 
1+022454.. 
3+021640.. 
2+023421.. 
7+023114.. 
6+022258.. 



450.00 
450.01 
460.02 
460.03 
460.04 
460.05 
460.06 
460.07 
460.08 
460.09 
460.10 
460.11 
460.12 
460.13 
460.14 
460.15 
460.16 
460.17 
460.18 
460.19 
460.20 
460.21 
460.22 
460.23 
460.24 
460.25 
460.26 
460.27 
460.28 
460.29 
460.30 
460.31 
460.32 
460.33 
460.34 
460.35 
460.36 
460.37 
460.38 
460.39 
460.40 
460.41 
460.42 
460.43 
460.44 
460.45 
460.46 
460.47 
460.48 
460.49 
460.50 
460.51 
460.52 
460.53 
460.54 
460.55 
460.56 
460.57 
460.58 
460.59 
460.60 
460.61 
460.62 
460.63 
460.64 
460.65 
460.66 
460.67 
460.68 
460.69 
460.70 
460.71 
460.72 
460.73 
460.74 
460.75 
460.76 
460.77 



);2S.6 
);33.3 
-;09.0 
);08.0 
);19.7 
):23.6 
);39.2 
):14.2 
);25.1 
);16.6 
);22.2 
):25.4 
1:57.0 
);01.7 
);44.9 
);04.2 
):57.2 
):17.2 
):07.2 
);48.1 
):04.5 
);17.2 
);30.S 
04.6 
50.1 
28.5 
16.6 
42.9 
56.7 
00.5 
48.3 
08.4 
10.2 
21.3 
45.2 
14.4 
;01.0 
01.5 
10.3 
11.2 
33.1 
06.6 
36.9 
28.8 
54.5 
56.1 
32.4 
25.2 
09.4 
00.0 
52.2 
44.3 
07.8 
10.4 
26.8 
23.8 
59.3 
32.4 
36.1 
19.4 
08.1 
18.7 
47.3 
20.0 
10.7 
27.2 
04.6 
00.3 
45.7 
47.5 
52.6 
46.3 
50.5 
06.1 
;41.3 
34.2 
21.7 
62.6 



!:32:01 
!:25:59 
!;22:55 
iL26:ll 
!;32:05 
!;21:56 
!;22:21 
!:30:17 
!:15:08 
!:26:37 
!:23:64 
!;25:42 
iL21:53 
!;24:26 
!;19:19 
!;20:59 
!;20:09 
!;25:21 
!;1S:03 
!l20:14 
!;30:42 
iL21:38 
!:31:04 
i;lS:20 
!;21:16 
iLl9:27 
!;20:00 
!:21:44 
!;20:14 
!;25:03 
iL29:26 
!;22:41 
!;26:24 
!:30:55 
!:22:63 
!:17:06 
!:28:00 
!;19:39 
!:22:23 
!:15:54 
!:30:46 
!;15:27 
!;19:38 
!:33:36 
!;19:19 
!;18:31 
!;21:48 
!:19:30 
!;22:23 
!l25:24 
iL21:33 
!;23:13 
!:23:06 
!;20:26 
!:23:18 
!;17:51 
!:30:00 
!;1S:02 
!l21:62 
!;20:24 
!;30:59 
!l2S:13 
!;20:49 
!:21:29 
!:27:07 
!:24:48 
!;26:33 
!l29:02 
!:1S:37 
i;25:20 
!;27:11 
!;29;31 
i;20:16 
!;24:54 
!:16:40 
!;34:21 
i;31:14 
!;22;68 



7.80 
6.89 
6.15 
5.73 
5.62 
5.62 
5.43 
5.32 
6.16 
4.94 
4. 86 
4.77 
4.76 
4.76 
4.74 
4.60 
4.59 
4.59 
4.49 
4.47 
4.44 
4.40 
4.34 
4.33 
4.30 
4.30 
4.27 
4.26 
4.24 
4.22 
4.16 
4.11 
4.10 
4.04 
4.04 
4.01 
4.00 
4.00 
4.00 
3.99 
3.9S 
3.96 
3.94 
3.SS 
3.SS 
3.82 
3. 81 
3.80 
3. SO 
3.79 
3.7S 
3.77 
3.77 
3.74 
3.74 
3.74 
3.73 
3.73 
3.72 
3.71 
3.71 
3.70 
3.70 
3.70 
3.6S 
3.6S 
3.66 
3.65 
3.65 
3.64 
3.64 
3.64 
3.63 
3.63 
3.61 
3.61 
3.61 
3.61 



40.58±6.20 
2S.43±4.13 
47.20±7.6S 
23.47±4.10 
29.15 + 5.19 
23.14±4.12 
22.29+4.11 
23.94+4.50 
31.25 + 6.06 
20.41+4.14 

20.09 + 4.14 
19.57+4.11 

22.03 + 4.63 
19.96 + 4.20 
21.73 + 4.59 

19.92 + 4.33 
26. 18 + 6. 4S 
19.01 + 4.14 

22.63 + 5.04 
2S. 25 + 6. 33 

23.28 + 5.24 

15. 14 + 4. 12 

21.10 + 4.86 

21.95 + 6.07 

19. 72 + 4. 5S 

15. 04 + 4. 20 

17.81 + 4.17 
29.66 + 6.97 
22.76 + 5.37 
IS. 01 + 4. 27 
20.99 + 6.06 

16.96 + 4.13 
16.86 + 4.11 

15. 93 + 4. 69 

24.73 + 6.13 
20.52 + 6.12 
24.27+6.07 
19.20 + 4.80 
16.47+4.12 
23.76 + 6.96 

19.15 + 4.81 
26.42 + 6.67 

16.71 + 4.24 
24.62 + 6.34 
20.99 + 5.41 
23.57+6.16 

16.72 + 4.13 
16. 90 + 4. IS 
16.72 + 4.13 

16.29 + 4.30 

19.82 + 6.26 

16.64 + 4.16 
16.57+4.13 
16.69 + 4.19 
16.46 + 4.13 
17.06 + 4.57 
20.17+6.40 
16.75 + 4.49 

16.32 + 4.12 

16.33 + 4.13 
IS. 82 + 6. 07 
16. 48 + 4. IS 

16.48 + 4.46 
16.23 + 4.12 
16.18 + 4.12 
16.13 + 4.12 
23.50 + 6.42 
17.79 + 4.87 
2S. 64 + 7. 84 
16.52 + 4.26 
IS. 86 + 5. IS 

25. 16 + 7. 74 
21.50 + 6.92 

23. 49 + 6. 4S 
20.01 + 6.54 

25. 05 + 7. 76 

17.34 + 4.80 
17.85 + 4.96 



37.54 + 6.58 
26.62 + 4.96 

41.70 + 9.00 
20.46 + 4.78 
26.35 + 6.04 
20.12 + 4.80 
19.32 + 4.76 

20.65 + 6.19 

26.71 + 6.93 

17.26 + 4.68 
16.91 + 4.67 
16.40 + 4.62 

15. 45 + 6. 21 

16.72 + 4.73 
IS. 18 + 6. 17 
16.53 + 4.88 

20.89 + 6.19 
16.76 + 4.67 

15. 66 + 6. 68 
23.24 + 7.14 
19.11 + 6.91 
14.83 + 4.64 
17.15 + 6.47 
17.81 + 6.71 

16.97 + 6.16 
14.60 + 4.73 
14.38 + 4.69 
23.91 + 7.84 

15. 32 + 6. 04 

14.46 + 4.80 
16.76 + 6.69 
13.46 + 4.69 
13.37 + 4.67 

14.90 + 6.36 
19.46 + 7.01 
16.10 + 6.87 
19.03 + 6.96 
16.05 + 6.51 

12.91 + 4.73 
IS. 60 + 6. 84 

14.98 + 6.53 
20.62 + 7.67 
13.01 + 4.89 
19.00 + 7.36 
16.18 + 6.28 

17.92 + 7.22 
11.90 + 4.86 
12.03 + 4.91 

11.89 + 4.86 

12.27 + 6.07 

14.90 + 6.19 
11.74 + 4.90 

11.68 + 4.88 

11.69 + 4.97 
11.51 + 4.90 

12.70 + 6.42 
14.98 + 6.42 
12.44 + 6.34 
11.35 + 4.90 

11.33 + 4.92 
13.90 + 6.04 
11.42 + 4.98 
12.15 + 5.30 
11.23 + 4.90 

11.14 + 4.91 
11.08 + 4.92 

17.15 + 7.68 
12.96 + 6.83 
20.86 + 9.38 
11.27 + 6.11 
13.68 + 6.21 
20.42 + 9.28 
16.58 + 7.09 
17.00 + 7.76 
14.44 + 6.64 
20.23 + 9.31 
12.50 + 6.76 
12.85 + 5.94 



1.06 
1.51 
1.48 
1.45 
1.48 
1.48 
1.56 
1.69 
1.67 
1.76 
1.78 
1.80 
1.81 
1.81 
1.81 
1.83 
1.84 
1.84 
1.86 
1.86 
1.87 
l.SS 
1.90 
1.90 
1.91 
1.91 
1.92 
1.92 
1.93 
1.93 
1.95 
1.97 
1.97 
1.99 
1.99 
2.00 
2.00 
2.00 
2.01 
2.01 
2.01 
2.02 
2.02 
2.06 
2.06 
2.13 
2.16 
2.16 
2.16 
2.18 
2.19 
2.20 
2.20 
2.23 
2.23 
2.24 
2.25 
2.26 
2.26 
2.27 
2.27 
2.28 
2.28 
2.29 
2.31 
2.32 
2.34 
2.34 
2.36 
2.36 
2.36 
2.37 
2.37 
2.38 
2.39 
2.39 
2.39 
2.40 



<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
<0.01 
0.02 
0.02 
0.03 
0.03 
0.03 
0.03 
0.03 
0.03 
0.04 
0.04 
0.04 
0.04 
0.05 
0.05 
0.06 
0.07 
0.07 
0.07 
0.07 
0.07 
0.08 
0.08 
0.08 
0.08 
0.08 
0.09 
0.09 
0.09 
0.09 
0.09 
0.09 
0.09 
0.09 
0.10 
0.10 
0.10 
0.10 
0.10 
0.11 
0.11 
0.11 
0.11 
0.11 
0.11 
0.11 
0.11 
0.11 



860.07 
860.02 

860.00 
850.03 
860. OS 
850.12 
860. 9S 
860.46 
860.34 
850.135 
850.41 

850.32 

850. 2S 

850.42 

850.24 
850.149 



50.116 
^60. 06 



860. 6S 
860.14 
850.90 
860. 4S 



860.47 
860.100 

860.109 
860.96 
850.06 



860.21 
860.101 



9.94+1.06 
9. 47+0. SO 
9.67+1.71 
16. 15 + 0. SO 
12.11 + 1.03 
7. 19 + 0. SO 
6. 49 + 0. SO 
3.1S + 0.SS 
6.62 + 1.16 
4. 41 + 0. SO 
2. 60 + 0. SO 
4. 12 + 0. SO 
1.71 + 0.S9 
4.66 + 0.S2 

1.28 + 0.91 
4.S6 + 0.S4 
8.99 + 1.14 
4. 12 + 0. SO 
-0.76 + 0.96 
0.77+1.16 
6.26 + 1.01 
0.72 + 0. SO 
3.04 + 0.96 
-0.36 + 0.96 

1.29 + 0.91 
-0.82 + 0.82 
-0.60 + 0.81 

4.26 + 1.26 
12.24 + 1.13 

0.94 + 0. S3 

2.23 + 1.04 
1.17+0. SO 
4. 87+0. SO 
0.59 + 0.93 

3.30 + 1.13 
-0.55 + 0.97 
5.67+1.31 
-0.98 + 0.91 
6. 25 + 0. SO 
-0.70+1.11 
3.52 + 0.96 
1.39 + 1.26 
3. 99 + 0. S3 
1.62 + 1.2S 

-0.15+1.12 
-1.21 + 1.30 
0.11 + 0. SO 
1.60 + 0. SI 

2. 73 + 0. SO 
-0.72 + 0.84 

4.74 + 0.9S 
2. 91 + 0. SI 

-0.46 + 0.80 
2.84 + 0.S2 
2. 93 + 0. SO 
9.30 + 0.S9 

-0.89+1.03 
0.56 + O.SS 
5. 29 + 0. SO 
2. 87+0. SO 
1.81 + 0.99 

1.24 + 0. SI 
-0.18 + 0.89 
0.05 + 0. SO 
0.86 + 0. SO 
0.00 + 0. SO 
4.29 + 1.36 
1.08 + 0.94 

-0.56+1.46 
0.72 + 0.S4 
1.28 + 0.99 
1.44 + 1.46 
0.12 + 1.09 
-1.96+1.39 
-0.38+1.10 
-1.80+1.56 
0.64 + 0.96 
0.11 + 0.96 



9.21 + 1.46 

5. 22 + 0. 70 
<6.13 

10.86 + 0.54 
10.54 + 0.90 
6.60+1.12 
6.83+1.13 
2.43+1.07 
4.64+1.40 
3.68+1.01 
1.87+1.01 
3.43 + 0.9S 

<2.68 
3.81 + 1.04 

<2.74 
4. 07+1. OS 

<3.43 
3.43 + 0.99 
<2.86 
<3.47 
6.32+1.34 

<2.40 
2.16+1.21 
<2.88 
<2.74 
<2.45 
<2.44 
3.14+1.56 
10.92+1.02 
<2.49 
<3.11 
<2.40 
4.12+1.06 
<2.79 
<3.39 
<2.91 
4.60+1.57 

<2.74 
6.56+1.11 

<3.32 
2.71 + 1.16 

<3.74 
3.30+1.00 
<3.83 
<3.37 
<3.89 
<2.41 
<2.44 
<2.40 
<2.52 
3.92+1.19 
2.22 + 0.9S 

<2.40 
2.13+1.00 
2.24 + 0.97 
S. 42 + 0. 92 
<3.10 
<2.65 
4. 60+1. OS 
2.19 + 0.9S 
<2.96 
<2.43 
<2.66 
<2.40 
<2.41 
<2.40 
3.06+1.70 
<2.82 
<4.38 
<2.53 
<2.97 
<4.34 
<3.26 
<4.17 
<3.31 
<4.65 
<2.86 
<2.84 



Table Notes. The >3.6cr 450/im— detected sources we extract within the central 394'^^ of our COSMOS map. The 3.6cr detection threshold is chosen based 
on an estimated 3—5% contamination rate. The list is ordered by detection signal-to-noise ratio (S/N). The "Raw" flux densities are those measured directly 

from our map. Th e "D eboosted" flux densities are those given after correction for confusion and Ejddington boosting as a function of detection S/N, as 
described in section |3.1| Wc also measure a 90% confidence interval for positional uncertainties and estimated probability of contamination, fcontam, from 
the results of our Monte Carlo tests as functions of detection S/N. In the last three columns, wc list the corresponding 850/Am counterparts given in Table [5] 

and fluxes: non-detection deboosted fluxes arc quoted at 3(T limits. 



to the new uncertainty in the sources' deboosted flux den- 
sity. As a result, the sources' signal-to-noise ratios appear to 
be lower, which is true for the uncertainty on the flux density 
measurements but is not true for the detections' signal-to- 
noise. The net effect of deboosting is a shift towards lower 
flux densities (in the x-direction on the number counts plot), 
seen in the middle panels of Figure [7| 



The final effect which needs to be accounted for in our 
computation of the cumulative number counts is incomplete- 
ness and contamination in our samples. As shown in Fig- 
ure [5] and Figure |6j our samples are not complete down to 
our nominal ~15 mjy«3. 60-450 or ~2.9 mJy«3.6cr85o detec- 
tion threshold. Furthermore, 3-5% of our sources are likely 
contaminants. To correct for sources missing and contami- 
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nant sources we construct a correction factor, Ci, for each 
source, based on its flux density and signal-to-noise, such 
that Ci = (1 - fcontam)/ f complete- We correct the number 
counts accordingly. Not surprisingly, this correction for in- 
completeness and contamination has an effect on the low 
flux density bins and a minimal effect on the high flux den- 
sity bins. Our deboosted and corrected number counts are 
given in Table [2] 

In functional form, we can assume that the differential 
number counts follow a few different forms. The types we 
adopt in this paper are a broken double power law (as most 
recently used in Chen et al.^2 012J, taking the form: 



dN 

lis 



S<So 
S>So 



_. , ... o (1) 

5o ^ So ^ 

where So represents the turnover flux density, No a normal- 
isation for the function in source density per mjy, a is the 
faint end slope and /3 is the bright end slope (presumed to 
be steeper). Alternatively we can assu me the number coun ts 
follow a Schechter function, as used in Geach et al. (2012|^ 



dN 

Is 



No 
S 



\SoJ 



(2) 



A Schechter function fit is inadequate at representing the 
high-flux density end, particularly when we consider the 



500/im results from Oliver et al. ( 2010 1 and Clements et al 



(2010a I. With this in mind, a double power law is more ap- 



propriate across a wide range of flux densities, especially at 
>40 mJy. However, a Schechter function might be appropri- 
ate in the flux density range addressed in this paper. The 
sources at higher flux densities could be dominated by a dif- 
ferent underlying population or process (e.g. gravitational 
lensing). In this paper we offer fits to both Schechter and 
a double powerlaw. With a more intimate understanding of 
the underlying population at high flux densities one might 
be able to disentangle which fit provides a more appropri- 
ate physical interpretation. Until then, it is not realistic to 
claim we can constrain the underlying shape of the num- 
ber counts; here we offer analysis and solutions with both 
parametrisations . 

We perform a optimisation to determine best-fit pa- 
rameter values at 450pim and 850/im respectively in equa- 
tions [T] and [2] which are given in Table [2] These best-fit 
parameter values are tested for consistency with the Monte 
Carlo simulations; in other words, we adjust the fit parame- 
ters, populate the jackknife maps accordingly, then extract 
sources and re-measure the raw number counts. The parame- 
ters which reproduce the raw number counts most accurately 
are used in the Monte Carlo tests for boosting factor, posi- 
tional uncertainty and estimating contamination and com- 
pleteness rates, as given in the previous section. Once we 
arrive at the corrected number counts, we then remeasure 
our best-fit parameters with a test which agrees with our 

^ Note that Equation 1 of Geach et al. should be written 
dN/dS = iN'/S')(S/S'y-"exp(-S/S'), and the best-fit param- 
eter N' from their paper should be quoted as A*^' = 4900 it 
1040deg-2 mjy-l rather than N' = 490 ± 104 deg-^. Also note 
we use an exponent of —a instead of 1 — o. By our equation, the 
parameters of the Geach et al. data are: No = 4900 deg~^ mjy~^. 
So =10mJy, and a =2.0. 



Table 2. True Number Counts at 450/xm, 850/xm, and best-fit 
double power-law/Schechter function parameters. 



450^m-Number Counts 

S450 N dN/dS 

(mJy) (mJy-1 dog" 



850^m-Number Counts 



N(> S) 
(dcg-2) 



SgSO N 
(mJy) 



dN/dS 



N{> S) 



(dog 



10.59 15 
14.13 27 
18.84 29 



129 
36.1 



1430+370 
679+;!= 



+ 15.0 
12.2 



62 



25.12 
33.50 



-2.58 
_+1.80 
-0-63 



215 

48+f. 
18 



+ 32 
-12 



2.65 
3.37 
4.29 
5.46 
6.96 
8.86 
11.28 
14.37 



579 
412 
306' 
78.2 
40.4 



+ 157 
124 
+ 104 



83 



10403 
719+ 



+ 70 
-57 
+ 31.0 
22.2 
+ 19.7 
13.2 

13 3+" + 
-^•^"^ 6.0 

13.9^" " 



□ + 9.6 
-5.7 



2.7' 



+ 5.1 



r,+300 

-230 
+ 189 
-150 
380+lf 

issieo 

f. + 101.5 
""■'^-33.2 

38-0-20.3 
9 5+29-6 



Best-Fit Paramaters to Equation [T] 

450/im Parameters 850^m Parameters 

No = (1.4±0.4)xl0^ (1.5±0.4)xl03 

So = 18±8 4.3±1.8 

a = 1.91+0.40 1.34+0.09 

/3 = 5.5+1.0 4.50+0.38 

Best-Fit Paramaters to Equation [2] 

450^m Parameters 850^im Parameters 



JVo = 
So = 



(1. 1+0.4) xlO*' 
6.1 + 1.0 
0.5+1.3 



(3.3+1.4) X 10-^ 
3.7+1.1 
1.4+1.1 



Table Notes. Number counts are given in difli'crcntial and 
cumulative form. The number of sources contributing to each 
bin is given as the number "N." The best-fit parameters to the 
number counts are given at bottom and are fit to the differential 
counts. 

initial injected counts within uncertainty. We find the power- 
law fit produces a better flt to our data than a Schechter fit, 
so we use the double power-law Monte Carlo predictions 
throughout. 

We measure a counts deficit in this ~400'^ region of 
COSMOS relative to other nearby regions (e.g. AzTEC mea- 
surements just to the north; Scott et al. 20081. Although 



unexpected, this deficit is easily explained; it is also seen in 



the smaller yet deeper coverage of Geach et al. (20121. The 



Chandra Deep Field South (CDFS) has a similar deficit of 
870 Aim sources, as found by LABOCA ( |Weifi et al!1[2009l ). 
Furthermore, the deficit only seems to be prominent at the 
highest fiux densities where sources are intrinsically more 
rare and prone to cosmic variance on small scales. 



3.3 Source Lists 

The 450/im sources detected >3.6(t are summarised in Ta- 
ble [l] and the 850fim sources detected at >3.6ct are sum- 
marised in Table [s] (given at the end of the paper due to 
its length) . Both raw and deboosted flux densities are given 
in the tables, along with a 90% confldence interval on po- 
sitional uncertainty, A{a,S), the probability of contamina- 
tion or spuriousness for that source (Pcontam), and the cor- 
responding source detected at the opposite wavelength if it 
exists. 

Of the 78 450^m sources, 24 have 850/im counterparts 
detected above 3.6(j (a further 7 sources are detected at 
850fj,m above 3a). Such a low fraction of 450/im corre- 
sponding to 850/im sources ( ^ 31/78~39%) suggests that 
indeed, the 450/im population is intrinsically different than 
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Figure 7. The cumulative number counts at 450/im (left column) and 850/xm (right column) from our COSMOS data. The top row 
shows the raw number counts, with extracted flux densities from our map, the middle row corrects the sources' flux densities for confusion 
and Eddington boosting, and the bottom row corrects the deboosted number counts for incompleteness and contamination (as discussed 
in the text). The best-fit double power law parametorisation of the number counts are shown in gray in both bottom panels, with 
associated light gray uncertainties. Wc compare our measured counts (black open and closed circles) to those in the literature. At 450/xm, 
we compare to the recent results of |Chen et al.H2012[ | (blue asterisks, blue dashed line) and |Geach et ar] | |2012[ | (red squares), and we 
note that the latter work is data taken in the same sky region as our COSMOS map although covering only ~l/4 the area. We also 
compare 450/im data to 500/xm data from the HerschelSY'lKE, instrument from |01iver et aT| f2010) (green triangles) and |CIements et al.| 
| |2010b| l (orange squares). The |Bethermin et al. (2012 work is not from direct 500/im extractions but rather stacks on 24/im sources at 
500/im. At 850/im, we overplot the [Chen et aL 20121 results and the best-fit models from [Knudsen et al.| l [2b08[ l; we omit many data 
measurements taken at 850/im to avoid clutter on the plot. 
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the 850^tm population, either probing a different redshift 
regime, luminosity regime, or dust temperature regime than 
the 850/im population. 



[Chen et al. (20121 show that sources detected at both 
450/im and 850/xm at low-S/N (say 3-4(t) have effective 
S/N higher than the individual S/N in either band, in other 
words, that a source which is 3a" at 450/xm and Scr at 850/im 
is effectively a Aa detection. We use this technique to ex- 
tract 450/im sources below our nominal detection threshold. 
The additional 3< a <3.6 450/im sources we detect with 
850/im counterparts are given in Table [3] We exclude these 
sources from our 450/im number count computation due to 
the unquantified bias which requiring a 850/im counterpart 
enforces on the sample, but we can use this sample for com- 
parisons to the independently-detected >3.6ct 450/im sample 
in redshift, luminosity and dust temperature space. 



3.4 Counterpart Matching 

Determining the optical/near-infrared counterparts to sub- 
millimeter sources is complex, but thanks to the low posi- 
tional uncertainty on our 450/im sources, it is a much more 
straightforward process than in other submillimeter maps at 
other wavelengths or taken at other observatories. A hand- 
ful of sources in our sample (five sources at 450/im, eight 
at 850/im) have interferometric observations either with the 
SMA ( [Younger et al.|2007||2009[) or the Plateau de Bure In- 
terferometer ( Smolcic et al.|2012 1 where the counterparts to 
submm observations are known, but the rest of the sample 
needs careful counterpart analysis. Our counterpart match- 
ing procedure at 450/im and 850/im is the same at both 
wavelengths, but we note that the added uncertainty at 
850/im generates many more potential counterparts and is, 
therefore, naturally more uncertain. Here we describe our 
counterpart matching procedure and discuss the selection of 
the best-guess multiwavelength counterpart (s). 

All optical (i-band), near-infrared (3.6/im), 24/im, and 
1.4 GHz radio sources which fall within the 90% confidence 
positional uncertainty are considered as possible counter- 
parts. Optical sources are taken from the Ilbert et al. cat- 
alogue, which identifies individual sources down to i ~27. 
Although we consider sources within the nominal positional 
uncertainty found by our simulations (as shown in Figure|4|, 
additional accommodation must be made for uncertainty in 
the positions measured at other wavelengths or source size. 
Since ScuBA-2 450/im positional uncertainties are on the or- 
der of 1-2", a source's size or positional uncertainty at other 
wavelengths is not negligible. At optical and near-infrared 
wavelengths we are able to estimate sources' 'positional un- 
certainty' by considering its measured size; we fix the po- 
sitional uncertainty for each source individually by taking 
the area (measured in arcseconds squared) and elongation 
(both measured using SExtractofQ to convert to a semi- 
major axis. Most optical/near- infrared sources have approx- 
imated semi-major axes ~0.5", with 90% of sources being 
<1". Therefore, when seeking potential matches to 450/im 
sources, we search within a radius of the summed 450/im 



source's positional uncertainty and the optical/near-infrared 
sources' approximated semi-major axis. 

We estimate the positional uncertainty of 24/im and 
1.4 GHz sources using a more statistical method (since the 
resolution is not adequate to resolve source sizes) by measur- 
ing the global offset of 24/im and 1.4 GHz sources to A'-band 
and j-band positions. We find that 90% of 24/im sources are 
within 1.05" of their associated optical/near-infrared coun- 
terpart and 1.4 GHz are within 0.85". Therefore, we add 
1.05" and 0.85" to our search radius for 24/im and 1.4 GHz 
counterparts, respectively. To test the validity of this 'ex- 
tension' of the positional uncertainty, we observe many in- 
stances of 24/im and 1.4 GHz counterparts falling in the ~l" 
beyond the 450/im positional uncertainty but not within; 
the likelihood of this happening is slim, especially at radio 
wavelengths where the sources density is so rare and chance 
alignment for multiple sources is unlikely. 

In addition, radio and mid-infrared sources are more 
likely to be infrared-luminous than random field galaxies. 
The advantage of identifying radio counterparts comes from 
the well-known FIR/radio correlation for starburst galaxies 
( Helou et al.|1985,^Condon|1992 1, whereby a radio-detected 
galaxy is likely to be FIR-bright, and thus correspond to the 
source generating the FIR emission detected in our S CUBA- 2 
maps. MIPS 24/tm-emission also correlates with FIR emis- 
sion, although the relationship is more complex due to the 
variation of SED types in the mid-infrared (PAH emission 
features drifting in and out of the band, and AGN gener- 
ated power law emission) and the increased source density 
at 24/im as opposed to 1.4 GHz. The existence of a mid- 
infrared or radio counterpart within the beamsize of 450/tm 
observations can help distinguish that source from others as 
the likely counterpart. 

A similar correlation has been studied between near- 



infrared luminosity and far-infrared luminosity (e.g. Wang 
|et al.|2006| , which implies a correlation between stellar mass 
and star formation rate (this is of course observed in most 
'normal' star-forming galaxies and has been dubbed the 
'main sequence,' as in [Noeske et al.|[2007l [Daddi et al.| 
|2007[ ; however it is not clear that this has any bear- 
ing on starbursting systems). Pursuing an exclusive near- 
infrared-driven counterpart selection has been used for other 
submillimeter populations (a very thorough analysis of its 
success is given in [Alberts et alT 20131. Here we pursue 



both optical and near-infrared matches since many of our 
sources lack 3.6/im counterparts, implying that a poten- 
tial bias is introduced when only near-infrared counterparts 
are considered. As we find later in our analysis and discus- 
sion, 450/im sources are on average less infrared luminous 
than their 850/im or 1.1-1.2 mm-luminous cousins so near- 
infrared counterpart matching might have a different bias 
or be less effective (it also has to do with the depth of both 
catalogues); we leave a detailed analysis of 450/im counter- 



part matching strategy (as was done for 1.1mm in Alberts 
et al.|20T3 1 to a future work 



^ SExTRACTOR: a Source Extractor program developed in |Bertinl 
& Arnoutsl (119961. 



A summary of counterparts (24/tm, 1.4 GHz and 
optical/near-IR) which fall in the appropriate search radius 
is found in Table [t] (for >3.6ct 450/im sources). Table [s] (for 
3< (7 <3.6 450/im sources), and Table|9](for >3.6(t 850/tm 
sources which are not 450/im- luminous), all given at the end 
of this paper. Photometric redshifts ( [Ilbert et al.||2010 l are 
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Table 3. Marginal 3< a <3.6 450^m-identified sources with 850//m >3cr Counterparts. 



Name 850/im- 450Mm- RA450 DEC450 S/Ngso S/N450 S450 S450 Ssso Ssbo 

Shore Short RAW Deboosted Raw Deboosted 

Name Name [mjy] [mjy] [mjy] [m Jy] [ ] 



SMM J 100031. 0+022751.. 


850 


62 


450 


78 


10 


00 


31 





02 


27 


51 


4.29 


3 


60 


14 


94+4 


15 


10 


.72+4 


99 


3 


46+0.81 


2.80+0.96 


2.42 


SMM J 100034. 3+022 121.. 


850 


10 


450 


81 


10 


00 


34 


3 


02 


21 


21 


8.60 


3 


58 


14 


71+4 


11 


10 


.47+4 


98 


6 


90+0.80 


6.24+1.13 


2.45 


SMM J095950. 8+022745.. 


850 


55 


450 


86 


09 


59 


50 




02 


27 


45 


4.51 


3 


55 


19 


80+5 


58 


13 


.95+6 


80 


4 


77+1.06 


3.89+1.27 


2.50 


SMM J 100025. 1+021847.. 


850 


09 


450 


87 


10 


00 


25 


1 


02 


18 


47 


8.89 


3 


55 


15 


24+4 


30 


10 


.74+5 


24 


7 


47+0.84 


6.83+1.17 


2.50 


SMM J100022. 2+022842.. 


850 


23 


450 


94 


10 


00 


22 


2 


02 


28 


42 


6.24 


3 


50 


14 


67+4 


19 


10 


.14+5 


17 


5 


06+0.81 


4.31+1.07 


2.58 


SMM J 100052. 8+02 1906.. 


850 


103 


450 


95 


10 


00 


52 


8 


02 


19 


06 


3.54 


3 


50 


18 


59+5 


31 


12 


.83+6 


56 


3 


88+1.10 


2.94+1.34 


2.59 


SMM J100049. 8+022448.. 


850 


133 


450 


96 


10 


00 


49 




02 


24 


48 


3.25 


3 


49 


15 


25+4 


36 


10 


.50+5 


39 


2 


85+0.88 


2.06+1.09 


2.60 


SMM J095950. 7+022823.. 


850 


33 


450 


99 


09 


59 


50 


7 


02 


28 


23 


5.54 


3 


47 


20 


15+5 


SO 


13 


.76+7 


20 


6 


20+1.12 


5.18+1.41 


2.63 


SMM J 100018. 6+022356.. 


850 


128 


450 


103 


10 


00 


18 


6 


02 


23 


56 


3.29 


3 


44 


14 


20+4 


13 


9 


58+5. 


15 


2 


63+0.80 


1.92+1.00 


2.69 


SMM J100025. 0+022757.. 


850 


85 


450 


105 


10 


00 


25 





02 


27 


57 


3.73 


3 


44 


14 


24+4 


14 


9 


60+5. 


16 


2 


99+0.80 


2.32+0.97 


2.69 


SMM J 100041. 8+022358.. 


850 


92 


450 


106 


10 


00 


41 


8 


02 


23 


58 


3.67 


3 


44 


14 


19+4 


13 


9 


57+5. 


15 


2 


94+0.80 


2.27+0.97 


2.69 


SMM J 100036. 9+022 124.. 


850 


108 


450 


110 


10 


00 


36 


9 


02 


21 


24 


3.48 


3 


42 


14 


07+4 


12 


9 


41+5. 


15 


2 


79+0.80 


2.10+0.98 


2.73 


SMM J100032. 0+023049.. 


850 


90 


450 


113 


10 


00 


32 





02 


30 


49 


3.67 


3 


40 


16 


37+4 


81 


10 


.88+6 


03 


3 


52+0.96 


2.71 + 1.16 


2.75 


SMM J 100025. 6+023051.. 


850 


159 


450 


126 


10 


00 


25 


6 


02 


30 


51 


3.05 


3 


36 


15 


79+4 


70 


10 


.32+5 


91 


2 


89+0.95 


2.02+1.19 


2.83 


SMM J 100029. 5+022131.. 


850 


131 


450 


133 


10 


00 


29 


5 


02 


21 


31 


3.27 


3 


32 


13 


67+4 


12 


8 


78+5. 


19 


2 


64+0.81 


1.91 + 1.00 


2.90 


SMM JlOOOll. 8+022935.. 


850 


52 


450 


134 


10 


00 


11 


8 


02 


29 


35 


4.73 


3 


31 


14 


45+4 


36 


9 


26+5. 


50 


4 


09+0.86 


3.37+1.04 


2.91 


SMM J095953. 0+022641.. 


850 


163 


450 


135 


09 


59 


53 





02 


26 


41 


3.00 


3 


31 


16 


78+5 


07 


10 


.76+6 


39 


2 


91+0.97 


2.02+1.22 


2.91 


SMM JIOOOIO. 1+021758.. 


850 


S3 


450 


166 


10 


00 


10 


1 


02 


17 


58 


3.77 


3 


24 


15 


87+4 


90 


9 


79+6. 


20 


3 


52+0.93 


2.74+1.12 


3.08 


SMM J095957.5+022726.. 


850 


01 


450 


168 


09 


59 


57 


5 


02 


27 


26 


13.00 


3 


23 


15 


24+4 


72 


9 


36+5. 


97 


11 


.98+0.92 


10.02+0.77 


3.10 


SMM J100005. 4+022516.. 


850 


104 


450 


173 


10 


00 


05 


4 


02 


25 


16 


3.51 


3 


22 


13 


31+4 


13 


8 


11+5. 


23 


2 


82+0.80 


2.12+0.98 


3.13 


SMM J100022. 6+023023.. 


850 


88 


450 


179 


10 


00 


22 


6 


02 


30 


23 


3.69 


3 


21 


14 


51+4 


52 


8 


80+5. 


73 


3 


31+0.90 


2.56+1.08 


3.15 


SMM JIOOIOO. 5+021946.. 


850 


105 


450 


185 


10 


01 


00 


5 


02 


19 


46 


3.50 


3 


19 


20 


06+6 


30 


12 


.00+7 


97 


4 


48+1.28 


3.37+1.57 


3.21 


SMM J100014. 1+022836.. 


850 


114 


450 


189 


10 


00 


14 


1 


02 


28 


36 


3.40 


3 


17 


13 


37+4 


22 


7 


92+5. 


33 


2 


81+0.83 


2.08+1.02 


3.25 


SMM J100005. 0+021718.. 


850 


49 


450 


193 


10 


00 


05 





02 


17 


18 


4.79 


3 


16 


17 


40+5 


50 


10 


.27+6 


95 


4 


91+1.02 


4.05+1.24 


3.27 


SMM JlOOOll. 1+021507.. 


850 


53 


450 


206 


10 


00 


11 


1 


02 


15 


07 


4.68 


3 


14 


20 


69+6 


60 


12 


.04+8 


32 


5 


72+1.22 


4.70+1.47 


3.33 


SMM J100023. 6+021916.. 


850 


27 


450 


215 


10 


00 


23 


6 


02 


19 


16 


6.08 


3 


11 


13 


05+4 


20 


7 


48+5. 


28 


5 


00+0.82 


4.23+1.08 


3.39 


SMM J100032. 0+023324.. 


850 


67 


450 


240 


10 


00 


32 





02 


33 


24 


4.07 


3 


06 


19 


01+6 


21 


10 


.61+7 


75 


5 


15+1.26 


4.12+1.50 


3.51 


SMM J100041. 3+022534.. 


850 


151 


450 


247 


10 


00 


41 


3 


02 


25 


34 


3.12 


3 


05 


12 


58+4 


13 


6 


97+5. 


14 


2 


52+0.81 


1.79+1.02 


3.54 


SMM J100013. 4+022224.. 


850 


40 


450 


252 


10 


00 


13 


4 


02 


22 


24 


5.20 


3 


04 


12 


57+4 


14 


6 


91+5. 


15 


4 


16+0.80 


3.47+0.99 


3.57 


SMM J100049. 4+023028.. 


850 


76 


450 


253 


10 


00 


49 


4 


02 


30 


28 


3.83 


3 


04 


16 


73+5 


51 


9 


20+6. 


85 


4 


36+1.14 


3.42+1.36 


3.57 



Table Notes. Sources extracted down to 3cr in the 450^m map using 850^m >3fT positional priors. The 450^tm >3(T list (with 274 sources) was checked 
against the 850^m >3(T list (with 164 sources) and 30 sources were found. The positions here arc given by their 450^m detection with associated positional 
uncertainty. The short names are numbered in order of decreasing S/N for all sources extracted in the map, e.g. from 450.00 at S/N= 7.80 to to 450.273 

with S/N— 3.00, hence the high numbers given at both 450^m and 850/^m. 



quoted where computed. In the case that a candidate coun- 
terpart is detected in the near-infrared and not the opti- 
cal, its i-band magnitude is quoted as "— " and position 
given as the IRAC channell (3.6 /im) position. Also, the 
handful of sources which have interferometric observations 
have their known counterparts in bold in Tables [7j [S] and |9] 
The interferometrically-observed subset of ScuBA-2 sources 
(along with all ScuBA-2 sources which were detected with 
other submillimeter instruments) is given in Table |4] 

In the case of 450/.tm sources, the best-guess 
optical/near-infrared counterparts are those which are clos- 
est to the Scuba- 2 position and those are the counterparts 
we assume for the rest of this paper. While assuming one sin- 
gle counterpart puts us at a disadvantage for understanding 
the nature of submillimeter-luminous galaxy pairs (or pro- 
jected pairs) as found by [Barger et al.| ( |2012[ ) and |Karim| 
et al. (2012 I, these observations have no way of constraining 



tor of 1.5x larger (this being a prominent ongoing issue 
with submillimeter-selected galaxy samples; e.g.|Smail et al. 



multiples. In addition, our 450/.im sources are far less likely 
to be multiples than the larger beamsize 850/im sources. 
For lack of a better assumption, we assume single, closest 
counterparts for the whole sample. 



Figure 18 (given at the end of this paper) show 12" x 12" 
postage stamp cutouts of the 450jUm sources in the optical 
{Biz tricolor) and near-infrared (IRAC 3.6//m) with radio 
and 24/im counterparts identified for reference. 

Our 850/im counterpart matching is more uncertain 
given the larger positional uncertainties at 850/im than at 
450/tm. The number of optical and near-infrared counter- 
parts falling into the search radius is on average a fac- 



2006 



19971 jBlain et alT] [20021 [Pope et al.| [20061 | Wang et al 



2011[ [Barger et al.|[2012[ [Karim et al.[[2012[ ). Since 



850/im counterparts are more numerous and uncertain than 
450/im counterparts, we take advantage of the correlation 
between FIR emission and radio/mid-infrared wavelengths 
and choose the best-guess counterpart as the nearest 24/tm 
or 1.4 GHz source, if it exists. Otherwise we select the near- 
est optical/near-infrared source. We tested the difference be- 
tween these samples (redshift distribution, mean luminosity, 
dust properties) and found it does not greatly impact the 
bulk result; of course more detailed analysis on individual 
sources requires interferometric observations (as in jWang 
et al.|2011[|Barger et al.|2012[[Karim et al.|2012[ ) to unequiv- 
ocally identify the multiwavelength source(s) generating the 
submm emission. 



A significant fraction of sources, both at 450/im and 
850/im lack 24/im or radio counterparts. At 450/im, only 
40% of sources have 24/im or radio counterparts and at 
850/im, only 60%. This statistic should improve with deeper 
radio data, which should be available in the coming years in 
the COSMOS field, but it would not fix the bias introduced 
in submillimeter samples which are reliant on these counter- 
parts. In fact, studies which are reliant on identification of 
mid-infrared and radio counterparts at similar luminosities 
are likely significantly biased. 
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Table 4. Identification of ScuBA-2 Sources at 870/jm-1.2 mm 
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450/jm- 


SSO^m- 


Alternate 




Sa 


INTERFER- 


Inst./ 


Reference 


Name 


Name 


Name 


[mm] 


[mjy] 


OMETRIC? 


Tel 





450.00 
450.01 



850.07 
850.02 



450 


07 


850.98 


450 


08 


850.46 


450 


16 




450 


20 


850.24 


450 


27 


850.115 


450 


28 


850.05 


450 


38 


850.14 


450 


55 


850.06 


450 


66 


850.146 


450.168 


850.01 



AzTEC .1100028.94+023200. 3 1.1 3.8±0.9 

AzTEC/C80 1.1 4.1±0.9 

COSLA-47 0.87 9.0±2.8 

COSLA-47 1.3 3.11±0.59 

AzTEC .J100008. 03+022612.1 1.1 8.3±1.1 

AzTEC-2 0.89 12.4±1.0 

AzTEC/C3 1.1 10.5±1.0 

AzTEC .1100019.73+023206.0 1.1 6.5±1.3 

AzTEC-5 0.89 9.3±1.3 

AzTEC/C42 1.1 4.8±1.1 

AzTEC/C38 1.1 5.1±1.0 

COSLA-35 0.87 8.2±2.2 

COSLA-35 1.3 2.15±0.51 

AzTEC/C169 1.1 3.1±0.8 

COSLA-8 0.87 6.9±1.6 

COSLA-8 1.3 2.65±0.62 

C0SB03** 1.2 7.4±1.1 

AzTEC/C6** 1.1 9.6±1.0 

AzTEC .1100004.54+023040.1 1.1 3.3±0.8 

AzTEC/C150 1.1 3.3±1.2 

AzTEC/C65 1.1 4.4±1.0 

C0SB03** 1.2 7.4±1.1 

AzTEC/C6** 1.1 9.6±1.0 

AzTEC/C24 1.1 5.7±1.0 

C0SB07 1.2 5.0±0.9 

AzTEC/C160 1.1 3.1±1.2 

AzTEC/C66 1.1 4.3±0.9 



Y 
Y 



AzTEC 

AzTEC 

LABOCA 

PdBI 

AzTEC 

SMA 

AzTEC 

AzTEC 

SMA 

AzTEC 

AzTEC 

LABOCA 

PdBI 

AzTEC 

LABOCA 

PdBI 

MAMBO 

AzTEC 

AzTEC 

AzTEC 

AzTEC 

MAMBO 

AzTEC 

AzTEC 

MAMBO 

AzTEC 

AzTEC 



Scott et al 



]]2008[ 

Aretxaga et aLM^Olll 
M _^ „ l * L 'A . 



Navarette et a l . in pr ep 
'Smolcic ot al. 720121 

Scott ot al. |200S 



Younger ct al. 


20071 |2009l 


Arctxaga ct al. 


jpit ' 


Scott ct al. |200S^ 


Younger ct al. 


2007 2009[ 


Arctxaga ct al. 


120111 


|Aretxaga et alT 


poir| 



'Smolcic ct al. 72012^ 

Arctxaga ct al. (2011) 
Navarette ct al 
Smolcic ct al, 



n prep 
(2012:i 



Bcrtoldi et 


al. 


2-m\ 


Arctxaga c 


al. 




Oil 


Scott ct al. 


120 


oi 




Arctxaga c 
Arctxaga c 


al. 
al. 


I2011I 

12011 


Bcrtoldi et 


al. 


20071 


Arctxaga ct al. 


12011 


Arctxaga ct al. 
Bcrtoldi et al. 
.Arctxaga ct al. 


12011 
20071 
12OII1 


|Aretxaga et al. 


jiOll 



850.13 
850.15 

850.18 
850.20 
850.22 



850.33 
850.35 
850.50 
850.57 
850.62 
850.63 



AzTEC .1095957.22+022729.3 1.1 5.8±1.0 

AzTEC-9 0.89 13.5±1.8 

AzTEC/C18 1.1 7.9±1.5 

COSBOl 1.2 6.2+0.9 

AzTEC/C7 1.1 8.9±1.1 

AzTEC/C114 1.1 3.7±0.9 

AzTEC .1100025.23+022608.0 1.1 1.9±0.6 

AzTEC/C30 1.1 5.5±1.1 

AzTEC .1100023.98+022950.0 1.1 2.6±0.7 

AzTEC .1100026.68+023128.1 1.1 2.8±0.8 

COSLA-50 0.87 5.6±1.6 

COSLA-50 1.3 

AzTEC/C33 1.1 5.3±1.1 

C0SB019 1.2 3.0±0.8 

COSLA-38 0.87 5.8±1.6 

COSLA-38 1.3 8.19±1.85 

AzTEC .1095950.69+022829.5 1.1 3.6±0.9 

AzTEC/C74 1.1 4.2+0.9 

AzTEC/C35 1.1 5.2+1.0 

AzTEC/C45 1.1 4.8±1.0 

AzTEC J100031. 06+022751. 5 1.1 2.7±0.8 

COSB036 1.2 5.7±1.3 

AzTEC/C71 1.1 4.3±1.1 

COSB029 1.2 3.1±1.0 

AzTEC/C162 1.1 3.1±1.2 



AzTEC 

SMA 

AzTEC 

MAMBO 

AzTEC 

AzTEC 

AzTEC 

AzTEC 

AzTEC 

AzTEC 

LABOCA 

PdBI 

AzTEC 

MAMBO 

LABOCA 

PdBI 

AzTEC 

AzTEC 

AzTEC 

AzTEC 

AzTEC 

MAMBO 

AzTEC 

MAMBO 

AzTEC 



ISeott ct al.'('2008'l 
Younger ct al. 12007 2009[ 
Arctxaga ct al. 120111 



Bcrtoldi ct al. 


(2007) 


Arctxaga ct a. 
Arctxaga ct a. 


. (201l| 

. (2'nTI] 


Scott ct al. (20081 

Arctxaga ct al. CzOll 
'-:f-r-,++ ^1 / '?n(iijl 



^Scott et al.|(|2O0 
Navarette et a l . in pr ep 
■Smolcic et al.l<2012( 

Arctxaga ct al. (2011) 
Bcrtoldi et al. 120071 
Navarette et a l . in pr ep 
ISmolcic et al. 1120121 

Scott ct al. I 2003 

Arctxaga ct al. C 20 1 1 
Arctxaga et al. a3 
Arctxaga ct al. fl^ 
Scott et al. I 2008 
Bcrtoldi et al. 
Arctxaga ct al. 
Bcrtoldi et E 
Arctxaga ct al. 



Table Notes. Counterparts of ScuB.'l-2 sources at 870Aim-1.2 mm as identified by LABOCA, PdBI, AzTEC, SMA and MAMBO with given references. 
Sources arc divided into 450/^m >3.6fr samples and 850^m >3.6fT {but not 450/im-detected) samples. None of the marginal 3< cr <3.6 450^m sources were 
reported as detections in the given references. Sources marked with ** correspond to multiple ScuBA-2 sources in the Table, i.e., 450.16, 450.28 and 850.05 

all correspond to COSB03/AzTEC/C6. 



3.5 Reliability of Counterpart Identifications 

Identifying the correct multiwavelength counterpart (s) for 
each submillimeter source is the most important, but also 
the most difficult aspect of characterising the submillimetre 
galaxy population. Recently, it has become clear that direct, 
far-infrared interferometry is the only guaranteed method of 
identifying counterparts correctly (e.g. recent results from 
SMA, PdBI and ALMA targeting 850/im-l.lmm sources; 

[20TT| [Barger et aL]|20T2l [Smolcic et ai:][20T2l 
2012 I. For 850/xm sources detected with a ~15" 



Wang et al. 



Karim et al. 



beamsize, this work suggests that best-guess counterpart 
matching to radio and/or 24/im counterparts will fail ~30% 
of the time. This provides a good estimate to the accuracy of 



the counterparts given in Table [5] for our 850/im sample, but 
what is the reliability of our 450/im source counterparts? 

Measuring the reliability of the 450/im source counter- 
parts also requires far-infrared interferometric observations. 
In their absence, we can speculate that the reduced beam- 
size with respect to 850/im would see a reduction in coun- 
terpart contamination proportional to the reduction in sky 
area searched for potential counterparts. The positional un- 
certainties of 450/im range from 1.0-2.4" at >3.6o- (with a 
median of ~2.l"), while 850/im positional uncertainties vary 
from 2-6" at >3.6(t (median of 3.6"). By rough comparison, 
the ratio in searchable sky area at 450/im to 850/im is «0.17. 
If 850/tm has failure rates of ~30%, this suggests the same 



© 0000 RAS, MNRAS 000, 000-000 



14 C. M. Casey et al. 




850/u,m 



Measured S/N 

Figure 8. The probability of correct optical counterpart iden- 
tification at 450/im (green) and 850/xm (magenta) as a function 
of detection signal-to-noise measured by our Monte Carlo simula- 
tions. The vertical dashed lines mark the 3.6cr detection threshold. 



•^03 



50. 



A50. 



450.: 



counterpart identification strategy (when applied to 450/im 
samples) should fail only ~5% of the time. 

While this thought experiment gives a good ballpark 
figure for counterpart correctness, we also use the results 
from our Monte Carlo simulations in section [3?T] to estimate 
counterpart correctness. We take each source realisation and 
map it to an arbitrary optical source in the COSMOS field. 
We then offset the centre of the measured 450/im source 
from its input according to the offsets shown in Figure [4] 
with arbitrary angle; we then perform nearest-neighbour op- 
tical matching to that offset position. Figure[8]illustrates the 
fraction of all simulated sources which recovered the correct 
input optical counterpart as a function of measured signal- 
to-noise. In other words, at 4(t, the nearest optical coun- 
terpart of a 450/im source has a 75% probability of being 
the correct counterpart (65% chance for an 850/im source), 
whereas at 5cr that increases to 90% (70% chance at 850/im). 

We can also use these individual source probabilities 
to compute the likely number of incorrect counterparts in 
our 450/tm source list given in Table [T] The likelihood of all 
78 nearest-neighbour counterparts being correct ~1%. The 
probability of 3% contamination is 46%, of 5% contamina- 
tion is 25%, and above 5% contamination is 30% (above 10% 
contamination is only 3% likely). When analysing the bulk 
properties of the 450% sample (i.e. redshift distribution, lu- 
minosities, SFRs, etc), this will only impact our results by 
3-5%. 



3.6 Multiplicity 

In confused submillimeter maps with large beamsizes, mul- 
tiple sources can contribute flux within the same beam, mis- 
leading the observer into thinking there is perhaps only 
one bright source instead of two or more less luminous 
sources. Recent work have confirmed this by following up 



bolometer-identified sources interferometrically ( Wang et al 
2011[|Barger et al.||2012| |Karim et"aI1|2012[ ). The improved 
resolution from interferometry makes it clear that, in a sig- 
nificant fraction of cases, bright SMGs are actually mul- 
tiples. Whether or not these sources are associated pairs 
on a collision course (as suggested by |Hayward et aL]|2012[ ) 
or just projected pairs is unclear, but knowing that such 



Figure 9. A 35x40" 450/im signal-to-noise map cutout of 450.16 
and 450.28. The 450/tm sources are circled and labelled in green, 
while C0SB03, the MAMBO-detected 1.2 mm source is shown in 
red and the 850/tm ScuBA-2 source, 850.05, is shown in magenta. 
The beamsizes at 450/tm, 850/tm, and 1.2mm are 7", 15", and 
11" respectively. 



submillimeter-multiples exist is quite important to the even- 
tual evolutionary interpretation of these galaxies. 

Although ScuBA-2 is a bolometer array and not an in- 
terferometer, the high resolution at 450/im can detect multi- 
ple sources when longer-wavelength observations only mea- 
sure one SMG. This phenomenon is found once in our sam- 
ple, surrounding sources 450.16 and 450.28; the system is 
illustrated in Figure |9] There is only one long- wavelength 
submillimeter source detected in this area: 850.05 in our 
sample, as well as C0SB03, detected at 1.2 mm ([Bertoldi 
et al.|2007[), and AzTEC/C6, detected at 1.1 mm ([Aretxaga 



et al. 20111. However, we find two 450/im sources within 
the ScuBA-2 850/im beam: 450.16 with a S/N of ~4.6 and 
450.28 with a S/N of ~4.2. Although this is the only un- 
equivocal identification of a multiple at 450/im in our sam- 
ple, there is the potential for many more, at slightly lower 
signal-to-noise. When pushed deeper, closer to the confusion 
limit, it is likely that more 850/im-1.2mm sources will be 
shown to break into multiples at 450/im. This work demon- 
strates that shorter-wavelength, smaller beamsize bolome- 
ters can be used to probe submillimeter galaxy multiplicity 
in blank-field surveys at potentially lower observational cost 
than interferometric follow-up. 



4 DERIVED SOURCE CHARACTERISTICS 

This section explores the inferred physical characteristics of 
the ScuBA-2 detected galaxies and assumes that the best- 
guess counterparts identified in § |3.4| are correct. We have 
shown that contamination in this sample is minimal (on the 
3-5% level for both spurious identifications and incorrect 
counterpart identifications) , so the conclusions we draw here 
on the population should be widely applicable without sig- 
nificant bias. 
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Figure 10. The distribution in photometric redshift for ScuBA-2 
sources identified at 450/im (green) and 850/im {magenta). The 
median redshifts of the two samples are {z) = 2.0 and {z) = 2.2 
respectively. At 450/im, we illustrate two distributions, the first 
consisting of just >3.6fT 450/im sources (solid green shaded re- 
gion), and the second of the 3< a <3.6 sample (dashed green 
line) . We also indicate the number of sources without photomet- 
ric redshifts at right. Due to the larger beamsize at 850^m, our 
counterpart matching is done first using 24^m or radio identifica- 
tion (solid magenta shaded region). The remaining 850^m SMGs 
without 24^m/radio counterparts are shown in the hashed region. 
The dashed magenta line corresponds to the 3< a <3.6 sample. 



4.1 Redshift Distributions 

Figure [lO] shows the distributions of optical/near-infrared 
photometric redshifts ( Ilbert et al. | |20l "0l for the best-guess 
counterparts identified in Tables [7| and ^ Photometric red- 
shifts span < 2 < 5, with most sources at z < 3. The 
median redshift for the 450/im sample isz = 2.0±l.l while 
the median for the 850/im population is z = 2.2 ± 1.6, in 
agreement with previous work at 850pim ( [Chapman et al.] 



2003[ |2005[ ). It is expected that the 450/im-selected popu- 
lation peaks at lower redshift than the 850/xm population, 
although by how much is uncertain. [Geach et al.| ( |2012[ ), 
who reference Roseboom et al. in prep, find a 450/im redshift 
distribution which peaks at (z) = 1.3, which is statistically 
consistent with, but offset from, our finding of (z) = 2.0. 
This offset is likely due to the difference in depths of cov- 
erage between our surveys; the Geach et al. and Roseboom 
etal. work is a factor of ~3x deeper at its centre, while 
ours is a factor of Ri4x larger, so our sample naturally picks 
up more luminous sources that might sit at higher redshifts 
than the fainter 450^m-population. One might think that, 
if the fainter 450/im sources are less-luminous (and as a 
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Figure 11. The FIR colour properties of the ScuBA-2 sam- 
ple, measured as log(S'85o/S'45o). At top, we plot against the 
optical/near- infrared photometric redshift. The 450/im sample 
(green) has, on average, 450/tm flux densities which exceed 850/tm 
flux densities by a factor of 5, while the 850/xm sample (magenta) 
has naturally 'redder' colours. No redshift evolution is measured 
in either 450/im or 850/tm sample. At bottom, we plot the dis- 
tribution in FIR colour for both samples with the subset of the 
samples which do not have any photometric redshift. The distri- 
bution of colours for the samples without redshifts resembles the 
parent population, indicating no obvious redshift-colour bias. 



consequence at lower redshift), then perhaps including our 
marginally detected sample at 3 < a < 3.6 would shift our 
distribution towards lower redshifts. This isn't the 
can be seen by the dashed-line distributions on Figure |10[ 
and is likely because 850/im-priors were required for the 
marginal 450/im sample, thus swaying that low-luminosity 
end of the 450/im distribution towards higher redshifts. 

At 850/tm, we split our sample by its potential biases. 
The increased beamsize of 850/im observations means that 
850/im sources which lack 24/im or radio identifications have 
more uncertain counterparts. We illustrate the distribution 
of redshifts for the 24/im/radio identified subset, then those 
without (dashed area) and finally those from the 3 < a < 3.6 
marginal catalogue. Note that the number of sources lacking 
redshift identifications is ~17-18 (amounting to 17-23%) of 
either sample. 
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Figure 12. Example far-infrared SEDs for select ScuBA-2 sources, some 450/im-selected and some 850/im-selected. The assumed 
optical/near-infrared photometric redshift is specified in the upper right corner of each panel, as well as the measured Lir, Tjiusti and 
Afdust- If the SED fit used a fixed dust temperature due to a lack of >2 infrared photometric detections, then no dust mass is measured. 
In cases where there is no 24/jm detection, the mid-infrared portion of the SED is shaded to indicate its uncertainty. The uncertainty 
in LiR reflects this uncertainty. The radio portion of the SEDs are extrapolated from assuming the FIR/radio correlation holds (and 
evolves as measured inllvison et al.|2010a|b|l. 



Are those sources without photometric redshifts intrin- 
sically different than those with photometric redshifts? We 
investigate this by seeing if the FIR properties, i.e. Ssso to 
5*450 colour, (a) evolves with redshift, or (b) differs between 
the sample of galaxies with photometric redshifts and with- 
out. Figure [TT] plots the FIR colour of both >3.6(7 samples 
against photometric redshift and in histogram form. As ex- 
pected, 450/^m sources are much brighter at 450/^m than 
850/im and vice versa. What is interesting is that there does 
not seem to be any redshift evolution of the log(585o/S'45o) 
colour; this does not necessarily bode well for work which 
makes use of FIR photometric redshifts which assume a 
certain FIR template with fixed dust temperature and use 
FIR photometry to estimate the redshift, as we find that 



sources display a wide range of FIR properties irrespective 
of redshift. Furthermore, sources which do not have photo- 
metric redshifts show roughly the same distribution in FIR 
colour as the parent population of galaxies they are selected 
from. This implies that the systems lacking photometric red- 
shifts are not intrinsically biased with FIR properties, and 
by proxy, probably not biased with redshift. 

Although we did not find a evolution of FIR colour with 
redshift, one might think that the 850/im sources not de- 
tected at 450/im sit at higher redshifts on average. This is 
not the case. Sources with or without 450/.im detections are 
found to peak at 2; = 2.0; in fact, it is only after the inclusion 
of several 850/im sources detected at 450/im that the median 
redshift of the 850^tm sample shifts redward to {z) = 2.2. 
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4.2 Infrared luminosities, dust temperatures, dust 
masses 

We can go beyond basic FIR colour comparisons with 
our ScuBA-2 data and use SED fits to measure infrared 
luminosities, dust temperatures and dust masses. Simple 
greybody spectral energy distribution (SED) fitting com- 
bined with a mid-infrared powerlaw gives us this informa- 
tion with a minimum number of free parameters. We fit 
all available infrared photometry, from 24/.im MIPS, 70/im 
MIPS, 450^m ScuBA-2, SSO^m Scuba-2, 870^m LABOCA, 
1.1mm AzTEC, 1.2 mm MAMBO, and 1.3 mm PdBI where 
available. 

The combination of a greybody and mid-infrared pow- 
erlaw accounts for both galaxy-wide cold dust emission and 
smaller-scale warm dust emission. The fitting method is de- 



scribed in full in Casey (20121 and is represented by the 
equation 



5. = iVb 



^hcv/kT _ ^ 



+ iVpi(c/i.)"e 



(3) 



where Sv is flux density (in units of mjy), T is greybody dust 
temperature, /3 is the emissivity index of the greybody, a is 
the mid-infrared powerlaw slope, and A^'bb and A^pi are the 
normalisations of the greybody and powerlaw, respectively 
(A^pi is a fixed function of A^bb, T, and a). Note that the 
dust temperature from the complete SED, Tdust, is diflerent 
from the dust temperature measured from the greybody, T; 
Tdust is measured from the SED's peak in Si, via Wien's 
law, and in this case is ~7-10 K cooler than the input dust 
temperature T (see |Casey|2012[ for more details). Since the 
number of far-infrared photometric data points is limited to 
2-5 (from Spitzer-MYPS, ScuBA-2, AzTEC and MAMBO 
measurements), we only fit SEDs with two free parameters, 
Lib,, and Tdust and do not attempt to constrain emissivity or 
mid-infrared powerlaw index. The emissivity index, /3 = 1.5 
is fixed (e.g . |Pope et al.j|2008l [Younger et"ar][2009l |Chapin| 
et al ||2010[), as is the mid-infrared powerlaw index, a = 2.0 
(e.g. Blain et al.|2003 Casey|2012 I. A representative sample 
of best-fit SEDs, including both high and low signal-to-noise 
sources, is illustrated in Figure [12] 

Tables [To[ and [TT] give the basic derived properties from 
the SED fits, including Tir, the implied SFR from Tir, to 
one significant figure, and the dust temperatures and masses, 
Tdust and Mdust- We do not attempt to constrain sources 
without photometric redshifts, but there is no reason to be- 
lieve the FIR SEDs of those sources are intrinsically different 
than those in this sample. Sources with just two FIR pho- 
tometric points (and upper limits at other bands) are con- 
strained to have the mean dust temperature of their parent 
sample, either 450/im-selected galaxies or 850pim-selected 
galaxies. In the next section we measure these temperatures 
as 27.0±4.4K and 27±13K, respectively. 



4.3 Bulk Infrared Properties of the Population 

Figure[T3]illustrates stacked FIR photometric measurements 
of ScuBA-2 samples: the overlapping 450/im- and 850/im- 
selected sample, the entire 450/xm-selected and the 850/im- 
selected samples. Photometric data points are shifted in 
wavelength by photometric redshift and in flux density by 
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Figure 13. Stacked FIR photometric measurements for three 
subsamples of ScuBA-2 sources: the overlapping 450/jm and 
850/im identified sample (top panel), the entire 450^m selected 
sample (middle panel), and the entire 850/im selected sample 
(bottom panel). Sources are red- or blue-shifted to a redshift of 
z = 2, the approximate mean of the whole sample. We take me- 
dian photometric measurements to construct a more fully sampled 
SED representative of each population and fit it to a 4-parameter 
SED as described in the text. We do not measure significant dif- 
ferences in the dust temperatures or infrared luminosities of the 
subsamples. From this, we determine that sources without con- 
straining photometry should be fit with dust temperature fixed 
to =27 K, the mean of both parent 450/^m and 850/^m samples. 



the implied luminosity distance; all sources are shifted to a 
redshift of z — 2. Although there is substantial variation in 
the photometry caused by variation in dust temperature and 
luminosity, we can measure median photometric points to 
construct a more thoroughly sampled FIR SED. We then fit 
FIR SEDs to these constructed SEDs, as done for individual 
sources in the previous section. This SED fitting reveals that 
there are not pronounced differences between subsets of the 
ScuBA-2 population. The uncertainties on measured FIR lu- 
minosity and dust temperatures make all subsets of the data 
consistent with one another. Although this stacking gives the 
advantage of 'filling in' a poorly sampled SED, it washes out 
any variation which we might observe within samples which 
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Figure 14. Histograms of infrared luminosity and dust temper- 
ature (the two free parameters of our SED fitting teclinique) for 
various subsamples of ScuBA-2 sources: tlie overlapping >3.6(t 
450 and 850 sources (detected in both bands, blue), 450/im >3.6(t 
sources not detected at 850/^m (green), 850/xm >3.6o" sources 
not detected at 450/im (magenta), and the marginal 3< cr <3.6 
450 and 850 source catalogue (light blue). Vertical lines mark 
the median value for each subset. This illustrates that those sys- 
tems detected in both bands are the warmest galaxies (although 
the marginal catalogue sits at lower overall luminosities), while 
sources detected in only one band are less luminous and colder. 
Although the 450^m sources seem to be the coldest subsample 
here, we note that there are several sources found at >30 K, which 
is not seen in the 850/^m subsample. 



we can probe by comparing measured infrared luminosities 
to dust temperatures, dust masses and redshift. 

Figure [M] plots the distributions in infrared luminosity 
and dust temperature (the two free parameters of our SED 
fitting technique) to assess basic differences between subsets 
of the population. The subset of Scuba-2 sources detected at 
both 450/^m and 850^m above >3.6(t are the most luminous 
and warmest subset. The marginal 3< a <3.6 catalogue of 
sources sit at lower luminosities (by ~0.3dex) than those at 
>3.6cr, but equally warm dust temperatures. The subset of 
sources which are only detected at either 450/im or 850/im 
exhibit remarkably similar properties, despite some being 
only 450/im-detected and some only 850/im-detected. Both 
subsamples have lower infrared luminosities by ~0.2 dex and 
dust temperatures by ~3 K than those which are detected 
in both bands. Although 450/xm-detected galaxies might be 
thought to sit at warmer dust temperatures than galaxies 
detected at 850/im, we note that here the effects of dust 
temperature and luminosity are intertwined. There is an ob- 
served relation between dust temperature and infrared lu- 
minosity, so naturally, the samples at lower luminosities are 



cooler. Although the 450/im subsample has a cooler median 
dust temperature than 850/im sources, there are several out- 
liers at warm temperatures (>30K) that are not present in 
the 850/im subsample. 

Figure[l5]illustrates infrared luminosities against photo- 
metric redshift for both 450/tm- and 850/im-selected galax- 
ies, as described in Tables [lO] and |11| Both 450/tm and 
850/im samples by and large sit at 1 < 2 < 4 with 
10^^ < LiR < 10^"^ L0. This is very similar to the lumi- 
nosity and redshift range probed by the original ScuBA 
850/im surveys; the distribution of the [Chapman et alT] 
(2005) spectroscopically-confirmed sample (with improved 
luminosity estimates using photometry from |Kovacs et al.| 



2006 1 is shown for comparison. We also compare to the re- 



cent spectroscopically-confirmed sample of Herschel-SPiRE 



sources analysed in Casey et al. (2012b I and Casey et al 



(2012a I; the SPIRE sources (selected at 250-500/im) peak in 
density at lower redshifts, z ~ 1, thus are primarily lower 
luminosity sources. However, the SPIRE sources at z ^ 2 are 
similar to ScuBA-selected SMGs, and overall, slightly more 
luminous than our ScuBA-2-selected galaxies. 

Understanding the distribution of sources in Lir, — z 
requires an understanding of the source selection, which de- 
pends on source luminosity as well as dust temperature of 
the greybody portion of the SED. Selection at 850/im is well 
known to benefit from a negative K-correction, by which 
the luminosity detection threshold is fairly constant across 
1 < z < 8. Selection at 450/im, on the other hand, does 
not have the same benefit, since it probes the peak of the 
SED (with typical temperatures 20-60 K) at much lower 
redshifts; therefore, its luminosity detection threshold is a 
much stronger function of redshift than 850/im-selection. We 
can see this manifested in the rough detection thresholds 
seen in Figure [15] at different assumed dust temperatures 
at the two selection wavelengths. While overall the 450/im- 
and 850/im-sources occupy the same space, there are fewer 
450/im-selected ~10^^ Lq sources at 3 < « < 4 than 850/im- 
selected sources. 

The middle panel of Figure [15] illustrates Lir against 
dust temperature, Tdust, for sources with constrained tem- 
peratures. Here, the difference between Scuba/Scuba-2- 
selected sources against SPlRE-selected sources is more strik- 
ing. SPIRE sources sit at lower luminosities but compara- 
ble dust temperatures, which is primarily a function of the 
shorter selection wavelength (250-500/tm) than Scuba-2 
(450, 850/im). Scuba-2 detection is biased against warmer- 
dust SEDs, when sources of equal luminosity are considered; 
this is more pronounced at longer wavelengths. This poses a 
problem for ScuBA-2-selected samples since they are missing 
a potentially significant fraction of the infrared-based star 
formation activity in this redshift range, 1 < z < 4. This has 
already been demonstrated with _ffersc/ie/-SPIRE selected 



samples in Casey et al. ( 2012b| ) who find that a large fraction 
(~60-80%) of 250-500/im galaxies would not be detectable 
at 850/im. This leads us to ask how ScuBA-2-selected sources 
fit into the larger parent population of infrared-luminous 
galaxies at high-redshift? What can we learn from this sub- 
sample that is different from other submillimeter-selected 
samples? 

A different way of thinking about infrared-selected 
galaxies is by dust mass and not infrared luminosity. While 
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Figure 15. LEFT: Infrared luminosity against photometric redshift for 450/im (green) and 850/im-selected (magenta) sources. Detection 
limits at both wavelengths are estimated assuming SEDs of set dust temperatures, determined by the innerquartile and median dust 
temperatures measured in our sample (23— 35K, (T^ust) = 27 K). Sources with fewer FIR photometric measurements have much larger 
uncertainties on Lia. MIDDLE: Infrared luminosity against fitted SED dust temperature. Detection limits are presented here at fixed 
redshifts z = 1.2 (solid lines) and z = 3.0 (dashed lines), representing the lower limit and last quartile limit of our sample. Large circular 
points represent the median values of the similarly-coloured subsamples in Figure |14| RIGHT: Dust mass against dust temperature. 
Luminosity limits exclude warm dust sources with high dust masses. Besides the lowest redshift sources at low dust temperatures and 
dust masses, both 450/im and 850/im sources occupy similar parameter space, suggesting no strong evolutionary distinction in dust 
properties. Boundaries represent lower limits on dust mass at z = 1.2 (solid lines) and z = 3.0 (dashed lines). Herschet-SPIRE sources 
from [Casey et aL] (|2012b|a^ are overplottcd in all three panels for comparison. 



the latter probes the instantaneous star formation rate, the 
former is an indication of the existing reservoir of dust and 
should effectively trace the stellar mass, or the integrated 
star formation rate over a longer period of time than in- 
frared luminosity. We estimate dust masses for those SCUBA- 
2 galaxies which can be constrained well when redshifts and 
multiple photometric points exist. The right panel of Fig- 
ure [15] illustrates dust mass against dust temperature for 
the ScuBA-2 sample, as well as the Casey et al. ( 2012b|a" l 
Herschel-SPIRE selected sample. Scuba-2 sources are typi- 
cally more massive than SPlRE-selected sources on account 
of their selection wavelength (dust mass scales strongly with 
flux density at long wavelengths). So although the Spire 
sources often lie at higher integrated infrared luminosities 
than ScuBA-2 sources, they are less massive and perhaps at 
a stage of evolution which is statistically earlier than that 
of ScuBA-2 sources. 



5 DISCUSSION 

The ultimate goal of this work is the characterisation of 
galaxies which emit at 450/im and 850/tm, how they re- 
late to other similarly-selected galaxy populations, how they 
relate to the more extensively-studied optical/UV-selected 
populations, and how important they are in the context 
of total star formation in the Universe. After initial as- 
sessment of their universal 'importance,' follow-up studies 
can pursue source characterisation in more detail— from op- 
tical and near-infrared spectroscopic observations, resolved 
interferometric observations of dust, molecular gas and radio 
synchrotron emission, to kinematic characterisation. These 
initial observations from ScuBA-2 are also incredibly valu- 
able in providing predictions for observations from the next. 



large-scale ground-based infrared observatory, the Cornell- 
Caltech Atacama Telescope (CCAT). 



5.1 Distinction between 450/im and 

850/im— 1.2mm Selected Populations 

Many submillimeter-selected galaxy populations have been 
presented in the literature in the past decade, selected and 
characterised at a range of far-infrared wavelengths with 
different detection thresholds and confusion limits. This 
includes (but is not limited to) the original 850/im sur- 
veys from ScuBA ( [Hughes et al.[|l998| [Barger et ai.||1998| 
Small et al. 20021, 1.1 mm-selected AzTEC sources (Scott 
et al.|2008[|Aretxaga et al.|2011[), 1.2mm-selected MAMBO 
sources ( Bertoldl et al.[[2007 Lindner et al.|[2bll l, observa- 
tions of various sources at 350/im using sharc-II at the 



Caltech SubmlUlmeter Observatory (Kovacs et al 



250-500/im-selected blast sources (Dunlop et al. 



20061, 



20101, 



and 70-500//,m-selected Herschel sources. What do SCUBA- 
2 450/im-selected samples teach us that we could not learn 
from studying these other SMG populations? 

The primary distinction between the ScuBA-2 450/im- 
selected population and all other SMG populations is the 
accuracy with which multiwavelength counterparts can be 
identified without the need for interferometric follow-up. As 
discussed in the previous section, counterpart identification 
is difficult and often done incorrectly using rough scaling 
relations between 24/tm or 1.4 GHz flux density and the far- 
infrared. Identifying galaxies at 450/im within a beamsize 
~7" allows nearly unequivocal identification of optical/near- 
infrared counterparts, without bias towards mid-infrared or 
radio samples. 

Only 40% ScuBA-2 450/im-identified sources have 24/tm 
or 1.4 GHz counterparts, while 60-70% of 850/im sources 
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have counterparts. This means that work on submilhmeter 
galaxies which rely heavily on the existence of 24/im and/or 
radio counterparts should be taken with a grain of salt as 
selection biases could significantly impact the results of the 
study. If one assumes that the 450/im and 850/im popu- 
lations are intrinsically similar physically, then this could 
mean that a significant fraction of current 850/xm popula- 
tion analysis is based on misidentified counterparts. 

When contrasting the physical characteristics of 450/im 
with 850/im galaxies, like redshift distribution, luminosity, 
dust temperature and mass, this work suggests that the 
two populations are quite similar, perhaps surprisingly. The 
peak of the redshift distributions differs by Az «0.2 (for 
those which have optical/near-IR photometric redshifts). 
The luminosities of both populations ranges from ~lx 10^^- 
7x 10^^ L© and the samples overlap significantly (40% of the 
450/im sample is 850/im-detected and vice versa), yet the 
overlap is not complete. 

Dust temperature has been thought to cause a signif- 
icant selection difference between 850/im-lmm and popu- 
lations and those selected at shorter wavelengths, due to 
the strong dependence of ~850/im flu x density on dust t em- 
2, 



Blain et al. 



20041. Is 



perature (at z ~ 2, Ssso ~ LmT^ 
dust temperature the primary difference between the 450/im 
galaxies not luminous at 850/im and those that are? Are 
450/im galaxies warmer than their 850/im-selected cousins? 

Although the median dust temperatures of both 450/im- 
and 850/im-selected samples is the same (~27K), dust tem- 
perature does vary if we consider systems detected only at 
450/im or 850/im. As Figure \T4\ shows, systems detected in 
only one band or the other are colder than galaxies detected 
in both bands. While it might seem puzzling that 450/im 
galaxies are colder than systems detected in both bands, 
the selection effect has to do with luminosity. Since there's a 
known correlation between galaxy luminosity and dust tem- 
perature, those galaxies detected in both bands are more 
luminous, thus warmer. The 850/im-only detected systems 
are expected to be cooler and less luminous; the 450/im-only 
detected systems are mostly cooler (due to their lower lumi- 
nosities) although a few systems are hotter than the median 
of all subsamples >30K. 

On the surface, this observation— that lower luminosity, 
cooler galaxies are less likely to be 850/im detected— sounds 
contrary to the missing 'warm-dust' ULIRG hypothesis 



(Chapman et al. 2004 Casey et al. 20091. However, in re- 



ality the situation is more complex. The systems detected 
either only at 450/im or only at 850/im seem to be the same 
class of galaxy, but with slight differences in SED causing 
detection in one band and not the other. Besides the hand- 
ful of systems with warm 450/im temperatures, there's no 
distinction between the dust temperatures of the two sub- 
samples. It is likely that these systems sit near the detection 
limit in the opposite band and more observations, particu- 
larly at 450/im, would shed more light on the nature of the 
two populations. 

Due to its similarities to the 850/im population, the 
450/im population can be used as a probe of the entire 
850/im-detected population, since 450/im does not suffer 
from the same 24/im/radio counterpart-matching bias that 
our 850/im data does (and that most ~lmm data does, 
until full operation of AzTEC at the LMT or until con- 



struction of CCAT). Deeper observations and follow up of 
450/im-detected galaxies will provide essential clues to the 
underlying population biases in current, detailed SMG stud- 
ies. Furthermore, it is important to recognise that the lack 
of one-to-one overlap in the 450/tm and 850/im-samples im- 
plies that systems selected at any one wavelength in the FIR 
are not representative of all ultraluminous infrared star for- 
mation at high redshift. This becomes an issue when studies 
try to compare the relative importance of normal galaxies to 
infrared galaxies in the context of universal star formation. 



5.2 Relation to 'Normal' Galaxies 

Recent work has proposed that most submilhmeter or 
infrared- luminous systems are a different, more extreme 
class of galaxy than most 'normal' star forming galaxies in 
the Universe (e. g. |Daddi et al.|2009| |Genzel et al.|2010| |Tac- 
|coni et al.|2010[ ). This distinction is likely valid when speak- 
ing of the most luminous subset of SMGs. SMGs at «10^^ Lq 
are nearly all found to be major merger starbursts, much 
more extreme and much more rare than normal, more mod- 



est star forming galaxies at comparable redshifts (e.g. Engel 



et al. 2010 1. To place these starbursts in context with normal 
galaxies, many recent works have raised the argument that 
starbursts nearly always fall off (above) the galaxy "main se- 



quence," or at high specific star formation rates (see 


Noeske 


et al.l2007 


IDaddi et al.|2007||Elbaz et al.|2011||Roc 


ighiero 


et al.poil 


1. This is based on the observation that infrared 



luminous systems have extreme star formation rates yet ex- 
hibit stellar masses of typical star forming galaxies, so it 
is thought they are undergoing some rare star formation 
event that is atypical of the normal galaxy population. In- 
deed, most 'normal' galaxies which are not infrared lumi- 
nous have much more modest star formation rates, typically 
at the LIRG or sub-LIRG level ( < 150 Mq yr'^). 

Although many SMGs exhibit extreme properties not 
seen in typical star forming systems, this is not true of 
the entire infrared-luminous population. Here we measure 
a mean infrared luminosity of ScuBA-2 galaxies as «1- 
2xlO'^^L0, and the average star formation rates in our 



sample are 
(>101°Mq) 



:i200M 



©yr 



which is typical of high mass 



16 



shows 
shown 



normal' galaxies at z ~ 2. Figure 
where ScUBA-2 galaxies sit on the "main sequence 
here as star formation rate against stellar mass. COSMOS 
galaxies with comparable photometric redshifts follow a se- 
quence where stellar mass traces star formation rate. Both 
quantities are lifted from the Ilbert et al. catalogue and are 
estimated with the Le Phare code which makes use of 



the Bruzual & Chariot (20031 templates. When consider- 



ing the measured star formation rates from our infrared 
SEDs, several of these Scuba-2 sources sit near the main 
sequence, particularly at the high-mass end, but overall we 
see that most SMGs do fall above the relation. However, is 
the comparison of infrared-derived star formation rates to 
UV/optical/near-IR derived star formation rates a fair one? 

Tables [To] and [TT] also give the Le PHARE-estimated 
quantities for star formation rate and stellar mass, as de- 
rived from the UV through to the near-infrared. The right 
panel of Figure [16] compares the star formation rates from 
the UV/optical/near-infrared (dubbed "SFRuv") directly 
with the star formation rates in the far-infrared. Despite the 
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Figure 16. LEFT: Stellar mass against star formation rate for galaxies in the COSMOS field (contours), with our ScuBA-2 sources 
ovorplottcd as blue points. Two different measurements of star formation rate are given and connected for each source by a gray line: 
SFRiR (light blue asterisks) from our far-infrared SED fitting, and SFRuv (navy blue points) from Ilbert ct al. (2010). RIGHT: The 
infrared-derived star formation rate, SFRjj^ against the ratio of SFRjj^ to SFRyy. Green points are sources detected at 450/xm while 
magenta points here are sources detected at 850/im but not 450/im. On average, the star formation rates calculated from templates in the 
UV/optical/near-infrared underestimate the star formation rate of these SMGs by ~13 times, but could be underestimates by factors up 
to ~100 times. Also consider that SFRuv is also already supposed to be corrected for dust extinction. The left panel "main sequence" 
plot also illustrates this severe discrepancy. 



fact that SFRuv is supposedly corrected for the effects of 
dust extinction, there is a huge discrepancy between SFRuv 
and SFRiR. On average, SFRuv underestimates the star for- 
mation rate (which can be bolometrically approximated as 
SFRiR in infrared-luminous galaxies) by a factor of 10. 
When placed on the "main sequence" with SFRuv instead of 
SFRiR, the ScuBA-2 sources sit directly in the heart of the 
main sequence. While these ScuBA-2 sources are known to 
be unusual in that the are directly detected in the infrared, 
and therefore are the most extreme subset, the discrepancy 
between star formation rate measurements should be alarm- 
ing, especially if galaxies' placement on the "main sequence" 
is to be interpreted as an evolutionary trend. More work on 
understanding this discrepancy, especially at lower infrared 
luminosities (which can be measured through stacking), is 
desperately needed, especially with so many extragalactic 
evolutionary studies making use of this galaxy "main se- 
quence." 



others, all with typical statistical uncertainties of ~0.3dex 
( Borys et a"L]|2005| [Hainline et al.||2009[ [Michalowski et al 
20121. Higher stellar masses could mean that SMGs sit on 



the main sequence even when using their infrared-derived 
star formation rates. 

The magnitude of these discrepancies (~ldex in SFR 
and at least ~0.3dex in M*) illustrates that the placement 
of SMGs on or off the "main sequence" is not so straightfor- 
ward. While it might be tempting to generalise evolutionary 
scenarios based on specific star formation rate, this com- 
parison highlights that there is no overwhelming reason to 
think that most SMGs are such outliers, besides being very 
dusty. Whether most ^lO^'^ Lq SMGs are merger dominated 
or secularly evolving galaxies is unclear (since statistics on 
evidence either way has been too small), although break- 
ing down their evolutionary histories will be essential, given 
their relatively significant contribution to cosmic star for- 
mation, which we discuss in the next subsection. 



The discrepancy is not only present in star formation 
rate measures, but also stellar mass. Although we leave a rig- 
orous measurement of stellar mass in these Scuba-2 galaxies 
to a future work, it is important to highlight that the mea- 
surement of stellar mass is highly dependent on an accurate 
understanding of star formation history. The stellar masses 



calculated via Le Phare and the Bruzual & Chariot ( 2003 1 



templates includes a range of star formation histories which 
does not accommodate for extreme evolutionary scenarios 
like major mergers. Even when the probable merging his- 
tory of ultraluminous-infrared galaxies is taken into account, 
there are still discrepancies in the literature. Some SMG stel- 
lar mass estimates are systematically lower or higher than 



5.3 Contribution to Universal Star Formation 

Comparing populations can be done in terms of their net 
contribution to the star formation rate density (SFRD) 
of the Universe, or in other words, how significantly they 
impact net star formation in the Universe at any given 
epoch. We compute SFRD estimates for the ScuBA-2 pop- 
ulation in this paper by assuming uniform field depth over 
the 394arcmin^ survey area. The maximum accessible vol- 
ume for a source is computed given that sources' selection 
wavelength, flux density, and dust temperature (either mea- 
sured explicitly or assumed to be 27 K if unconstrained). 
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We split the sample into five redshift bins spanning the 
range of our data: 0<2<1, 1<2;<2, 2<2;<3, 
3 < z < 4, and 4 < z < 5. Figure [iT] shows the results; 
we see that overall, Scuba-2 sources follow expectations 
of 850/im-selected samples of similar depths, contributing 
~0.01 Mq yr"^ Mpc~^ across < z < 5. When split up 
into 450/im or 850/xm samples, the agreement with the lit- 



erature sources selected at 850/im-1.2mm 


Chapman et al. 


2005 Wardlow et al.|2011 Roseboom et al. 


2012 1 improves. 



This is likely because the combined 450/im- and 850/im- 
ScuBA-2 sample is a more complete subset of 10^^ Lq 
activity at these epochs than populations just selected at 
850/im-1.2mm. 

The contribution from these infrared-luminous galaxies 
at z > 2 does seem to be significantly lower than the SFRD 
contribution from more normal, optical and UV selected 



galaxies (Hopkins & Beacom 20061. However, it should be 



noted that the ScuBA-2 contribution noted here is only rep- 
resentative of the galaxies detected and analyzed in this 
paper. When deeper and wider-area coverage is available 
in the future, and when there is higher confidence in the 
sources' redshifts (either from spectroscopic identification 
or analysis of the quality of photometric redshifts in dusty 
galaxies), this should be revisited. With substantial depth 
and confident redshifts, luminosity functions can be con- 
structed and the contributions of infrared-based star for- 
mation can be estimated down to fixed luminosity limits 
(as in Le Floc'h et al. 2005 out to 2 ~ 1). Similarly, the 
contribution of UV/optically bright systems to the SFRD 
should be reassessed; given the huge discrepancy between 
predicted dust extinction and measured extinction in this 
sample, UV/optical estimates to SFRD should not incorpo- 
rate correction for dust. Then the SFR from infrared can be 
directly compared to SFR from optical, over a range of lu- 
minosities, and provide the best insight into the role of dust 
obscuration in cosmic star formation. 



5.4 Expectations from CCAT 

These ScuBA-2 observations represent one of the first for- 
ays into deep, uniform ~mjy mapping in ~200-500/im sub- 
millimeter regime, but many future observations will push 
deeper towards the confusion limit (nearing ~0.5mJy RMS 
at 450/im with the ScuBA-2 Cosmology Legacy Survey, 
S2CLS). This work has demonstrated the overall (and per- 
haps interesting or unexpected) similarity between 450/im 
>12mJy and 850/im >2mJy populations, and drawn com- 
parisons to populations selected at shorter wavelengths. 
What can we anticipate for the fainter 450/im sources, and 
for the future 200-800j(im sources detected by CCAT which 
will reach into the sub-mjy regime, without the additional 
aid of gravitational lensing? 

The chief improvement CCAT will bring with its 25 m 
dish is similar high-resolution (small beamsize) mapping. Al- 
though similar to the capabilities of LMT at ~1 mm, the res- 
olution CCAT will achieve this at shorter wavelengths (200- 
850/im) and will be at least a factor of two improved over 
JCMT. Situated in the Atacama near the ALMA site also 
means dryer conditions than Mauna Kea and more opportu- 
nity to map the sky at wavelengths where the sky is typically 
opaque. Although Scuba-2 is groundbreaking in opening 
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Figure 17. The estimated contribution of ScuBA-2 sources in 
this paper to the cosmic star formation rate density (SFRD). 
We split up the contribution estimates by selection wavelength, 
showing both the contribution from 450/im-selected (dark green) 
and 850/im-seIected sources (magenta), and we compare the to- 
tal (black points) to literature SFRD values for similarly selected 
populations. This includes the [Chapman et ah ( 2005 ' 850/jm- 



selected SCUBA sources (blue), [Wardlow et al".| ||2011| 870/im- 
selected LABOCA sources (light sea green), the [Roseboom et al.[ 
||2012[l 1.2m m-sclected MAMBO sources (light green), and the 
[Casey et al.||2012b|a| Herschel-SPIRE 250-500/^m-selected sources 
(orange). Note that our SFRD estimates are only representative 
of the sources summarised in this paper, and do not extend far- 
ther down the luminosity function. It also does not attempt to 
correct for incompleteness in the (U/Hy)LIRG selection biases in 
the submillimeter, so these estimates should be regarded as lower 
limits on the contribution from ultraluminous infrared activity to 
cosmic star formation. 



up our window to shorter wavelength far-infrared selection, 
CCAT will reveal sources orders of magnitude more abun- 
dant over much larger areas of sky. CCAT will be the tele- 
scope needed to link the studies of ultraluminous infrared 
activity to normal star-forming galaxies at high-redshift and 
solve issues like source multiplicity and accurately-sampled 
infrared star formation rates efficiently. 



6 CONCLUSIONS 

With simultaneous mapping at 450/im and 850/im, the 
Scuba-2 instrument makes substantial headway in the de- 
tection and characterisation of infrared emission in the dis- 
tant Universe. The ~7" high-resolution blank-field mapping 
at 450/tm opens up a unique parameter space in high-z 
submillimeter science, and the efficient mapping at 850/tm 
means that the confusion limit is reached in far less time 
than previous observations. This paper has presented initial 
Scuba-2 mapping of a uniformly-covered 394 arcmin'^ area 
in the COSMOS field at both 450/im and 850/im. We reach 
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the following conclusions about our maps and the identifi- 
cation of point sources: 

• The root-mean-square (RMS) noise values of our maps 
are (7450 =4.13mJy and asso =0.80 mjy at our map cen- 
tres and find largely uniform RMS over a 22.4' diameter 

(394 arcmin^), the threshold where the RMS exceeds two 
times the central RMS value. 

• We extract point sources down to a 3.6(t detection 
threshold at 450/im and 850/im, where we find 78 and 99 
sources respectively. These extraction limits are determined 
via Monte Carlo simulations and an expected sample con- 
tamination rates of 3-5%. 

• Source number counts at 450/L(m and 850/xm are mea- 
sured and we compare them to previous results at both 
wavelengths. This area of the COSMOS field exhibits a 
slight deficit at both wavelengths, which we attribute to 
cosmic variance. We provide best-fit double-power law and 
Schechter function fits. 

• Positional uncertainties of 1-2.5" at 450/im and 2-6" 
at 850/jm are estimated using Monte Carlo tests; we use 
this uncertainty to match both 450//m and 850/im sources 
to multiwavelength counterparts. The low positional uncer- 
tainties at 450/xm allow direct matching to counterparts in 
the optical/near-infrared, without relying on detection at 
24/im or radio wavelengths. Our 850/Lim counterpart match- 
ing is done first by identifying 24/Ltm and radio counterparts, 
then if absent, taken as the nearest neighbour optical source. 

• 60% of all 450/im galaxies and 40% of all 850/im galaxies 
lack both 24/im and radio counterparts down to the deep 
field detection limits (524 « 80/tJy and 51.4 «50/tJy). This 
suggests that submillimeter samples relying on identification 
at either wavelength could be significantly biased, producing 
~l/2 misidentifications. 

• We observe a case of a long-wavelength identi- 
fied (850/tm-1.2mm) SMG splitting into multiple 450/im 
sources. Although SMGs have been resolved as multiples 
with far-infrared/submm interferometry, this is the first 
time a SMG has been resolved into multiples using high- 
resolution bolometer observations. 

After analysis of field number counts and counterpart 
matching techniques, we also analyse the population char- 
acteristics of 450/tm and 850/tm galaxies, how they differ 
from one another and how they relate to other high-redshift 
galaxy populations. We reach the following conclusions on 
their physical attributes: 

• The redshift distributions for both populations is mea- 
sured using the extensive COSMOS ancillary optical and 
near-infrared data. The 450/im population peaks at (z) —2.0 
while the 850/im population peaks at (z) =2.2. 

• The far-infrared colour of ScuBA-2 sources (measured 
as log(S85o/S45o)) is found not to evolve with redshift, al- 
though 450/im identified galaxies have statistically 'bluer' 
far-infrared SEDs than 850/im identified galaxies. Similarly, 
850/tm-detected sources which are 450/im-dropouts peak at 
z m 2 and do not sit at higher redshifts as would be sus- 
pected if colour did evolve with redshift. 

• Infrared luminosities, dust temperatures and dust 
masses are estimated for both populations, and we deter- 
mine that 450/tm galaxies and 850/tm galaxies are physi- 



cally similar, with slight differences being the Az =0.2 dif- 
ference in redshift distribution. Galaxies detected in only one 
of the two bands are less luminous by ~0.2dex than those 
detected in both, although there is no significant difference 
between the subset of 450/im-detected/850/im-dropout and 
850/tm-detected/450/tm-dropout galaxies. This contradicts 
some previous work which anticipated that 450/im-identified 
850/tm-dropouts would be warmer, but equally luminous to 
their 850/(m counterparts. 

• The star formation rates measured directly in the 
far-infrared exceed those predicted from UV/optical/near- 
infrared photometry by a factor of ~13, on average. This 
implies that template-derived star formation characteristics 
should not be used for infrared-luminous samples and should 
likely be revisited for high-star formation rate systems that 
are perhaps just under current far-infrared detection limits. 

• Placing ScuBA-2 sources on a galaxies "main sequence" 
plot, our systems (450/im and 850/im) fall well above the 

galaxies main sequence, with very high specific star forma- 
tion rates. However, when using the star formation rates 
from the UV/optical/near-infrared, they lie on the relation. 
Notwithstanding the large uncertainty in stellar mass which 
we do not address in this paper, this work emphasises the 
need to proceed with caution when using a galaxy's place- 
ment on the main sequence alone as evolutionary interpre- 
tation. 

• The contribution of these ScuBA-2 sources to the cos- 
mic infrared luminosity density, or the star formation rate 
density, is measured. We find the contributions of 450/tm 
and 850/tm peaks at « « 2, at ~l/5 the total measured 
value at that epoch in the literature. Over < z < 5, 
we estimate the contribution from these sources is fiat at 
RiO.Ol Mq yr^^ h'* Mpc^"^. Note however, this only repre- 
sents the contributions from sources detected in this pa- 
per, and does not represent a fixed luminosity sample or 
is it representative of all ultraluminous infrared activity in 
the early Universe. We anticipate that the contribution of 
infrared-identified systems to cosmic star formation will be 
measured to increase when more complete, deeper samples 
are formed in the coming years with more sophisticated far- 
infrared observations. 



ScuBA-2 high-resolution 450/im mapping has allowed the 
first detailed look at an infrared-luminous population not 
surveyed before. This work has revealed the necessity to 
exercise caution when identifying rrmltiwavelength counter- 
parts of submillimeter sources (especially those with a large 
beamsize) and the biases and limitations of previously anal- 
ysed 850/im-lmm samples. Direct far-infrared interferom- 
etry still proves the best method for unequivocally iden- 
tifying multiwavelength counterparts, however future work 
from CCAT will enable large field-of-view mapping with res- 
olution significantly improved over JCMT, both at 850/im 
and 450/(m. This work provides a stepping stone to link 
infrared-luminous systems to the more ubiquitous, 'normal' 
star forming galaxies across a wide range of epochs, shed- 
ding light on galaxy evolution and the nature of cosmic star 
formation. 
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Figure 18. 12"xl2" cutouts for 450/im-detected sources in optical tri-colour (SUBARU B, I, and Z) and 5pjfaer-IRAC 3.6/im. Wiiite 
solid circles indicate the position of the 450/xm source with associated positional uncertainty, varying 1.5—3". The outter ring (white 
dashed circle) corresponds to the 95% confidence interval on position, while the inner is the 90% interval. Spitzer-WllPS 24fim source 
positions are marked with red circles and 1.4 GHz source positions are marked with green circles. Optical counterparts are marked with 
blue circles. 
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Table 5. 850/im-idcntificd >3.6cr point sources in COSMOS. 
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< 12.85 


SMM J100024 


0+022947. . . 


850 


18 


10 


00 


24 




02 


29 


47 


6.85 


5 


87±0 


86 


5 


13±1 


16 




45 


<0 


01 




6.37+4.36 


< 13.09 


SMM J100015 


7+022445. . . 


850 


19 


10 


00 


15 


7 


02 


24 


45 


6.74 


5 


40±0 


80 


4 


71±1 


08 




48 


<o 


01 




2 . 13 + 4.13 


< 12.39 


SMM J100027 


0+023137. . . 


850 


20 


10 


00 


27 





02 


31 


37 


6.64 


6 


69±1 


01 


5 


83±1 


36 




51 


<0 


01 




3.93 + 5.00 


< 15.00 


SMM J100035 


9+022151. . . 


850 


21 


10 


00 


35 


9 


02 


21 


61 


6.60 


5 


29±0 


80 


4 


60±1 


08 




52 


<o 


01 


450.58 


15 . 32 + 4. 12 


11.35 + 4.90 


SMM J100018 


7+021655. . . 


850 


22 


10 


00 


18 


7 


02 


16 


65 


6.25 


6 


05±0 


97 


5 


15±1 


28 




69 


<0 


01 




10.34 + 5. 10 


< 15.29 


SMM J100022 


1 +022842. . . 


850 


23 


10 


00 


22 


^ 


02 


28 


42 


6.25 


5 


06±0 


81 


4 


31±1 


07 




69 


<o 


01 


450.94 


14.67+4. 19 


10 . 14+ 5.17 


SMM J100004 


8+023045. . . 


850 


24 


10 


00 


04 


8 


02 


30 


45 


6.20 


6 


26±1 


01 


5 


32±1 


34 




69 


<0 


01 


450.20 


23 . 28 + 5.24 


19.11 + 5.91 


SMM J100012 


5+021444. . . 


850 


25 


10 


00 


12 


5 


02 


14 


44 


6.16 


7 


93±1 


29 


6 


72±1 


70 




69 


<o 


01 




14.85 + 6.90 


<20.71 


SMM JlOOOlO 


3+022627. . . 


850 


26 


10 


00 


10 


3 


02 


26 


27 


6.09 


4 


87±0 


80 


4 


12±1 


06 




69 


<0 


01 


450.32 


16.86 + 4. 1 1 


13.37+4.67 


SMM J100023 


5+021918. . . 


850 


27 


10 


00 


23 


5 


02 


19 


18 


6.09 


5 


00±0 


82 


4 


23±1 


08 




69 


<o 


01 


450.215 


13.05 + 4.20 


7.48 + 5.28 


SMM J100004 


3+022059. . . 


850 


28 


10 


00 


04 


3 


02 


20 


69 


5.79 


4 


86±0 


84 


4 


07±1 


08 




76 


<0 


01 


450.15 


19.92 + 4.33 


16.53 + 4.88 


SMM J100050 


1+022615. . . 


850 


29 


10 


00 


50 


^ 


02 


26 


15 


5.70 


5 


13±0 


90 


4 


28±1 


15 




79 


<0 


01 




10.89 + 4.50 


< 13.49 


SMM J100006 


7+022209. . . 


850 


30 


10 


00 


06 


7 


02 


22 


09 


5.69 


4 


58±0 


80 


3 


82±1 


03 




79 


<o 


01 




3.90 + 4. 15 


< 12.45 


SMM J100026 


1+021741. . . 


850 


31 


10 


00 


26 




02 


17 


41 


5.57 


5 


02±0 


90 


4 


19±1 


14 




84 


<0 


01 




1 2 . 34 + 4.62 


< 13.85 


SMM JlOOOOl 


5+022429. . . 


850 


32 


10 


00 


01 


5 


02 


24 


29 


5.56 


4 


56±0 


82 


3 


81±1 


04 




84 


<o 


01 


450.13 


19.96 + 4.20 


16.72 + 4.73 


SMM J095950 


7+022827. . . 


850 


33 


09 


59 


50 


7 


02 


28 


27 


5.55 


6 


20±1 


12 


5 


18±1 


41 




85 


<0 


01 


450.99 


20. 15 + 5.80 


13.76 + 7.20 


SMM J100016 


4+022638. . . 


850 


34 


10 


00 


16 




02 


26 


38 


5.49 


4 


41±0 


80 


3 


68±1 


01 




86 


<o 


01 


450.09 


20.41 + 4. 14 


17.26 + 4.68 


SMM J100105 


1+022150. . . 


850 


35 


10 


01 


05 




02 


21 


60 


5.47 


7 


60±1 


39 


6 


35±1 


74 




87 


<0 






0.38 + 6.54 


< 19.62 


SMM JlOOOOl 


3+021745. . . 


850 


36 


10 


00 


01 


3 


02 


17 


45 


5.39 


5 


60±1 


04 


4 


68±1 


30 




90 


<o 


01 




0.09 + 5.62 


< 16.86 


SMM J100035 


6+022826. . . 


850 


37 


10 


00 


35 


6 


02 


28 


26 


5.38 


4 


45±0 


83 


3 


72±1 


03 




91 


<0 


01 




1 . 74 + 4. 25 


< 12.74 


SMM J100023 


5+021536. . . 


850 


38 


10 


00 


23 


5 


02 


15 


36 


5.22 


5 


67±1 


09 


4 


72±1 


34 




99 


<o 


01 




13.39 + 5.74 


< 17.21 


SMM J100012 


1 +022310. . . 


850 


39 


10 


00 


12 




02 


23 


10 


5.21 


4 


17±0 


80 


3 


47±0 


99 




99 


<o 


01 




10.53 + 4. 13 




SMM J100013 


5+022227. . . 


850 


40 


10 


00 


13 




02 


22 


27 


5.20 


4 


16±0 


80 


3 


47±0 


99 




00 


<0 


01 


450. 252 


1 2 . 57+4. 14 


6.91 + 5.15 


SMM J100025 


4+022544. . . 


850 


41 


10 


00 


25 




02 


25 


44 


5.17 


4 


12±0 


80 


3 


43±0 


98 


3 


01 


<o 


01 


450.11 


19.57+4. 11 


16 . 40 + 4.62 


SMM J100017 


2+022519. . . 


850 


42 


10 


00 


17 




02 


25 


19 


5.14 


4 


12±0 


80 


3 


43±0 


99 


3 


02 


<0 


01 


450. 17 


19.01 + 4.14 


15.76 + 4.67 


SMM J100026 


3+021528. . . 


850 


43 


10 


00 


26 


3 


02 


15 


28 


5.02 


5 


55±1 


11 


4 


60±1 


35 


3 


10 


<o 


01 




2 . 32 + 5.82 


< 17.46 


SMM J100018 


2+022249. . . 


850 


44 


10 


00 


18 




02 


22 


49 


5.00 


3 


99±0 


80 


3 


31±0 


97 


3 


10 


<0 


01 




3 . 93 + 4. 1 2 


< 12.35 


SMM J100006 


9+022047. . . 


850 


45 


10 


00 


06 


9 


02 


20 


47 


4.99 


4 


12±0 


83 


3 


42±1 


01 


3 


11 


<o 


01 




3 . 79 + 4.25 


< 12.76 


SMM J100025 


4+021506. . . 


850 


46 


10 


00 


25 


^ 


02 


15 


06 


4.85 


5 


62±1 


16 


4 


64±1 


40 


3 


20 


<0 


01 


450.08 


31.25 + 6.06 


26.71 + 6.93 


SMM J095952 


3+022137. . . 


850 


47 


09 


59 


52 


3 


02 


21 


37 


4.83 


4 


74±0 


98 


3 


92±1 


19 


3 


21 


<o 


01 


450.50 


19.82 + 5.25 


14.90 + 6 . 19 


SMM J100037 


0+021942. . . 


850 


48 


10 


00 


37 





02 


19 


42 


4.83 


3 


99±0 


83 


3 


30±1 


00 


3 


22 


<0 


01 


450.42 


16. 71 + 4.24 


13.01 + 4.89 


SMM J100005 


0+021719. . . 


850 


49 


10 


00 


05 





02 


17 


19 


4.80 


4 


91±1 


02 


4 


05±1 


24 


3 




<o 


01 


450. 193 


17.40 + 5.50 


10.27+6.95 


SMM J100009 


1 +022023. . . 


850 


50 


10 


00 


09 




02 


20 


23 


4.78 


3 


94±0 


82 


3 


25±1 


00 


^ 


25 


<o 


01 




4.68 + 4.23 


< 12.69 


SMM J095956 


4+021854. . . 


850 


51 


09 


59 


56 




02 


18 


54 


4.74 


4 


94±1 


04 


4 


06±1 


26 


3 


27 


<0 


01 




9.29 + 5.61 


< 16.83 


SMM JlOOOll 


9+022937. . . 


850 


52 


10 


00 


11 


9 


02 


29 


37 


4.74 


4 


09±0 


86 


3 


37±1 


04 


3 


28 


<o 


01 


450.134 


14. 45 + 4.36 


9.26 + 5.50 


SMM JlOOOll 


4+021508. . . 


850 


53 


10 


00 


11 




02 


15 


08 


4.68 


5 


72±1 


22 


4 


70±1 


47 


3 


32 


<0 


01 


450.206 


20.69 + 6.60 


1 2 . 04 + 8.32 


SMM J095948 


9+022748. . . 


850 


54 


09 


59 


48 


9 


02 


27 


48 


4.56 


5 


10±1 


12 


4 


17±1 


35 


3 


41 


<o 


01 




0.47+5.90 


< 17.70 


SMM J095950 


9+022742. . . 


850 


55 


09 


59 


50 


9 


02 


27 


42 


4.52 


4 


77±1 


06 


3 


89±1 


27 


3 


43 


<0 


01 


450.86 


19.80 + 5.58 


13.95 + 6.80 


SMM J100005 


1 +021526. . . 


850 


56 


10 


00 


05 


^ 


02 


15 


26 


4.42 


5 


69±1 


29 


4 


63±1 


54 


3 


58 


0.02 




-4.70 + 6.84 


<20.51 


SMM J100006 


8+023307. . . 


850 


57 


10 


00 


06 


8 


02 


33 


07 


4.35 


5 


94±1 


36 


4 


S2±l 


63 


3 


68 


0.02 




3.66 + 7. 17 


<21 .50 


SMM JlOOlOl 


3+022800. . . 


850 


58 


10 


01 


01 


3 


02 


28 


00 


4.33 


5 


67±1 


31 


4 


60±1 


57 


3 


71 


0.03 


450.36 


24.27+6.07 


19.03 + 6.96 


SMM JIOOOOO 


7+022001. . . 


850 


59 


10 


00 


00 


7 


02 


20 


01 




3 


92±0 


91 


3 


18±1 


08 


3 




0.03 




2.88 + 4.76 


< 14.28 


ID IVl iVl J ± u u u u u 






















4 32 


3 


79±0 


88 


3 


07±1 


05 




73 


0.03 




11 26 + 4 49 




SMM J100002 


0+022820... 


850 


61 


10 


00 


02 





02 


28 


20 


4.31 


3 


83±0 


89 


3 


10±1 


06 


3 


73 


0.03 




11.06 + 4.53 


<13.59 


SMM J100031 


0+022751... 


850 


62 


10 


00 


31 





02 


27 


51 


4.29 


3 


46±0 


81 


2 


80±0 


96 


3 


76 


0.03 


450.78 


14.94 + 4.15 


10.72 + 4.99 


SMM J095953 


3+021850... 


850 


63 


09 


59 


53 


3 


02 


18 


60 


4.28 


4 


77±1 


11 


3 


86±1 


33 


3 


78 


0.03 




3.99 + 6.07 


<18.20 


SMM JIOOOOO 


6+022137... 


850 


64 


10 


00 


00 


6 


02 


21 


37 


4.18 


3 


59±0 


86 


2 


89±1 


02 


3 


95 


0.04 




4.27+4.41 


<13.24 


SMM J100059 


2+022108... 


850 


65 


10 


00 


59 


2 


02 


21 


08 


4.15 


4 


71±1 


14 


3 


79±1 


35 


4 


04 


0.04 




2.93 + 5.46 


<16.38 


SMM J100024 


4+022831... 


850 


66 


10 


00 


24 


4 


02 


28 


31 


4.09 


3 


30±0 


81 


2 


64±0 


96 


4 


19 


0.05 




8.36 + 4.17 


<12.52 


SMM J100032 


4+023321... 


850 


67 


10 


00 


32 


4 


02 


33 


21 


4.08 


5 


15±1 


26 


4 


12±1 


50 


4 


21 


0.05 


450.240 


19.01 + 6.21 


10.61 + 7.75 


SMM J095955 


3+021954... 


850 


68 


09 


59 


55 


3 


02 


19 


64 


4.06 


4 


05±1 


00 


3 


24±1 


19 


4 


27 


0.05 




11.46 + 5.35 


<16.05 


SMM J100043 


8+022857... 


850 


69 


10 


00 


43 


8 


02 


28 


67 


4.03 


3 


74±0 


93 


2 


98±1 


10 


4 


34 


0.06 




-2.18 + 4.64 


<13.91 


SMM J100020 


0+023020... 


850 


70 


10 


00 


20 





02 


30 


20 


4.02 


3 


57±0 


89 


2 


85±1 


06 


4 


36 


0.06 




6.36 + 4.49 


<13.47 


SMM J095944 


0+022105... 


850 


71 


09 


59 


44 





02 


21 


05 


4.01 


5 


07±1 


26 


4 


04±1 


50 


4 


38 


0.06 




6.21 + 6.92 


<20.76 


SMM J100013 


6+021731... 


850 


72 


10 


00 


13 


6 


02 


17 


31 


3.99 


3 


77±0 


95 


3 


00±1 


13 


4 


45 


0.06 




3.00 + 4.96 


<14.89 


SMM J100020 


6+022251... 


850 


73 


10 


00 


20 


6 


02 


22 


61 


3.89 


3 


11±0 


80 


2 


46±0 


95 


4 


69 


0.07 




-0.35 + 4.11 


<12.34 


SMM J100005 


5+021448... 


850 


74 


10 


00 


05 


5 


02 


14 


48 


3.88 


5 


50±1 


42 


4 


33±1 


69 


4 


77 


0.08 




-1.85 + 7.61 


<22.82 


SMM J100026 


7+022230... 


850 


75 


10 


00 


26 


7 


02 


22 


30 


3.83 


3 


08±0 


80 


2 


41±0 


96 


4 


95 


0.09 




1.12 + 4.14 


<12.43 


SMM J100048 


9+023023... 


850 


76 


10 


00 


48 


9 


02 


30 


23 


3.83 


4 


36±1 


14 


3 


42±1 


36 


4 


95 


0.09 


450.253 


16.73 + 5.51 


9.20 + 6.85 


SMM J100005 


5+022958... 


850 


77 


10 


00 


05 


5 


02 


29 


68 


3.83 


3 


58±0 


94 


2 


81±1 


12 


4 


97 


0.09 




0.41 + 4.78 


<14.33 


SMM J100008 


7+021912... 


850 


78 


10 


00 


08 


7 


02 


19 


12 


3.82 


3 


33±0 


87 


2 


60±1 


04 


5 


00 


0.09 




7.53 + 4.52 


<13.55 


SMM J100102 


2+022234... 


850 


79 


10 


01 


02 


2 


02 


22 


34 


3.82 


4 


57±1 


20 


3 


57±1 


43 


5 


04 


0.10 




16.35 + 5.70 


< 17.09 



Table Notes. The >3.6cr 850/Am— detected sources we extract within the central 394'^ of our COSMOS map. The 3.6(T detection threshold is chosen based 
on an estimated 3—5% contamination rate. The list is ordered by detection signal-to-noise ratio (S/N). The "R,AW" flux densities are those measured directly 

from our map. Th e "D EBOOSTED" flux densities are those given after correction for confusion and Eddington boosting as a function of detection S/N, as 
described in section |3.l| We also measure a 90% confidence interval for positional uncertainties and estimated probability of contamination, J^contami from 
the results of our Monte Carlo tests as functions of detection S/N. In the last three columns, we list the corresponding 450/J.m counterparts given in Table [T] 

and fluxes; non-detection deboosted fluxes arc quoted at 3(T limits. 
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ScuBA-2 Galaxies in COSMOS 



Table 6. 850^m-idcntificd >3.6f7 point sources in COSMOS. 



Short Name RAgsQ DECgSQ S/^ SgsQ SgSO -Pcontam ■iZOfim- S450 S450 

Raw Deboosted Sol'rce Raw Deboosted 

[mjy] (mjy] ["] [mjy] [mjy] 



SMM J100014 


1+022704... 


850 


80 


10 


00 


14 


1 


02 


27 


04 


3 


81 


3 


06±0 


80 


2 


39±0 


96 


5 


06 





10 




0.71±4.13 


<12.40 


SMM J095969 


0+022441... 


850 


81 


09 


59 


69 





02 


24 


41 


3 


80 


3 


21±0 


84 


2 


61±1 


01 


5 


10 





10 




-1.34±4.32 


<12.97 


SMM J100054 


8+021945... 


850 


82 


10 


00 


54 


8 


02 


19 


46 


3 


78 


4 


14±1 


09 


3 


22±1 


32 


5 


20 





11 




7.10 + 6.30 


<15.91 


SMM JlOOOlO 


2+021759... 


850 


83 


10 


00 


10 


2 


02 


17 


69 


3 


77 


3 


52±0 


93 


2 


74±1 


12 


5 


22 





11 


460.166 


16.87+4.90 


9.79 + 6.20 


SMM J100041 


8+022111... 


850 


84 


10 


00 


41 


8 


02 


21 


11 


3 


75 


3 


07±0 


82 


2 


39±0 


99 


5 


32 





11 




S.93±4.19 


<12.67 


SMM J100025 


0+022753... 


850 


85 


10 


00 


25 





02 


27 


63 


3 


73 


2 


99±0 


80 


2 


32±0 


97 


5 


40 





12 


460.106 


14.24 + 4.14 


9.60 + 6.16 


SMM JlOOlOl 


3+022439... 


850 


86 


10 


01 


01 


3 


02 


24 


39 


3 


73 


4 


23±1 


13 


3 


2S±1 


37 


5 


41 





12 




-0.42 + 5.41 


<16.22 


SMM J100040 


3+021758... 


850 


87 


10 


00 


40 


3 


02 


17 


68 


3 


70 


3 


49±0 


94 


2 


70±1 


14 


5 


56 





13 




4.66 + 4.74 


<14.22 


SMM J100022 


2+023026... 


850 


88 


10 


00 


22 


2 


02 


30 


26 


3 


70 


3 


31±0 


90 


2 


66±1 


08 


5 


56 





13 


460.179 


14.61 + 4.52 


8.80 + 6.73 


SMM J100020 


2+021727... 


850 


89 


10 


00 


20 


2 


02 


17 


27 


3 


69 


3 


41±0 


92 


2 


63±1 


11 


5 


60 





13 




10.84 + 4.77 


<14.32 


SMM J100032 


8+023049... 


850 


90 


10 


00 


32 


8 


02 


30 


49 


3 


68 


3 


62±0 


96 


2 


71±1 


16 


5 


64 





13 


460.40 


19.16 + 4.81 


14.98 + 6.53 


SMM J100027 


8+022554... 


850 


91 


10 


00 


27 


8 


02 


25 


64 


3 


67 


2 


93±0 


80 


2 


26±0 


96 


5 


66 





13 




6.82 + 4.10 


<12.30 


SMM J100041 


8+022404... 


850 


92 


10 


00 


41 


8 


02 


24 


04 


3 


67 


2 


94±0 


80 


2 


27±0 


97 


5 


66 





13 


460.106 


14.19 + 4.13 


9.67 + 5.15 


SMM J100056 


7+022945... 


850 


93 


10 


00 


66 


7 


02 


29 


46 


3 


67 


4 


80±1 


31 


3 


69±1 


58 


5 


68 





14 




2.09 + 6.13 


<18.3S 


SMM J100002 


6+021632... 


850 


94 


10 


00 


02 


6 


02 


16 


32 


3 


65 


4 


21±1 


15 


3 


23±1 


39 


5 


77 





14 




-1.62 + 6.22 


<18.67 


SMM J095959 


9+022705... 


850 


95 


09 


59 


69 


9 


02 


27 


06 


3 


65 


3 


18±0 


87 


2 


44±1 


05 


5 


77 





14 




6.96 + 4.46 


<13.3S 


SMM J100026 


4+022315... 


850 


96 


10 


00 


26 


4 


02 


23 


16 


3 


64 


2 


93±0 


80 


2 


24±0 


97 


5 


81 





14 


460.54 


16.46 + 4.13 


11.51 + 4.90 


SMM J100013 


4+021807... 


850 


97 


10 


00 


13 


4 


02 


18 


07 


3 


61 


3 


27±0 


90 


2 


60±1 


09 


5 


91 





15 




6.64 + 4.69 


<14.0S 


SMM ,7100014 


5+023008... 


850 


98 


10 


00 


14 


5 


02 


30 


OS 


3 


61 


3 


18±0 


88 


2 


43±1 


07 


5 


96 





16 


460.07 


23.94 + 4.50 


20.65 + 6.19 



Table [5]- Continued. 
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Table 7. Counterpart identifications and Multiwavelength Properties for ScuBA-2 450/xm-detectcd Sources 



NAME Counter- 


RA24 


DEC24 


524 


S24 


RAi 4 


DEC]^ 4 S\ 4 


R-Aopt 


DECopt 


t^opt i 


3.6Mm Zphot 


part # 






["1 






["1 If-iy] 






("] [AB-„ 


lag] [AB-mag] 



450 


00 


1/2 






2/2 


450 


01 


1/1 






2/2 


450 


02 


1/3 






2/3 






3/3 


450 


03 


1/1 


450 


04 


1/3 






2/3 






3/3 


450 


05 


1/1 


450 


06 


1/2 






2/2 


450 


07 


1/1 


450 


08 


1/2 






2/2 


450 


09 


1/1 






2/3 






2/3 


450 


10 


1/1 






2/2 


450 


11 


1/2 






2/2 


450 


12 


1/1 


450 


13 


1/1 


450 


14 


1/2 






2/2 


450 


15 


1/1 






2/2 


450 


16 


1/2 






2/2 


450 


17 


1/5 






2/5 






3/5 






4/5 






5/5 


450 


18 


1/3 






2/3 






3/3 


450 


19 


1/2 






2/2 


450 


20 


1/2 






2/2 


450 


21 


1/1 


450 


22 


1/2 






2/2 


450 


23 


1/1 


450 


24 


1/2 






2/2 


450 


25 


1/1 


450 


26 


1/2 






2/2 


450 


27 


1/1 


450 


28 


1/4 






2/4 






3/4 






4/4 


450 


29 


1/2 






2/2 


450 


30 


1/1 


450 


31 


1/2 






2/2 


450 


32 


1/1 






2/2 


450 


33 


1/1 






2/3 






3/3 


450 


34 


1/2 






2/2 


450 


35 


1/1 


450 


36 


1/2 






2/2 


450 


37 


1/2 






2/2 


450 


38 


1/4 






2/4 






3/4 






4/4 


450 


39 


1/2 






2/2 


450 


40 


1/1 


450 


41 


1/2 






2/2 


450 


42 


1/3 






2/3 






3/3 


450 


43 


1/2 






2/2 


450 


44 


1/3 






2/3 






3/3 


450 


45 


1/1 



10:00:33.35 +02:26:01.6 1.87 297±16 



10:01:09.15 +02:22:55.0 1.70 1492 + 248 



10:00:08.01 +02:26:11.0 1.13 287+15 10:00:08.02 +02:26:12.1 0.83 
10:00:19.76 +02:32:03.9 1.87 189±13 10:00:19.75 +02:32:04.3 1.42 



10:00:23.65 +02:21:55.3 1.54 236±16 10:00:23.67 +02:21:55.3 1.68 

10:00:39.25 +02:22:21.1 0.69 544+17 10:00:39.24 +02:22:21.0 0.77 

10:00:14.25 +02:30:19.0 1.26 467+100 10:00:14.20 +02:30:19.1 1.70 

10:00:16.59 +02:26:38.4 0.72 890±17 10:00:16.57 +02:26:38.4 0.92 

10:00:22.27 +02:23:53.5 1.46 95±18 

10:00:25.45 +02:25:44.2 1.53 116±15 

10:00:01.66 +02:24:27.9 1.47 287+12 10:00:01.66 +02:24:28.2 1.80 

10:00:45.01 +02:19:21.0 2.15 228±20 

10:00:04.36 +02:20:59.6 1.50 402 + 14 10:00:04.37 +02:20:59.7 1.62 

10:00:57.28 +02:20:12.0 2.37 132 + 16 10:00:57.26 +02:20:12.2 2.54 



43±11 
138±15 



56+12 
81±15 



10:00:17.23 +02:21:39.7 0.99 333±21 



10:00:50.23 +02:21:18.1 2.36 1421 + 216 

10:00:28.55 +02:19:28.2 0.49 243±16 10:00:28.57 +02:19:28.0 0.39 



09:59:42.98 +02:21:44.9 0.19 401±13 09:59:42.94 +02:21:45.0 0.59 
10:00:56.78 +02:20:12.9 2.17 90±27 



10:00:48.33 +02:29:26.3 0.63 185 + 14 

10:00:10.27 +02:26:25.4 1.29 143±13 

10:00:21.19 +02:30:55.5 2.32 255+18 

09:59:45.11 +02:22:55.1 2.50 726±14 09:59:45.10 +02:22:55.1 2.67 

10:01:01.24 +02:28:00.6 2.88 104+15 

10:00:10.38 +02:22:24.4 0.98 381±14 10:00:10.36 +02:22:24.8 1.14 



10:00:36.88 +02:19:41.3 2.85 87+15 
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Table Notes. All possible counterparts (24/Am, radio and optical/near-IR) for ScUBA-2 450^im-identified >3.6o" sources. The best-guess counterpart optical 
position is given as the first counterpart and is always the closest source. Counterparts with optical/near-IR positions in bold are submillimeter sources 
whi ch h ave been observed interferon! etrically in the literature, therefore the correct counterpart is known as the source at the given position. See the text 

(§[3751 for details on reliability of counterpart identifications. All positional offsets, 524, '^1.4 and (5opt are measured with respect to the ScUBA-2 450/J.m 
position. [-A-] The z-band and 3.6/J.m magnitudes for 450.03/850.00 is measured from the HST imaging, where it can be decoupled from the bright elliptical 
nearby (1.6"). [%] The 24/jm flux density for 450.07 is originally measured as 1.3 mjy; this is an obvious blend with a nearby bright source. We have 
approximated the contribution from 450.07 to be '-^38%, thus we arrive at 467/Ajy. Sources whose optical magnitudes arc given as " — " are dropouts and the 

position is measured from the Spztzer-IRAC 3.6/J.m. 



© 0000 RAS, MNRAS 000, 000-000 



ScUBA-2 Galaxies in COSMOS 



NAME COIMTEB- RA24 DEC24 ^24 S24 RAj 4 DECi 4 61 4 Si 4 RA„pt DEC„pt «„pt i 3.6f,m j!p|„t 

PART # ["] [MJy] ["] [fiJyl ["] [AB-mag] [AB-mag] 



450 


.46 


1/2 




























10 


:00: 


:32 


46 


+ 02 


21 


:49 


.1 





.74 


24 


.77 


22.16 
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:08 


.3 


2 


.16 


26 


.64 


- 


- 


450 


.53 


1/1 


10:00 


10.33 


+ 02:20:27.1 


1.19 


413±15 


















10 


;00: 


:10 


.32 


+ 02 


;20 


:27 


.4 


1 


.36 


22 


.96 


19.03 


1.01 


450 


.54 


1/2 




























10 


;00: 


:26 


86 


+02 


23 


:18 


.5 





.40 


24 


.82 


22.74 


1.90 






2/2 




























10 


;00: 


:27 


.01 


+ 02 


;23 


:20 


.7 


2 


.73 


27 


.04 


24.62 


2.01 


450 


.55 


1/2 


10:00 


23.96 


+ 02:17:50.1 


2.78 


162 + 17 


















10 


;00: 


:24 


01 


+02 


17 


:50 


.3 


3 


.26 








— 






2/2 




























10 


;00: 


:23 


.81 


+ 02 


;17 


:54 


.3 


2 


.60 


24 


.73 


22.64 


2.60 


450 


.56 


1/4 




























09 


;59: 


:59 


40 


+02 


30 


:02 


.0 


1 


.29 


26 


.49 




2.82 






2/4 




























09 


;59: 


:59 


34 


+02 


29 


:59 


.4 


1 


.54 


26 


.06 


— 


3.18 






3/4 




























09 


;59: 


:59 


.50 


+ 02 


;30 


:00 


.2 


1 


.69 


26 


.82 


— 


1.94 






4/4 




























09 


;59: 


:59 


46 


+02 


29 


:58 


.2 


2 


.69 


26 


.63 


— 


1.69 


450 


.57 


1/1 




























10 


;00: 


:32 


.46 


+ 02 


;18 


:04 


.4 


1 


.79 


28 


.69 


21.85 


3.48 


450 


.58 


1/1 


10:00 


36.05 


+ 02:21:51.1 


2.28 


194+17 


















10 


;00: 


:36 


03 


+02 


21 


:51 


.1 


2 


.53 


24 


.47 


20.79 


— 


450 


.59 


1/3 


10:00 


19.35 


+ 02:20:24.2 


1.08 


197±52 


















10 


;00: 


:19 


.44 


+ 02 


20 


:24 


.5 





.61 


20 


.61 


19.60 


0.47 






2/3 


10:00 


19.61 


+ 02:20:24.3 


3.05 


86±15 




























_ 








_ 






— 




450 


.60 


NA 














































— 








_ 


450 


.61 


1/1 


10:00 


18.76 


+ 02:28:13.5 


0.32 


136±13 


10 


:00: 


:1S 


76 


+ 02:28:13.7 





20 


92 + 11 


10 


:00: 


;18 


81 


+ 02 


28 


:14 


.1 


1 


.11 


25 


.04 


21.88 


4.58 


450 


.62 


1/2 




























10 


;00: 


:47. 


23 


+ 02 


20 


:48. 


.4 


1 


.71 


27 


.78 




1.96 






2/2 




























10 


;00: 


:47, 


.17 


+ 02 


;20 


:51 


.6 


2 


.74 


26 


.76 




1.97 


450 


.63 


1/3 




























10 


;00: 


:19 


.94 


+02 


21 


:29 


.2 


1 


.20 


25 


.09 


- 


2.36 






2/3 




























10 


;00: 


:20 


.07 


+ 02 


;21 


:31 


.8 


2 


.25 


24 


.40 


23.74 


2.84 






3/3 




























10 


;00: 


:19 


87 


+02 


21 


:28 


.1 


2 


.58 


25 


.19 


24.42 


1.56 


450 


.64 


1/1 




























10 


;00: 


:10 


63 


+ 02 


;27 


:09 


.2 


2 


.19 


26 


.70 


— 


— 


450 


.65 


1/2 




























10 


;00: 


:27 


23 


+02 


24 


:49 


.8 


1 


.29 


25 


.42 


— 


0.70 






2/2 




























10 


;00: 


:27 


.16 


+ 02 


;24 


:47. 


.0 


2 


.58 


27 


.67 


— 


— 


450 


.66 


1/4 


10:01 


04.63 


+ 02:26:34.0 


0.47 


586 + 24 


10 


:01 


:04.64 


+ 02:26:34.0 





38 


86 + 11 


10 


:01 


:04 


63 


+02 


26 


:33 


.2 





.61 


23 


.20 


20.64 


_ 






2/4 




























10 


;01 


:04 


.70 


+ 02 


;26 


:34 


.9 


1 


.45 


22 


.69 


20.44 


_ 






3/4 




























10 


:01 


:04 


53 


+02 


26 


:35 


.5 


2 


.60 


25 


.94 


22.14 


_ 






4/4 




























10 


;01 


:04 


.56 


+ 02 


;26 


:31 


.4 


2 


.64 


21 


.86 


20.60 


_ 


450 


.67 


1/2 


10:00 


00.37 


+ 02:29:03.9 


1.43 


538 + 15 


10 


:00:00 


35 


+ 02:29:04.4 


1 


76 


114+14 


10 


;00: 


:00 


.35 


+ 02 


;29 


:04 


.1 


1 


.47 


21 


.37 


19.06 


0.93 






2/2 




























10 


:00: 


:00 


42 


+02 


29 


:00 


.5 


2 


.64 


25 


.83 


23.47 


1.67 


450 


.68 


1/1 




























09 


;59: 


:45 


.92 


+ 02 


;18 


:40 


.9 


4 


.39 


22 


.94 


20.77 


1.12 


450 


.69 


1/2 




























10 


;00: 


:47, 


55 


+02 


25 


:21 


.6 


1 


.16 


26 


.82 


23.67 


2.04 






2/2 


10:00 


47.29 


+ 02:25:20.7 


3.20 


251 + 14 


















10 


;00: 


:47, 


.27 


+ 02 


;25 


:20 


.3 


3 


.45 


22 


.28 


19.36 


0.96 


450 


.70 


1/3 


09:59 


52.66 


+ 02:27:12.8 


1.13 263±15 


















09 


;59: 


:52 


.66 


+02 


27 


:14 


.0 


2 


.30 


23 


.77 


21.06 


3.03 






2/3 




























09 


;59: 


:52 


.59 


+ 02 


;27 


:10 


.3 


1 


.67 


26 


.17 


22.47 


2.89 






3/3 




























09 


;59: 


:52 


.62 


+02 


27 


:09 


.3 


2 


.47 


26 


.91 




2.34 


450 


.71 


1/3 




























09 


;59: 


:46 


.40 


+ 02 


29 


:32 


.0 





.42 


26 


.81 


- 


1.96 






2/3 




























09 


;59: 


:46 


.45 


+02 


29 


:30 


.1 


1 


.84 


26 


.60 


- 


2.12 






3/3 




























09 


;59: 


:46 


.40 


+ 02 


29 


:29 


.0 


2 


.75 


27 


.67 




2.25 


450 


.72 


1/4 




























09 


;59: 


:50 


.36 


+ 02 


;20 


:16 


.5 


2 


.32 


26 


.93 




2.06 






2/4 




























09 


;59: 


:50 


68 


+02 


20 


:16 


.3 


2 


.45 








0.91 






3/4 




























09 


;59: 


:50 


44 


+ 02 


20 


:18 


.9 


2 


.48 


25 


.61 


22.92 


1.43 






4/4 




























09 


;59: 


:50 


54 


+02 


20 


:19 


.6 


2 


.87 


24 


.92 


23.11 


1.72 


450 


.73 


NA 
























































450 


.74 


1/3 




























10 


:00; 


;41 


.39 


+02; 


:16: 


;39 


.9 





.99 


26 


.46 




0.96 






2/3 




























10 


;00: 


:41 


.26 


+02 


:16: 


;42 


.1 


1 


.96 


22 


.79 


21.36 


0.76 






3/3 




























10 


;00: 


:41 


.40 


+02 


16 


:38 


.4 


2 


.36 


26 


.44 




2.16 


450 


.75 


1/2 




























10 


;00: 


:34 


.20 


+ 02 


34 


:22 


.4 


1 


.48 


26 


.81 


24.79 


2.26 






2/2 




























10 


;00: 


:34, 


33 


+02 


34 


:24 


.0 


2 


.32 






23.01 




450 


.76 


1/2 




























10 


;00: 


:21 


76 


+02 


31 


:14 


.0 





.77 


31 


.02 


24.83 


1.96 






2/2 




























10 


;00: 


:21 


.93 


+02 


31 


:15 


.7 


2 


.92 


25 


.27 


23.61 


0.73 


450 


.77 


1/2 




























09 


;69: 


:52 


65 


+02 


22 


:57. 


.8 





.92 


25 


.12 


21.05 


1.43 






2/2 




























09 


:59: 


:52 


62 


+ 02 


23 


:00 


.0 


1 


.41 


25 


.64 


22.35 


1.66 



Table [7] Continued. 
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Table 8. Counterpart identifications and Multiwavelength Properties for marginal ScuBA-2 450/xm- and 850/^m-detected Sources 



NAME 


Counter- 


RA24 


DEC24 


^24 -524 


RA1.4 


DEC]^ 4 4 Si 4 


RAopt 


DECopt i5opt i 3.6^171 Zpi,pt 




part # 






["] [A^.iy] 




["1 [f'Jy] 




["] [AB-mag] [AB-mag] 



450.78 / 850.62 


1/2 




2/2 


450.81 / 850.10 


1/3 




2/3 




3/3 


450.86 / 850.55 


1/2 




2/2 


450.87 / 850.09 


1/2 




2/2 


450.94 / 850.23 


1/3 




2/3 




3/3 


450.95 / 850.103 


1/2 




2/2 


450.96 /850.133 


1/1 


450.99 / 850.33 


1/4 




2/4 




3/4 




4/4 


450.103 / 850.128 


1/2 




2/2 


450.105 / 850.85 


1/1 


450.106 / 850.92 


1/2 




2/2 


450.110 / 850.108 


1/2 




2/2 


450.113 / 850.90 


1/4 




2/4 




3/4 




4/4 


450.126 / 850.159 


1/4 




2/4 




3/4 




4/4 


450.133 / 850.131 


1/4 




2/4 




3/4 




4/4 


450.134 / 850.52 


1/1 




2/4 




3/4 




4/4 


450.135 / 850.163 


1/1 




2/4 




3/4 




4/4 


450.166 / 850.83 


1/2 




2/2 


450.168 / 850.01 


1/4 




2/4 




3/4 




4/4 


450.173 / 850.104 


1/1 




2/3 




3/3 


450.179 / 850.88 


1/4 




2/4 




3/4 




4/4 


450.185 / 850.105 


1/4 




2/4 




3/4 




4/4 


450.189 / 850.114 


1/3 




2/3 




3/3 


450.193 / 850.49 


1/6 




2/6 




3/6 




4/6 




5/6 




6/6 


450.206 / 850.53 


1/4 




2/4 




3/4 




4/4 


450.215 / 850.27 


1/3 




2/3 




3/3 


450.240 / 850.67 


1/5 




2/5 




3/6 




4/5 




5/5 


450.247 / 850.151 


1/5 




2/5 




3/5 




4/6 




5/5 


450.252 / 850.40 


1/6 




2/6 




3/6 




4/6 




5/6 




6/6 


450.253 / 850.76 


1/3 




2/3 




3/3 



10:00:34.33 +02:21:21.4 0.53 130±17 10:00:34.37 +02:21:21.6 0.21 517+23 10 



10:00:25.29 +02:18:45.9 2.76 204 + 29 10:00:25.28 +02:18:46.2 2.53 58 + 12 



10:00:22.12 +02:28:44.8 2.50 
10:00:52.91 +02:19:09.4 2.94 



128 + 15 
109 + 15 



09:59:50.90 +02:28:22.5 2.19 217 + 15 



10:00:25.19 +02:27:55.8 2.58 173 + 13 



10:00:32.06 +02:30:50.7 1.04 
10:00:31.84 +02:30:50.0 3.68 



291 + 60 
86+14 



10:00:29.67 +02:21:29.6 2.74 518 + 16 
10:00:11.87 +02:29:35.9 0.20 110 + 16 

09:59:52.96 +02:26:42.7 1.64 340 + 14 



10:00:05.38 +02:25:16.1 0.67 260 + 20 



10:00:22.77 +02:30:25.3 2.02 103 + 14 



10:00:14.10 +02:28:38.6 1.94 1345+17 10:00:14.07 +02:28:38.8 2.30 



10:00:05.09 +02:17:14.6 4.08 159 + 28 



10:00:23.69 +02:19:15.0 1.69 119 + 17 



10:00:13.58 +02:22:25.6 2.78 605 + 14 10:00:13.55 +02:22:25.7 2.43 88 + 15 



10 


00 


31 


05 


+ 02 


27 


51 


8 





40 


27 


35 


22.40 




_ 


10 


00 


31 


05 


+ 02 


27 


54 


2 


2 


55 


26 


26 


24.53 




_ 


10 


00 


34 


33 


+ 02 


21 


21 








76 


25 


34 


21.78 


1 


75 


10 


00 


34 


43 


+02 


21 


19 


6 


2 


36 


25 


18 


23.21 


1 


66 


10 


00 


34 


53 


+ 02 


21 


20 


2 


2 


94 


25 


61 


23.80 


1 


95 


09 


59 


50 


90 


+02 


27 


44 


6 


1 


34 


25 


86 


23.56 


3 


69 


09 


59 


60 


97 


+02 


27 


46 


7 


2 


07 


26 


80 


24.16 


3 


69 


10 


00 


25 


32 


+ 02 


18 


46 


3 


2 


98 


26 


22 


21.45 


2 


71 


10 


00 


25 


20 


+02 


18 


46 





1 


90 


26 


82 


21.87 


4 


85 


10 


00 


22 


26 


+ 02 


28 


43 


2 





83 


21 


82 


20.20 





71 


10 


00 


22 


04 


+02 


28 


43 


9 


2 


92 


24 


18 


20.44 


3 


34 


10 


00 


22 


32 


+02 


28 


45 


2 


2 


95 


24 


83 


24.95 





62 


10 


00 


62 


86 


+02 


19 


09 


2 


2 


56 






_ 






10 


00 


62 


74 


+ 02 


19 


06 


6 


1 


48 


20 


41 


20.97 





11 


10 


00 


49 


95 


+02 


24 


51 





2 


85 


26 


67 


22.93 


2 


34 


09 


59 


60 


86 


+ 02 


28 


22 


4 


1 


77 


26 


18 


21.21 


2 


16 


09 


59 


60 


69 
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28 


23 


8 


1 


30 


27 


43 
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41 
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23 
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99 


_ 


2 


14 


10 


01 


00 


48 


+02 


19 


45 


6 


1 


94 


27 


51 


_ 


2 


64 


10 


01 


00 


62 


+ 02 


19 


49 


6 


2 


86 


27 


10 


_ 


1 


66 


10 


00 


14 


08 


+02 


28 


38 


7 


2 


16 






18.34 


9 


99 


10 


00 


14 


25 


+ 02 


28 


35 


6 


1 


99 


26 


30 


22.83 






10 


00 


13 


90 


+02 


28 


35 


7 


3 


72 






_ 






10 


00 


05 


13 


+ 02 


17 


19 


4 


1 


14 


24 


62 


21.16 


3 


08 


10 


00 


04 


94 


+02 


17 


19 


3 


2 


03 


26 


59 


22.68 


3 


06 


10 


00 


04 


91 


+ 02 


17 


16 


7 


3 


06 


24 


95 


22.08 


2 


00 


10 


00 


04 


96 


+02 


17 


21 


4 


3 


13 


26 


67 


23.26 


3 


15 


10 


00 


05 


26 


+ 02 


17 


17 


2 


3 


31 


23 


99 


22.50 


3 


09 


10 


00 


05 


06 


+ 02 


17 


14 


1 


4 


63 


22 


61 


20.15 


1 


15 


10 


00 


11 


10 


+02 


16 


06 


4 


1 


47 


25 


86 


24.67 


1 


70 


10 


00 


11 


24 


+ 02 


16 


05 


6 


2 


58 


25 


27 


23.60 


1 


65 


10 


00 


11 


07 


+02 


16 


04 


7 


3 


15 


25 


49 


23.45 


1 


88 


10 


00 


10 


89 


+ 02 


15 


07 


1 


3 


80 


26 


08 


22.93 


1 


44 


10 


00 


23 


67 


+02 


19 


15 


1 


1 


67 


25 


19 


21.73 


2 


76 


10 


00 


23 


67 


+ 02 


19 


19 


6 


2 


77 


28 


16 


_ 


1 


66 


10 


00 


23 


66 


+02 


19 


13 


1 


3 


62 


26 


62 


22.74 


2 


87 


10 


00 


32 


14 


+ 02 


33 


22 


8 


2 


08 












10 


00 


31 


88 


+02 


33 


24 


4 


3 


OS 


27 


33 








10 


00 


32 


03 


+02 


33 


28 


1 


3 


53 


28 


88 


23.08 


2 


73 


10 


00 


31 


88 


+ 02 


33 


22 


6 


3 


69 


27 


48 








10 


00 


32 


23 


+02 


33 


27 


8 


3 


72 


25 


63 


22.48 


2 


64 


10 


00 


41 


36 


+ 02 


26 


34 


4 





95 


26 


95 




1 


74 


10 


00 


41 


42 


+02 


26 


33 


7 


2 


10 


26 


92 


25.04 


1 


78 


10 


00 


41 


48 


+ 02 


25 


35 


2 


2 


80 











92 


10 


00 


41 


46 


+02 


25 


37 


2 


3 


38 












10 


00 


41 


27 


+02 


26 


38 


4 


3 


68 


27 


85 




1 


79 


10 


00 


13 


67 


+02 


22 


25 


4 


2 


52 


23 


10 


19.95 


9 


99 


10 


00 


13 


21 


+ 02 


22 


23 


6 


3 


11 


25 


70 


23.68 






10 


00 


13 


41 


+02 


22 


21 


4 


3 


33 


23 


81 


23.66 






10 


00 


13 


21 


+02 


22 


26 


7 


3 


55 


24 


31 


23.57 






10 


00 


13 


31 


+ 02 


22 


21 


1 


3 


94 


26 


19 


24.18 






10 


00 


13 


34 


+02 


22 


28 


6 


4 


02 


26 


60 








10 


00 


49 


34 


+ 02 


30 


29 


1 


1 


55 


25 


34 


25.04 





25 


10 


00 


49 


64 


+02 


30 


28 


1 


3 


11 


26 


10 


23.20 


4 


OS 


10 


00 


49 


64 


+ 02 


30 


25 


3 


3 


77 


26 


87 


24.70 


3 


13 



Table Notes. All quantities are measured and given as in Table [7] 
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Table 9. Counterpart identifications and Multiwavelength Properties for ScuBA-2 850/xm-detectcd Sources 



NAME 


Counter- 
part # 


RA24 DEC24 524 
["] 


S24 

[MJy] 


RAi 4 DECi 4 5i 4 Si 4 

["] [MJy] 


RA 


opt 


DECopt 




f^opt 

["] 


i 

[AB-mag' 


3.6/^m 
1 [AB-mag] 




hot 


850, 


.04 


1/1 


10:00:15.47 +02:15:50.5 1.78 


123±18 




10:00: 


:15 


.49 


+02 


;15:50, 


.2 


1 


.43 


23.83 


20 


.72 


1 


.39 


850, 


.11 


1/2 








10:00: 


:49, 


.90 


+02 


;22:58, 


,8 


1 


.71 


25.83 


22 


.08 


2, 


.77 






2/2 








10:00: 


:49, 


.94 


+02 


;22:55 


,4 


1 


.82 


27.02 


24 


.07 


2, 


.68 


850, 


.13 


1/3 








10:00; 


:24, 


.12 


+02 


;20:05 


.0 





.28 


- 


24 


.09 


0, 


.86 






2/3 








10:00; 


:24, 


.27 


+02; 


;20:03, 


.8 


2 


.55 


26.14 


24 


.71 


1 


.03 






2/3 








10:00; 


:24, 


.08 


+02; 


;20:02, 


.8 


2 


.58 


24.63 


23 


.81 


0, 


.94 


850, 


.15 


1/1 








10:00; 


:25 


.25 


+02; 


:26:05, 


.9 


1 


.50 


21.70 


20 


.25 






850, 


.16 


1/1 


10:01:03.92 +02:24:48.2 1.08 


229±18 


10:01:03.91 +02:24:48.9 1.69 57+10 


10:01; 


:03 


.92 


+ 02; 


:24:48 


,0 





.89 


24.95 


20 


.82 










2/2 








10:01; 


:03, 


.83 


+ 02; 


;24:45 


.1 


2 


.32 


23.15 


20 


.70 






850, 


.17 


1/1 








10:00; 


:00 


.36 


+ 02; 


;22:56, 


,3 


1 


.26 


26.09 








- 


850, 


.18 


1/2 








10:00; 


:24, 


.02 


+ 02; 


;29:47, 


,3 





.62 


25.67 


23 


.71 


3, 


.01 






2/2 








10:00; 


:24 


.04 


+ 02; 


;29:50, 


,1 


2 


.90 


26.81 


23 


.61 


2, 


.83 


850, 


.19 


1/2 








10:00; 


:15 


.67 


+ 02; 


;24:45 


,5 


1 


.76 


26.87 


24 


.99 


2, 


.24 






2/2 








10:00; 


:15 


.92 


+ 02; 


;24:46, 


.8 


2 


.51 


27.43 


23 


.17 


3, 


.61 


850, 


.20 


1/1 








10:00; 


:27 


.13 


+ 02; 


;31:38, 


.5 


2 


.32 


24.77 


22 


.22 


2, 


.80 


850 


.22 


1/3 


10:00:18.69 +02:16:51.9 3.38 


300+17 




10:00; 


:18, 


.70 


+ 02; 


;16:52, 


.5 


2 


.72 


25.50 


20 


.40 


2, 


.06 






2/3 








10:00; 


:18 


.77 


+ 02; 


;16:58, 


.0 


2 


.83 


27.05 


24 


.11 


2, 


.29 






3/3 








10:00; 


:18 


.98 


+ 02; 


;16:54, 


,7 


3 


.94 


- 










850, 


.25 


1/5 


10:00:12.59 +02:14:44.1 0.20 


331 + 15 




10:00; 


;12 


.58 +02: 


14:44.0 





.20 


24.09 


20 


.31 


2, 


.33 






2/5 








10 


:00; 


:12 


.57 


+ 02 


:14:42 


.7 


1 


.52 


23.67 


20 


.79 


2, 


.78 






3/5 








10 


:00; 


:12, 


.71 


+02 


:14:45 


.6 


2 


.30 


26.28 


22 


.99 


2, 


.61 






4/5 








10 


:00; 


:12, 


.75 


+02 


:14:43, 


.9 


2 


.44 


26.27 


23 


.01 


2, 


.32 






5/5 








10 


:00; 


:12 


.44 


+02 


:14:43, 


.1 


2 


.45 


25.83 


22 


.14 


2, 


.75 


850, 


.29 


1/1 


10:00:50.18 +02:26:17.8 2.57 


155 + 18 




10 


:00; 


:50, 


.15 


+02 


;26:18, 


.3 


3 


.10 


- 


21 


.12 




- 






2/2 








10 


:00; 


:50 


.15 


+02 


:26:14, 


.4 





.77 


25.86 


23 


.42 






850, 


.30 


1/1 


10:00:06.82 +02:22:07.4 1.89 


186+14 




10 


:00; 


:06, 


.80 


+02 


:22:07, 


.5 


1 


.71 


20.94 


20 


.61 





.23 






2/2 








10 


:00; 


:06, 


.86 


+02 


:22:11 


.6 


2 


.64 


25.36 


23 


.80 






850, 


.31 


1/1 


10:00:26.17 +02:17:39.0 2.27 


296+44 




10 


:00; 


:26 


.18 


+02 


;17:39, 


.5 


1 


.89 


24.43 


20 


.46 


1 


.80 


850, 


.35 


1/3 








10 


:01 


:05 


.07 


+02 


:21:51 


.5 


1 


.35 


25.59 


23 


.02 


2, 


.36 






2/3 








10 


:01 


:05 


.02 


+02 


:21:47, 


.9 


2 


.60 


24.73 






1 


.06 






3/3 








10 


:01 


:05 


.15 


+02 


:21:52, 


.7 


2 


.65 


25.08 


22 


.90 


2, 


.63 


850, 


.36 


1/3 


10:00:01.23 +02:17:42.4 3.08 


127+49 




10 


:00 


:01 


.38 


+02 


:17:42 


.4 


2 


.96 


24.20 


20 


.42 


2, 


.05 






2/3 








10 


:00; 


:01 


.21 


+02; 


;17:43, 


,9 


1 


.94 


20.88 


21 


.44 


0, 


.22 






3/3 








10 


:00; 


:01 


.44 


+02; 


;17:45, 


.2 


1 


.97 


24.09 


21 


.70 


2, 


.04 


850, 


.37 


1/3 








10 


:00; 


:35 


.70 


+02; 


;28:27, 


,4 


1 


.84 


28.03 


23 


.31 


4, 


.57 






2/3 








10 


:00; 


:35 


.68 


+ 02; 


;28:28, 


.0 


2 


.08 


27.96 


23 


.31 


3, 


.99 






3/3 








10 


:00; 


:35 


.54 


+ 02; 


;28:29, 


.0 


2 


.92 


28.41 


23 


.98 


3, 


.53 


850, 


.38 


1/4 








10 


:00; 


:23, 


.45 


+ 02; 


;15:35, 


,5 


1 


.39 


24.90 


22 


.75 


2, 


.04 






2/4 








10 


:00; 


:23, 


.58 


+ 02; 


;15:37, 


.5 


1 


.50 


26.39 


23 


.01 


2, 


.72 






3/4 








10 


:00; 


:23 


.63 


+ 02; 


;15:35, 


.0 


1 


.99 


27.04 


24 


.92 


2, 


.58 






4/4 








10 


:00; 


:23 


.36 


+ 02; 


;15:36, 


.4 


2 


.60 


26.60 


23 


.31 


2, 


.43 


850, 


.39 


1/2 








10 


:00; 


:12, 


.21 


+ 02; 


;23:09, 


.8 





.62 


26.48 


24 


.32 










2/2 








10 


:00; 


:12 


.16 


+ 02; 


:23:07, 


.9 


2 


.31 


24.76 


23 


.20 




_ 


850, 


.43 


1/1 








10 


:00; 


:26, 


.31 


+ 02; 


;15:29, 


.0 





.77 


- 






0, 


.84 


850, 


.44 


1/3 


10:00:18.32 +02:22:50.0 1.24 


128+16 




10 


:00; 


:18, 


.29 


+ 02; 


;22:51 


,0 


1 


.81 


25.03 


22 


.78 


2 


.48 






2/3 








10 


:00; 


:18 


.16 


+ 02; 


;22:50, 


,3 


1 


.79 


23.27 


21 


.41 


2, 


.56 






3/3 








10 


:00; 


:18, 


.33 


+ 02; 


;22:47, 


.5 


2 


.06 


27.58 


23 


.56 


2, 


.03 


850, 


.45 


1/1 


10:00:06.75 +02:20:44.7 3.52 


105+33 




10 


:00; 


:06 


.78 


+ 02; 


;20:44, 


.8 


3 


.07 


24.22 


21 


.03 


1 


.32 


850, 


.50 


1/3 






10:00:08.97 +02:20:26.5 3.92 62+12 


10 


:00; 


:09, 


.11 


+ 02; 


;20:26 


.4 


3 


.20 


26.07 


23 


.73 


2, 


.38 






2/3 








10 


:00; 


:09 


.17 


+ 02; 


;20:23, 


,2 





.78 


25.06 


22 


.90 


1 


.97 






3/3 








10 


:00; 


:08 


.96 


+ 02; 


;20:24, 


,7 


2 


.71 


26.43 


23 


.47 


2, 


.87 


850, 


.51 


1/4 


09:59:56.41 +02:18:53.0 1.23 


193+15 




09 


:59; 


:56, 


.39 


+ 02; 


;18:52, 


.7 


1 


.66 


24.39 


20 


.82 


1 


.11 






2/4 








09 


:59; 


:56 


.32 


+ 02; 


;18:56, 


,2 


2 


.64 


25.18 


22 


.04 


1 


.09 






3/4 








09 


:59; 


:56, 


.56 


+ 02; 


;18:56, 


.5 


2 


.89 


21.06 


19 


.07 


0, 


.66 






4/4 








09 


:59; 


:56 


.63 


+ 02; 


;18:53, 


,3 


3 


.04 


25.46 


22 


.37 


1 


.02 


850, 


.54 


1/1 








09 


:59; 


:48 


.89 


+ 02; 


;27:48, 


,9 





.71 


26.82 


23 


.91 


3, 


.07 


850, 


.56 


NA 




































850, 


.57 


1/4 


10:00:06.90 +02:33:08.1 1.20 592+184 10:00:06.82 +02:33:07.9 0.75 82+14 


































2/4 








10 


:00; 


:06 


.81 


+ 02; 


;33:05, 


.4 


1 


.94 


23.93 


22 


.89 


0, 


.50 






3/4 








10 


:00; 


:06 


.68 


+ 02; 


;33:07, 


,2 


2 


.49 


20.13 


18 


.28 


0, 


.74 






4/4 








10 


:00; 


:07 


.01 


+ 02; 


;33:07, 


,2 


2 


.51 


24.78 


21 


.95 


1 


.51 


850, 


.59 


1/4 


10:00:00.80 +02:20:02.4 1.21 


85+13 




10 


:00; 


:00, 


.84 


+ 02; 


;20:02, 


.4 


1 


.60 


26.08 


24 


.06 


3, 


.15 






2/4 








10 


:00; 


:00 


.92 


+ 02; 


;20:01, 


,2 


2 


.13 


26.18 


24 


.01 


2, 


.64 






3/4 








10 


:00; 


:00, 


.77 


+ 02; 


;20:03, 


.4 


2 


.15 


24.22 


21 


.79 


2, 


.97 






4/4 








10 


:00; 


:00 


.60 


+ 02; 


;20:00, 


,8 


2 


.69 


24.59 


22 


.85 


0, 


.97 



Table Notes. All quantities are measured and given as in Table [t] Sources which have 450/im counterparts are excluded from this 
table as their 450/^m positional uncertainty is smaller than their 850/im uncertainty; the former is used to determine their counterparts 
which are given in either Table [7] or Table |8] Source 850.22 is in bold text in this table to indicate it was interferometrically observed by 

| |SmoIcic et al.|2012| l however no source was detected; it is likely that the reason for this is that the submillimeter emission eminates 
from multiple possible counterparts. For 850/im sources, we first match positions to 24/im/radio, and list that counterpart as the first if 
it exists, otherwise we provide the nearest-neighbor optical/near-infrared source. 
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NAME Counter- RA24 DEC24 ^24 £'24 RAi 4 DEC1.4 c5i 4 Si. 4 RAopt DECopt <5opt i 3.6fim Zphot 

PART # ["] [MJy] ["] [MJy] ["] [AB-mag] [AB-mag] 



850.60 


1/2 


10:00:00.86 +02:27:41.2 2.54 


239+14 


10: 


:00 


:00. 


,82 


+ 02: 


:27: 


:41. 


.3 


2.03 


24. 


.38 


19. 


.74 


1 . 


,44 




2/2 






10: 


:00 


:00. 


.81 


+ 02: 


:27: 


:39. 


.6 


1.69 


26. 


.20 


22. 


.12 


2. 


.75 


850.61 


1/1 






10: 


:00 


:01. 


.96 


+ 02: 


:28: 


:20. 


.8 


1.34 


23. 


.88 


22. 


.16 


0. 


,53 


850.63 


1/3 






09: 


:59: 


:53. 


.36 


+ 02: 


:18: 


:48. 


.9 


1.27 


25. 


.40 






2. 


.35 




2/3 






09: 


:59: 


:53. 


.32 


+ 02: 


:18: 


:51. 


.4 


1.43 
















3/3 






09: 


:59 


:53. 


.47 


+ 02: 


:18: 


:50. 


.6 


1.55 


25. 


.87 


23. 


.90 


1 . 


,89 


850.64 


1/3 






10: 


:00: 


:00. 


.72 


+ 02: 


:21: 


:37. 


.5 


1.10 


26. 


.44 


24. 


.01 


3. 


.39 




2/3 






10: 


:00 


:00. 


.66 


+ 02: 


:21: 


:36. 


. 1 


1.19 


24. 


.90 


23. 


.58 


2. 


,53 




3/3 






10: 


:00: 


:00. 


.84 


+ 02: 


:21 : 


:36. 


.6 


3.01 


27. 


.94 


23. 


.18 


4. 


.12 


850.65 


1/3 


10:00:59.10 +02:21:11.4 3.68 


172 + 15 


10: 


:00: 


:59. 


.04 


+ 02: 


:21 : 


ill. 


. 1 


4.09 


23. 


.33 


21. 


.76 


0. 


.48 




2/3 






10: 


:00 


:59. 


.13 


+ 02: 


:21 : 


:09. 


.6 


1.97 


23. 


,02 


21. 


.30 


1 . 


,27 




3/3 






10: 


:00: 


:59. 


.30 


+ 02: 


:21: 


:05. 


.7 


2.74 


29. 


,80 






1 . 


.94 


850.66 


1/5 






10: 


:00: 


:24. 


.46 


+ 02: 


:28: 


:30. 


.0 


1.60 


24. 


,98 


23. 


.96 


0. 


,94 




2/5 






10: 


:00: 


:24. 


.44 


+ 02: 


:28: 


:33. 


.5 


2.36 


25. 


.76 


23. 


.26 


2. 


.47 




3/5 






10: 


:00: 


:24. 


.56 


+ 02: 


:28: 


:30. 


.6 


2.56 


23. 


.58 


23. 


.62 


0. 


.36 




4/5 






10: 


:00: 


:24. 


.23 


+ 02: 


:2S: 


:31. 


.7 


2.54 


26. 


,76 


24. 


.77 


1 . 


,52 




5/5 






10: 


:00: 


:24. 


.32 


+ 02: 


:28: 


:34. 


.2 


3.16 


26. 


.58 


23. 


.31 


2. 


.26 


850.68 


1/5 


09:59:55.30 +02:19:51.1 3.13 


218+14 


09: 


:59: 


:55. 


.26 


+ 02: 


:19: 


:51. 


.2 


3.09 


25. 


,31 


20. 


.47 


1 . 


,95 




2/5 






09: 


:59: 


:55. 


.31 


+ 02: 


:19: 


:53. 


.5 


0.71 


26. 


,22 


21. 


.78 


2. 


.04 




3/5 






09: 


:59: 


:55. 


.31 


+ 02: 


:19: 


:56. 


.0 


1.78 


24. 


.09 


22. 


.27 


0. 


.93 




4/5 






09: 


:59: 


:55. 


.42 


+ 02: 


:19: 


:55. 


.0 


1.94 


27. 


,55 












5/5 






09: 


:59: 


:55. 


.10 


+ 02: 


:19: 


:53. 


.4 


3.14 


27. 


.59 


23. 


.11 


2. 


,30 


850.69 


1/1 


iU:U0:4d.DD +Uz:zC4:oo.b 3. Do 


142+15 


10: 


:00: 


:43. 


.64 


+ 02: 


:28: 


:55. 


.3 


4.15 


25. 


,60 


21. 


.28 


1 . 


,80 


850.70 


1/1 


10:00:20.01 +02:30:19.5 1.00 


119+15 


10: 


:00: 


:20. 


.07 


+ 02: 


:30: 


:19. 


.3 


0.96 


23. 


.46 


20. 


.56 


0. 


.78 




2/2 






10: 


:00: 


:20. 


.24 


+ 02: 


:30: 


:21. 


.5 


3.03 


26. 


.66 


23. 


.75 


2. 


.46 


850.71 


1/2 






09: 


:59: 


:44. 


.05 


+ 02: 


:21: 


:06. 


.3 


1.12 


27. 


,01 






2. 


,17 




2/2 






09: 


:59: 


:43. 


.99 


+ 02: 


:21: 


:03. 


.8 


1.51 


24. 


.55 


24. 


.13 


0. 


.90 


850.72 


1/6 


10:00:13.69 +02:17:34.2 3.06 


150+13 


10: 


:00 


:13. 


.70 


+ 02: 


:17: 


:34. 


.8 


3.62 
















2/6 






10: 


:00 


:13. 


.73 


+ 02: 


:17: 


:2S. 


.9 


2.55 






24. 


.84 


0. 


,93 




3/6 






10: 


:00: 


:13. 


.81 


+ 02: 


:17: 


:32. 


.6 


2.67 


25. 


.64 


23. 


.40 


2. 


.08 




4/6 






10: 


:00 


:13. 


.47 


+ 02: 


:17: 


:32. 


.2 


2.88 


25. 


,48 


24. 


.22 


1. 


,82 




5/6 






10: 


:00: 


:13. 


.54 


+ 02: 


:17: 


:28. 


.9 


2.88 


25. 


.59 






0. 


.69 




6/6 






10: 


:00 


:13. 


.63 


+ 02: 


:17: 


:2S. 


.4 


2.88 


26. 


,51 






2. 


,35 


850.73 


1/5 


10:00:20.52 +02:22:49.2 2.93 


112+26 


10: 


:00 


:20. 


.48 


+ 02: 


:22: 


:4S. 


.6 


3.80 


26. 


,48 


22. 


.48 


3. 


,91 




2/5 






10: 


:00: 


:20. 


.63 


+ 02: 


:22: 


:51. 


.2 


0.38 


25. 


,75 


22. 


.90 


1 . 


,82 




3/5 






10: 


:00 


:20. 


.70 


+ 02: 


:22: 


:51. 


.3 


0.67 


25. 


,07 


23. 


.14 


1. 


,65 




4/5 






10: 


:00: 


:20. 


.69 


+ 02: 


:22: 


:53. 


.2 


1.99 


25. 


.59 


22. 


.43 


2. 


,50 




5/5 






10: 


:00 


:20. 


.69 


+ 02: 


:22: 


:49. 


. 1 


2.16 


25. 


,13 


23. 


.59 


1 . 


,67 


850.74 


1/1 






10: 


:00 


:05. 


.41 


+ 02: 


:14: 


:47. 


.1 


2.78 


27. 


,42 


23. 


.87 


2. 


,76 


850.75 


1/4 






10:00:26.97 +02:22:30.5 3.64 82 + 16 10: 


:00: 


:26. 


.96 


+ 02: 


:22: 


:30. 


.8 


3.40 


25. 


.86 


21. 


.91 








2/4 






10: 


:00 


:26. 


.79 


+ 02: 


:22: 


:31. 


. 1 


1.20 


26. 


,69 


23. 


.28 








3/4 






10: 


:00: 


:26. 


.91 


+ 02: 


:22: 


:29. 


.3 


2.81 


25. 


.99 


23. 


.05 








4/4 






10: 


:00 


:26. 


.52 


+ 02: 


:22: 


:29. 


.4 


3.34 


25. 


,59 


23. 


.43 






850.77 


1/1 


10:00:05.47 +02:29:53.1 5.17 441 + 178 


































2/2 






10: 


:00: 


:05. 


.54 


+ 02: 


:29: 


:49. 


.6 


8.59 






17. 


.36 


0. 


.10 


850.78 


1/6 


10:00:08.77 +02:19:16.2 4.11 


207 + 14 


10: 


:00 


:08. 


.93 


+ 02: 


:19: 


:16. 


.2 


5.07 


22. 


,89 


21. 


.18 


0. 


,94 




2/6 






10: 


:00: 


:08. 


.66 


+ 02: 


:19: 


: 13. 


.2 


1.26 










0. 


.71 




3/6 






10: 


:00 


:08. 


.78 


+ 02: 


:19: 


:13. 


.8 


1.82 


25. 


,98 


24. 


.27 


2. 


,39 




4/6 






10: 


:00 


:08. 


.88 


+ 02: 


:19: 


ill. 


.4 


2.58 


29. 


,41 






3. 


,24 




5/6 






10: 


:00: 


:08. 


.73 


+ 02: 


:19: 


il5. 


.4 


3.19 


26. 


.46 


22. 


.75 


4. 


.35 




6/6 






10: 


:00: 


:08. 


.67 


+ 02: 


:19: 


i09. 


.0 


3.26 


27. 


,65 


23. 


.83 


2. 


.94 


850.79 


1/3 






10: 


:01: 


:02. 


.28 


+ 02: 


:22: 


i34. 


. 1 


0.20 


21 . 


.39 


20. 


.55 








2/3 






10: 


:01: 


:02, 


. 14 


+ 02: 


:22: 


i34. 


. 1 


2.29 


23. 


.34 


22. 


.20 








3/3 






10: 


:01 


:02, 


.27 


+ 02: 


:22: 


i37. 


.0 


2.83 


23. 


,99 


23. 


.40 






850.80 


1/4 






10: 


:00: 


: 14. 


.21 


+ 02: 


:27: 


i05. 


.3 


1.79 
















2/4 






10: 


:00 


: 14. 


.13 


+ 02: 


:27: 


i02. 


. 1 


2.15 


25. 


,23 


22. 


.72 


3. 


,54 




3/4 






10: 


:00: 


: 14. 


.26 


+ 02: 


:27: 


i03. 


.3 


2.28 


24. 


.07 


23. 


.74 


2. 


.40 




4/4 






10: 


:00: 


:13. 


.99 


+ 02: 


:27: 


iOl. 


.7 


3.19 


26. 


.74 


23. 


.94 


3 


.87 


850.81 


1/2 






09: 


:59' 


:59. 


.12 


+ 02: 


:24: 


i41. 


.9 


1.35 


24. 


,11 


23. 


.58 


0. 


,40 




2/2 






09: 


:59: 


:58. 


.97 


+ 02: 


:24: 


i39. 


.8 


1.73 


26. 


.81 






2. 


.81 


850.82 


1/1 


10:00:54.96 +02:19:47.4 2.46 


203+17 


10: 


:00 


:54. 


.95 


+ 02: 


:19: 


i47. 


.1 


2.09 


26. 


,03 


21. 


.11 


2. 


,87 


850.84 


1/1 






10: 


:00: 


:41. 


.67 


+ 02: 


:21: 


il2. 


.7 


2.58 


26. 


.85 


24. 


.03 


2. 


.42 


850.86 


1/1 


10:01:01.30 +02:24:39.4 1.04 


299+19 


10: 


:01: 


:01. 


.27 


+ 02: 


:24: 


i39. 


.9 


1.51 


23. 


.28 


20. 


.03 








2/2 






10: 


:01 


:01. 


.31 


+ 02: 


:24: 


i37. 


.6 


1.73 


21. 


,73 


19. 


.55 






850.87 


1/3 


10:00:40.23 +02:17:55.5 3.18 


284+15 


10: 


:00: 


:40. 


.26 


+ 02: 


:17: 


i55. 


.7 


2.85 


26. 


.99 


20. 


.78 


2. 


.24 




2/3 






10: 


:00: 


:40. 


.33 


+ 02: 


:18: 


iOO. 


.0 


1.75 


26. 


,48 


22. 


.91 


1. 


,45 




3/3 






10: 


:00: 


:40. 


.48 


+ 02: 


:17: 


i56. 


.6 


2.67 


26. 


.55 






1. 


.62 


850.89 


1/4 


10:00:20.26 +02:17:25.1 2.12 


401 + 16 


10: 


:00: 


:20. 


.25 


+ 02: 


:17: 


i25. 


.7 


1.53 


22. 


,53 


19. 


.94 


9. 


.99 




2/4 






10: 


:00 


:20. 


.19 


+ 02: 


:17: 


i27. 


.0 


1.00 


25. 


,42 


22. 


.16 


1. 


,87 




3/4 






10: 


:00: 


:20. 


.24 


+ 02: 


:17: 


i29. 


.0 


1.79 


25. 


.72 


23. 


.02 


1. 


,72 




4/4 






10: 


:00 


:20. 


.37 


+ 02: 


:17: 


i2S. 


.9 


2.43 


24. 


,88 


23. 


.21 


0. 


,92 
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NAME Counter- 
part # 


RA24 DEC24 S24 S24 RAi 4 DECi 4 <5i 4 

["] [fJy] ["1 


Si. 4 


RA 


opt 




DECc 


>pt 




r 

L 


Dpt 

J 


[AB- 


mag 


3 6Mm 


^phot 


850.91 


1/4 






10:00: 




76 


+02 


;25 


:55 






,69 


25 


,95 




2.68 




2/4 








27 




+02 


:25 








tn 












3/4 






in nn 


27 




+02 


:25 


51 








OK 


ar 




27 




4/4 








27 




+02 


:25 


51 


■ 








■ 


2491 


2 2 


850.93 


1/3 


10:00:57.06 +02:29:42.9 4.97 123±10 


10:00: 


57. 


06 


+02: 


:29; 


:42, 


.8 


5 


.08 


24, 


.47 


20.88 






2/3 






10:00: 


56. 


66 


+02: 


:29; 


:44, 


.7 


^ 


.59 


25 


,82 


24.76 






3/3 






10:00: 


56. 


69 


+02; 


:29; 


:48 


.3 


3 


.23 




.25 


20.79 


0.37 


850.94 


1/2 


10:00:02.62 +02:16:34.0 1.84 393±16 10:00:02.61 +02:16:33.8 1.65 


72+11 


10:00: 


:02, 


,61 


+ 02; 


:16; 


:34, 


.2 


2 


,03 


24, 


,62 


20.29 


2.04 




2/2 






10:00: 


:02, 


,62 


+ 02; 


:16; 


:32, 


.1 


0, 


.39 


25, 


,70 


22.38 


4.16 


850.95 


1/3 


09:59:59.82 +02:27:06.6 1.88 387+67 




09:59: 


:59, 


,80 


+ 02; 


:27; 


:07, 


.4 


2 


.74 


22 


,78 


19.60 


1.55 




2/3 






09:59: 


:59, 


,82 


+ 02; 


:27; 


:03, 


.7 


1 


.94 






19.82 


0.74 




3/3 






09:59: 


:59, 


,93 


+ 02; 


:27; 


:02, 


.0 


3, 


.21 






21.05 




850.97 


1/4 


10:00:13.60 +02:18:08.6 2.62 349+15 




10:00: 


13, 


,59 


+ 02; 


:18; 


:0S, 


.5 


2 


.40 


24, 


,02 


21.12 


1.48 




2/4 


10:00:13.57 +02:18:05.0 2.87 148+15 


10:00: 


13, 


,59 


+ 02; 


:18; 


:05, 


.1 


3, 


.00 


20, 


.46 


18.63 


0.73 




3/4 






10:00: 


13, 


,29 


+ 02; 


:18; 


:08, 


.7 


2 


.82 


25, 


.22 


22.91 


1.80 




4/4 






10:00: 


13, 


,49 


+ 02; 


:18; 


:10, 


.1 


2 


,94 


26, 


,71 




0.49 
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Table 10. Physical Characteristics of ScuBA-2450/im-selected Galaxies 



Short 




Estimates from IR.-SED Fitting ( Le Phare Pr 


opcrtics 


Short 




Estimates from IR-SED Fitting {^4.2| 


Le Phare Properties 


Name 




^IR ^PR-IR ^dust ^^dust SFRuv 
(Lq) (MQyr-1) (K) (Mq) (Mgyr-l) 


(Mq) 


Name 




{Lq) (Mq yr-1) (K) (Mq) 


SFRuv 
(Mq yr-1) (Mq) 



450 


00 


2 


86 


(2.1 + l-l> 


X 10I3 


4000 


45.9 + 24.4 S±10® 


50 


2xl0l0 


450 


.56 


2 


60 


(3.2 + 11,3 


)xl0l^ 


500 


= 27.1 - 


60 


Ixlfll" 


450 


01 


2 


88 




XI0I3 


2000 


37.0 + 3.2 


1 + 109 


70 


2 X loll 


450 


.57 


1 


69 


(1.0 + 3-^)xl0l2 


200 


= 27.1 - 


3. 


9x 10* 


450 


02 


2 


16 


(4.5 + f^6 


)xlol2 


800 


= 27.1 




10 


9x lol" 


450 


.58 


3 


48 


(9.6+13/ 


)xioi2 


2000 


45.5 + 8.0 4+10* 


0.0002 


Ixlflll 


450 


03 





34 


(5.6 + ^-«) 


xiolo 


10 


11.2 + 0.56 


9+10* 


0.3 


6x IqI'' 


450.61 


4 


58 


(1 3 + 3'0l 
'1-3-0.6' 


X 10I3 


2000 


56.4+18.6 3 + 10* 


500 


Ixiol'l 


450 


04 


3 


82 




X 10I3 


3000 


39.4+4.5 


i+io" 


300 


7xl0l" 


450.62 


1 


95 


(l-* + 04' 
(2.2 + ^^-*) 


X Ifll^ 


200 


= 27.1 - 


3. 


1x10^ 


450 


05 


3 


99 




X 10I3 


3000 


47.6 + 5.5 


7+10* 


100 


7xl0ll 


450.63 


2 


36 


X 10I2 


400 


= 27.1 - 


3. 


3x10^ 


450 


06 


2 


OS 


— 5.4 


X iqI^ 


1000 


37.0 + 3.4 


5 + 10* 


100 


ixioll 


450.65 





70 


(1.1+^1) 
(,.^+18_^0 


X Iflll 


20 


= 27.1 - 


0.5 


2x10* 


450 


07 


4 


70 


— 2. a 


X 10I3 


7000 


78.6 + 29.6 3+10* 


300 


4xloll 


450.67 





93 


) x loll 


100 


27.6 + 6.4 2 + 10* 


200 


5xl0l0 


450 


09 


1 


90 


(7.3 + 11) 


X iqI^ 


1000 


43.0 + 4.5 


3+10* 


200 


8x10" 


450.68 


1 


12 


(4.9+;y 


)xl0ll 


SO 


= 27.1 - 


10 


2x101^1 


450 


12 





93 




)xl0ll 


60 


= 27.1 




40 


3xl0l0 


450.69 


2 


04 


(i-s+ol) 


X Ifll^ 


300 


= 27.1 - 


5 


3x10^ 


450 


13 


2 


87 


— 0.8 


)xlol2 


1000 


43.4 + 5.5 


4+10* 


100 


2 x loll 


450 


.70 


3 


03 


(8.3+1*/ 


) X 10I2 


1000 


46.9+13.7 3 + 10* 


1000 


Oxlfll" 


450 


14 


1 


47 


X iqI^ 


300 


25.3 + 5.5 


7+10* 


70 


2 X loll 


450 


.71 


1 


95 


,,,+6.7, 

(l-* + 05) 


X 10I2 


300 


= 27.1 - 


5 


2x10^ 


450 


15 


2 


79 


)X1012 


2000 


46.3 + 5.5 


4+10* 


40 


4xl0l" 


450 


.72 


2 


06 


X1012 


300 


= 27.1 - 


s. 


4x10^ 


450 


16 


2 


32 




X iqI^ 


600 


21.1 + 2.9 


5 + 109 


20 


ixioll 


450 


.74 





96 


(2.7+9'*) 


X Iflll 


50 


= 27.1 - 


0.3 


3x10* 


450 


17 


1 


23 




)xl0ll 


100 


22.5+12.2 


6+10* 


0.4 


6x10* 


450 


.75 


2 


26 


(2.6 + 9-^,) 


X Ifll^ 


400 


= 27.1 - 


4. 


2x10^ 


450 


19 


1 


11 


(5.7+f^« 


)xloll 


100 


= 27.1 




4. 


8x10^ 


450 


.76 


1 


96 


(1.4+=-^) 


X Ifll^ 


200 


= 27.1 - 


0.9 


1x10^ 


450 


20 





76 


(2.0 + 3-^) 


X iflii 


30 


15.1 + 5.2 


1 + 10* 


0.9 


8x10* 


450 


.77 


1 


43 


(6.8+// 


) X loll 


100 


= 27.1 - 


20 


sxiflio 


450 


21 





84 


(4.8 + «-=) 


X iqIi 


80 


22.4 + 3.9 


3+10* 


90 


4xl0l" 


450. SI 


1 


75 


(1.5 + f^) 
(6.4+// 


X Ifll^ 


300 


22.9 + 4.5 1 + 10" 


40 


2xl0l" 


450 


22 


2 


11 




)xlol2 


500 


31.5 + 27.0 5 + 10* 


2. 


4x10^ 


450.86 


3 


59 


) X 10I2 


1000 


40.0 + 31.6 5 + 10* 


50 


Ixlfll" 


450 


23 





97 


(3.6 + 11/ 


)xl0ll 


60 


= 27.1 




3. 


8x10* 


450. S7 


2 


71 


(5.9 + */) 


X Ifll^ 


1000 


33.5 + 7.9 9+10* 


100 


Ixioll 


450 


24 





16 


(3.2 + =-^) 


xiolo 


5. 


19.9 + 4.7 


2 + 10'' 


4. 


6XI0IO 


450.94 





71 


(7.9+// 


)xioio 


10 


12.6 + 5.4 9+10* 


0.00004 


2xl0l" 


450 


25 





61 


(i.a+il) 

(3.6 + f/ 


X iqIi 


30 


18.2 + 3.1 


3+10* 


6. 


IXIO" 


450.95 


2 


68 


(3.3+f6) 


X Ifll^ 


600 


34.5 + 10.5 4+10* 


20 


2x101^1 


450 


26 


2 


68 


)X1012 


600 


= 27.1 




80 


7X109 


450.96 


2 


34 


(2.0 + -) 


X Ifll^ 


300 


31.8 + 27.1 4+10* 


10 


2x101^1 


450 


28 


1 


98 


^g.s+l%a 


)xloll 


200 


14.7+2.3 


9+109 


100 


ixiolo 


450.99 


2 


16 


(3.8 + ^1) 


X 10I2 


600 


30.1 + 3.8 9 + 10* 
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Some basic derived properties of 450^m sources using the infrared SED fitting described in § |4.2| The sources' optical/near-infrared 

photometric redshifts are given by Zp. Their integrated 8-1000/im infrared luminosity, Lir, infrared star formation rates SFRi^ 
(calculated via the relation given in |Kennicutt|1998^ , dust temperatures Tdugt, and approximated dust masses A/jugt are calculated 
assuming SEDs with fixed emissivity, /3 = 1.5 and mid-infrared powerlaw slope, a = 2.0. Sources with fewer than three IR photometric 
points have the dust temperature fixed to the mean of the remainder of the sample, measured to be (Tjust) = 32 K. 
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Table 11. Physical Characteristics of ScUBA-2-selected Galaxies 



Short Name Zp Estimates from IR-SED Fitting (^4.2| Le Phare Properties 
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Basic derived properties of 850^m sources (not detected at 450/im) using the infrared SED fitting described in § |4.2[ all values are 

given as in Table [Tol 
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